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The  numerical  results  of  the  second  triple-point  system  show  that  for  a 
given  incident  shock  Mach  number,  the  highest  pressure  is  achieved  through  a 
DMR  instead  of  a  RR..,  The  data  show  that  when  the  transition  of  RR«~»0MR 
occurs,  no  significant  pressure  change  results  in  the  region  behind  the 
reflection  point  P  in  RR,  nor  in  the  region  behind  the  second  triple  point 
T'.  Numerical  simulation  of  pressure  histories  as  functions  of  M  and  0W  at 
the  RR-*-K)MR  transition  boundary  is  required  to  substantiate  tne  present 
analytical^nd  experimental  results. 

An  application  of  reflections  in  pseudo-stationary  flow  to  the 
interaction  of  spherical  blast  waves  with  a  planar  surface  is  shown  and 
discussed.. 
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Summa  ry 


Pseudo-stationary  oblique-shock-wave  reflections  in  sulfur  hexafluoride 
were  investigated  experimentally  and  numerically.  Over  150  experiments  were 
conducted  in  the  UTIAS  10  x  18  cm  Hypervelocity  Shock  Tube  in  the  range  of 
incident  shock  wave  Mach  number  1.25  <  Ms  <  8.0  and  wedge  angle  4°  <  o  <47° 
with  initial  pressures  Pg  ranging  from  4  to  267  torr  (0.53  to  35.60  kPa)  and 
initial  temperatures  \Q  near  300  K.  The  four  major  types  of  shock-wave 
reflection,  ie.,  regular  reflection  (PR),  single-Mach  reflection  (SMR), 
complex-Mach  reflection  (CMR)  and  double-Mach  reflection  (DMR),  were 
observed.  These  were  studied  with  inf i n4 te-fringe  interferograms  using  a 
23-cm  dia  field  of  view  Mach-Zehnder  interferometer .  The 


obtained  and  the  density  distributions  along  the  wedge  surface  are 
for  the  various  reflection  processes.  Four  experimental  results 
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similar  Mach  numbers  are  compared. 
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The  analytical  transition  boundaries  between  the  four  types  of 
shock-wave  reflection  were  established  up  to  Ms  =  10.0  for  frozen  and 
vibrational  equilibrium  sulfur  hexafluoride.  The  examination  of  the 
relaxation  length  under  the  present  experimental  conditions  indicated  that  a 
vi brati ona 1 -equi 1 i bri urn  analysis  is  required.  In  general,  the 
equilibrium-gas  analytical  transition  boundaries  agree  with  the  experiments. 
However,  RR  persists  beyond  the  RR  termination  boundary  determined  by  the 
detachment  criterion,  and  several  Mach  reflections  (MR)  lie  outside  their 
analytically  predicted  domains.  Experiments  substantiate  that  the  CMR  to 
DMR  transition  line  approaches  the  SMR  to  CMR  transition  line  and  merges  at 
a  single  point  on  the  RR  termination  line.  The  analytical  results  are 
compared  with  the  experimental  data  for  the  reflection  angles  6  and  w' ,  and 
the  triple-point  trajectory  angles  x  and  x'>  and  it  is  verified  that  the 
shock-wave  reflection  systems  were  in  vibrational  equilibrium.  It  is  found 
that  negative  values  for  the  angle  w'  are  possible  under  certain  conditions 
of  Mach  reflection.  The  angle  between  the  two  triple-point  trajectories 
(x'-x)  is  also  discussed  and  the  pattern  of  DMR  with  x'  =  0  was  obtained 
experimental ly .  There  is  little  doubt  that  new  criteria  for  the  transition 
lines  are  required  in  order  to  improve  the  agreement  between  analysis  and 
experiment . 
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Notation 


a  speed  of  sound 

Ar  atomic  argon 

C  explosive  charge 

CMR  compl ex-Mach  reflection 

C02  carbon  dioxide 

DMR  double-Mach  reflection 

h  specific  enthalpy,  Planck  constant,  height  of  burst 

He  helium 

I  incident  shock  wave 

K  kink,  degrees  Kelvin,  61 adstone-Dale  constant,  Boltzmann  constant 

L  distance  between  reflection  point  and  wedge  corner, 

distance  between  Mach  stem  and  wedge  corner, 
depth  of  test  section 

L*  horizontal  distance  between  incident  shock  and  wedge  corner 

J l  horizontal  distance  between  kink  and  wedge  corner 

l*  horizontal  distance  between  incident  shock  and  kink 

M,  M‘  Mach  stems 

M..  flow  Mach  number  in  state  (i) 

M$  incident  shock  wave  Mach  number 

M 2 j,-  flow  Mach  number  in  region  (2)  relative  to  kink 

M2p  flow  Mach  number  in  region  (2)  relative  to  reflection  point 

M 2 1*  flow  Mach  number  in  region  (2)  relative  to  triple  point 

MR  Mach  reflection 

NR  no  reflection 

N2  molecular  nitrogen 


molecular  oxygen 
reflection  point 
pressure  in  state  (i) 

ratio  of  pressure  in  state  (i)  to  pressure  in  state  (j) 

universal  gas  constant 

reflected  shock  waves 

regular  reflection 

shock  radius 

si i pstreams 

sulfur  hexafluoride 

single-Mach  reflection 

triple  points 

tanperature  in  state  (i) 

time 

flow  velocity  in  state  (i)  with  respect  to  pseudo-stationary 
coordinates 

flow  velocity  in  region  1  relative  to  laboratory  frame 
velocity  of  incident  shock  wave  relative  to  laboratory  frame 
explosive  charge  weight 

distance  from  reflection  point  or  Mach  stem  along  wedge  surface, 
x-coordinate  in  physical  plane,  ground  range 

height  of  first  triple  point 

y-coordinate  in  physical  plane 

compressibility  factor 

angstrom,  10-10m 

specific  heat  ratio 

difference  between  successive  values 
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angle  between  incident  shock  wave  and  reflected  shock  wave  at  the 
triple  point 

displacement  angle  of  boundary  layer 

0.  angle  of  flow  deflection  in  state  (i) 

0W  actual  wedge  angle 

0'  effective  wedge  angle  {0  +  x) 

w  w 

0"  second  confluence  angle  (0  +  x' ) 

w  w 

02m  maximum  deflection  angle  from  state  (1)  to  state  (2) 

02s  sonic  deflection  angle  from  state  (1)  to  state  (2) 

X.  wavelength  of  light  source 

P  micro 

p-  density  in  state  (i) 

x  relaxation  time 

4^  incident  wave  angle  for  state  (i) 

<i>lK  incident  wave  angle  to  reflected  shock  wave  from  frame  attached  to 

kink  (upstream  side  of  kink) 

<t)jK  incident  wave  angle  to  reflected  shock  wave  from  frame  attached  to 

kink  (downstream  side  of  kink) 

X  first  triple-point  trajectory  angle 

x'  second  triple-point  trajectory  angle 

w'  reflection  angle 

Superscripts 

'  secondary,  effective 

n  normal  component 


initial  or  referenced  value 
thermodynamic  regions 
calculated  value 

according  to  detachment  criterion 
upstream  of  a  shock  wave 
downstream  of  a  shock  wave 
with  respect  to  the  kink 

according  to  mechanical  equilibrium  criterion,  maximum 
no  reflection 

with  respect  to  the  reflection  point 
shock,  according  to  sonic  criterion,  sonic 
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1.0  INTRODUCTION 


Snock  waves  appear  when  particles  in  a  fluid  approach  and  compress  one 
another  with  a  velocity  larger  than  the  local  sound  speed.  Tne  thickness  of 
a  shock  front  is  of  the  order  of  a  few  molecular  mean  free  paths.  The 
occurrence  of  shock  waves  is  found  in  processes  that  yenerate  abrupt 

compression  fronts  such  as  explosions,  electric  discharges  and  supersonic 
flight.  Shock  waves  become  important  in  the  practical  applications  of  fluid 
mechanics  due  to  the  nearly  discontinuous  changes  in  fluid  velocity  and 
pressure  they  create,  and  the  challenging  analysis  of  the  nonlinearity  of 
the  equations  of  motion  [1], 

Current  interest  in  the  problem  of  shock-wave  diffraction  is 

orientated  to  the  investigation  of  the  nuclear  blast  fields  around  aerospace 
vehicles  and  around  flush-mounted  surface  structures  for  their  vulnerability 
and  survivability  studies  [2],  The  valuable  accumulated  data  can  be 
utilized  to  determine  the  nonuni  form  shock  loading  to  be  applied  'n 
structural  analysis  for  the  design  of  aerospace  systems.  When  a  spherical 
blast  wave  interacts  with  a  planar  surface,  the  complete  range  of  shock 
reflections  result.  Regular  reflection  (RR)  occurs  from  0°  incidence  of 
the  blast  wave  with  the  surface,  and  at  higher  incidence  angle  Mach 
reflection  (MR)  occurs  up  to  9U°  incidence  (Fig.  1). 

The  study  of  oblique-shock-wave  reflection  has  been  carried  out  for 

more  than  a  century  by  many  experimentalists  and  theoreticians.  Although 
the  phenomenon  of  regular  reflection  and  sinyle-Mach  reflection  (SMR)  were 

first  observed  by  the  distinguished  physicist,  Ernst  Mach  [3],  back  in  1878, 
it  was  not  until  the  early  1940‘s  that  von  Neumann  |_4,5]  reinitiated  the 
problem  and  explored  it  systematical ly .  Under  the  supervision  of  Bleakney 
[6]  at  Princeton  University,  Smith  |_7  ]  discovered  compl ex-Mach  reflection 
( CMR )  in  1946,  and  White  [8]  identified  double-Mach  reflection  (DMR)  in 
19bl,  while  they  were  investigating  exper  iinenta  1  1  y  nonstat  i  ona  ry 
obi  i  que- shock-wave  reflections.  SMR,  CMR  and  DMR,  as  a  group,  are  called 
Mach  reflection.  The  four  types  of  reflections  are  illustrated 
schematical ly  in  Fig.  2. 

With  these  four  types  of  reflection  being  identified,  it  became 

necessary  to  establish  their  transition  criteria  and  delineate  their  regions 

on  a  transition  map.  The  simplest  laboratory  experiment  designed  for  the 

study  of  shock-wave  diffraction  consists  of  a  two-dimensional  sharp  wedge 
mounted  on  the  wall  of  a  shock  tube.  When  a  planar  shock  wave  collides  with 
a  sharp  compressive  corner  in  a  frozen  gas,  the  type  of  reflection  which 
occurs  as  a  result  of  the  shock-wave  reflection  and  the  deflection  of  the 
nonstationary  flow  behind  it  depends  on  the  incident  shock  Mach  number  Ms 
and  the  corner  wedge  angle  o  For  the  case  of  an  imperfect  gas  where 
nonequilibrium  and  equilibrium  flows  occur,  the  initial  pressure  PQ  and  the 
initial  temperature  T0  are  also  important.  An  analytical  transition 

criterion  for  the  termination  of  RR  in  a  two-dimensional,  pseudo-stationary , 
inviscid  perfect  gas  flow  was  suggested  by  von  Neumann  [4,5]  and  is 
referred  to  as  the  "detachment"  criterion. 

Smith  [7]  found  in  his  experimental  results  that  RR  persisted  beyond 
the  limit  of  detachment  criterion  and  MR  did  not  occur  immediately  when  the 
theoretical  RR  limit  was  exceeded.  Bleakney  and  his  students  [6,8,9] 
further  studied  the  disagreement  found  by  Smith  [7],  but  no  fruitful 
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explanation  was  given.  This  persistence  of  RR  in  the  MR  region  is  called 
the  "von  Neumann  paradox" .  Henderson  and  Lozzi  LlO]  investigated  the 
transition  problem  from  RR  to  MR  experimentally,  both  in  a  wind  tunnel  and 
in  a  shock  tube.  They  proposed  an  alternative  criterion  for  non-stationary 
flow,  which  is  referred  to  as  the  "mechanical  equilibrium"  criterion. 
Instead  of  reducing  the  disagreement  between  theory  and  experiment  found  by 
Smith  L 7 ] ,  tne  mechanical  equilibrium  criterion  made  it  even  worse  [see  11]. 
Hornung  and  Kychakoff  [12]  suggested  another  criterion  called  the  "sonic" 
criterion.  However,  the  transition  boundaries  predicted  by  the  sonic  and 
detachment  criteria  are  almost  identical  and  cannot  be  resolved 
experimental ly . 

When  RR  terminates,  three  types  of  Mach  reflection  can  occur  in 
nonstationary  flows,  and  they  are  SMR ,  CMR  and  DMk.  White  [8]  investigated 
the  reflection  patterns  of  CMR  and  DMR  in  addition  to  SMR,  and  suggested  a 
termination  criterion  for  SMR.  Henderson  and  Lozzi  [1U]  proposed  the 
criteria  for  the  termination  of  SMR  and  CMR.  However,  their  proposal  was 
neither  analytically  nor  experimentally  substant i at ed .  The  criteria  for 
transition  from  SMR  to  CMR  and  from  CMR  to  DMR,  wh’ch  were  initiated  by  Law 
and  Glass  [lb],  were  established  by  Ben-Dor  L 1 1 ]  and  Ben-Dor  and  Glass  [13, 
14].  An  additional  criterion  between  SMR  and  CMR  was  later  appended  by 
Shirouzu  and  Glass  [16], 

When  an  incident  shock  wave  impinges  a  convex  or  a  concave  wedge 
instead  of  a  sharp  wedge,  the  flow  is  truly  nonstationary  since  the  wedge 
angle  varies  along  the  wedge  and  the  entire  flow  process  is  transient. 
Ben-Dor  [11]  and  Ben-Dor  and  Glass  L 1 3 ]  hypothesized  that  in  truly 
nonstationary  flows,  the  transitions  from  RR  to  MR  and  from  MR  to  RR  are 
different.  Ben-Dor  et  al  [17]  and  Itoh  et  al  [18]  investigated  and 
confirmed  that  a  hysteresis  loop  exists  in  the  transition  phenomenon  of  RR 
and  MR  in  truly  nonstationary  flows.  The  existing  criteria  for  the 
transition  of  RR  and  MR,  which  are  for  pseudo-stationary  flows  based  on 
stationary  flow  analyses,  cannot  be  applied  to  the  truly  nonstationary  case. 
The  flow  is  pseudo- stationary  if  it  maintains  a  constant  shock  configuration 
and  grows  linearly  with  time  from  the  instant  when  the  incident  shock  wave 
encounters  the  leading  edge  of  the  wedge.  Jones  et  al  [19]  transformed  the 
equations  of  motion  of  an  unsteady  flow  into  those  of  a  steady  flow.  Since 
the  shock  waves  in  nonstationary  flow  are  usually  curved,  the  steady  state 
analysis  is  applicable  only  in  the  vicinity  of  the  reflection  point  for  RR 
and  the  triple  point  in  MR.  Bazhenova  et  al  [20]  were  the  first  attempting 
to  map  out  the  transition  boundary  lines  experimentally.  However,  only 
small  ranges  of  incident  shock  wave  Mach  number  and  wedge  angle  were 
covered.  Ben-Dor  [11]  was  the  first  who  provided  transition  maps  for 
perfect  and  imperfect  monatomic  and  diatomic  gases  (Ar  and  N2)  of  all  the 
major  types  of  reflection  in  the  range  of  1  <  M$  <  ID. 

Real  gas  effects  in  shock  waves  due  to  the  excitation  of  the  internal 
degrees  of  freedom  (vibration,  dissociation,  electronic  excitation  and 
ionization)  have  been  studied  by  some  researchers.  Mach  reflections  in 
perfect  argon  and  nitrogen,  oxygen  vibrational  ly  excited  air  and 
vibrational ly  excited  carbon  dioxide  were  studied  at  initial  pressures  of  12 
torr  (1.60  kPa)  by  Gvozdeva  et  al  [21].  They  found  that  while  the  incident 
shock  velocity  increases,  the  reflection  angle  decreases,  and  under  certain 
conditions,  becomes  negative.  Naboko  et  al  L22]  investigated  the  flow 
parameters  behind  a  moving  shock  wave  in  relaxing  gases.  Their  results 


completely  excited  for  4.5  <  M$  <  1 1 ;  nitrogen  is  not  vibrational ly  excited 
and  oxygen  molecules  are  in  equilibrium  for  3.7  <  M$<  7.5. 

Hornung  et  al  [23]  studied  the  transition  to  Mach  reflection  in  steady 
and  pseudo-steady  flow  with  and  without  relaxation  for  dissociating  nitrogen 
and  carbon  dioxide.  Besides  real-gas  effects,  they  explained  that  the  von 
Neumann  paradox  is  due  to  the  displacement  effect  of  the  boundary  layer  on 
the  wedge  surface  in  pseudo-stationary  flow.  Ben-Dor  [11]  and  Ben-Dor  and 
Glass  [13,  14]  analyzed  the  domains  and  boundaries  of  nonstationary  oblique- 
shock-wave  reflections  in  dissociating  nitrogen  and  ionizing  argon  under 
equilibrium  conditions.  They  showed  that  real-gas  effects  play  an 
important  role  in  the  size  of  the  regions  and  their  boundaries,  and  the 
analytical  imperfect  boundary  lines  are  different  from  the  perfect  boundary 
lines  depending  on  the  initial  thermodynamic  state  of  the  gas.  A  wide  range 
of  Mach  number  arid  wedge  angle  were  covered  in  thei  r  shock-tube  experiments, 
however  they  were  directed  well  in  the  regions  instead  of  to  the  transition 
boundaries.  From  the  comparison  between  their  experiments  and  calculations 
in  the  transi tion-boundary  maps,  Ben-Dor  and  Glass  erroneously  suggested 
that  nitrogen  should  be  considered  as  an  equilibrium  gas  and  correctly 
treated  argon  as  a  frozen  gas  in  the  range  of  their  experiments.  Ando  [24] 
computed  the  transi tion-boundary  maps  in  carbon  dioxide,  using  five 
different  thermodynamic  state  models  from  dissociational-vibrational 
equilibrium  to  a  perfect  (frozen)  gas.  Ando  concluded  from  his  experimental 
results  that  real-gas  effects  played  no  role  in  the  prediction  of  the 
transition  boundaries,  and  therefore  the  transition  boundaries  are  best 
determined  by  the  perfect  (frozen)  gas  model.  Shirouzu  and  Glass  [16] 
maintained  that,  although  the  perfect  (frozen)  C02  transition-boundary  map 
provides  very  useful  information  for  practical  applications,  the  agreement 
with  Ando's  results  is  accidental  and  the  perfect-gas  model  is  fictitious. 
In  their  studies,  Shirouzu  and  Glass  showed  that  C02  should  be  treated  as  an 
equilibrium  gas.  Lee  and  Glass  [25]  carried  out  an  analytical  and  numerical 
study  of  pseudo-stationary  obi i que- shock-wave  reflection  in  air  and 
established  the  transition  boundaries  up  to  Mach  number  2U  for  both  perfect 
and  imperfect  air.  No  crossings  of  transition  lines  were  found,  and  hence 
they  removed  the  conjecture  of  possible  triple-Mach  reflection. 
Deschainbault  [2b]  investigated  the  transition  boundaries  in  nitrogen  and  air 
for  1.1  <  Ms  <  in.  He  also  studied  the  pressure  histories  on  and  above  a 
specially  constructed  4U°  wedge  to  obtain  more  insight  of  the  complex  flow 
fields.  In  the  range  of  Ben-Dor's  L 1 1  ]  and  Deschambaul t 1 s  [26]  experiment, 
Lee  and  Glass  [2b]  reported  that  air  and  nitrogen  are  frozen,  or  nearly 
frozen,  regarding  vibrational  excitation  and  dissociation,  and  Shirouzu  and 
Glass  [16]  showed  that  for  their  range  of  experiments  and  initial 
conditions,  the  air  and  nitrogen  data  conclusively  support  the  frozen  gas 
model  at  lower  Mach  numbers  (M$  <  6)  and  tend  towards  equilibrium  at  higher 
Mach  numbers  (M$  >6). 

To  substantiate  the  claims  made  by  the  previous  researchers  that 
real-gas  effects  do  play  an  important  role  in  nonstationary  shock  wave 
systems,  an  investigation  of  a  polyatomic  gas  with  many  internal  degrees  of 
freedom  becomes  ideal.  Because  polyatomic  gas  sulphur  hexafluoride  SF^  has 
15  modes  of  vibrational  degrees  of  freedom,  a  low  specific  heat  ratio  of 
1.093  at  room  conditions,  and  the  merit  of  not  being  toxic  with  a  high  index 
of  refraction  useful  for  optical  studies,  it  is  thus  chosen  as  the  test  gas 
for  this  study. 


The  present  work  is  a  continuation  of  the  program  of  studying 
oblique-shock-wave  reflections  at  the  University  of  Toronto,  Institute  for 
Aerospace  Studies  (UTIAS)  by  Molder  [27],  Weynants  [.28],  Law  and  Glass  [15], 
Ben-Dor  [11],  Ben-Dor  and  Glass  [13,14,29],  Ando  [24],  Ando  and  Glass  [3D], 
Lee  and  Glass  [2b],  Shirouzu  and  Glass  [16],  Hu  and  Shirouzu  [31], 

Deschambaul t  [26],  and  Deschambault  and  Glass  [32].  The  objectives  of  the 
present  work  were  as  follows: 

(1)  To  predict  the  domains  and  transition  boundaries  of  the  various 

types  of  reflection  for  frozen  and  vibrational  equilibrium  sulfur 
hexafl uoride . 

(2)  To  investigate  the  various  transition  criteria  by  conducting 

experiments  close  to  the  transition  boundaries. 

(3)  To  study  the  density  distribution  of  the  flow  fields  using 

infinite-fringe  interferometric  technique. 

(4)  To  clarify  the  contribution  of  real-gas  effects  to  the  shock  wave 
system  by  comparing  the  analytical  results  with  the  experimental 
data. 

(5)  To  solve  the  second  triple-point  system  and  determine  the  type  of 
reflection  that  attains  the  highest  pressure. 


In  this  study,  over  150  experiments  were  performed  in  the  10  cm  x  18  cm 
UTIAS  Hypervelocity  Shock  Tube  in  sulphur  hexafluoride  at  initial  pressures 
ranging  from  4  to  267  torr  (0.53  to  35.60  kPa)  and  an  initial  temperature 
near  300  K  in  all  cases.  The  incident  shock  wave  Mach  number  range  was 
1.25  <  M  <  8,  over  a  series  of  sharp  steel  wedges  ranging  from  4°  to  47°. 
Dual  wavelength  laser  interferograms  as  well  as  shadowgrams  were  used  to 
record  the  reflection  process  and  were  obtained  with  a  23-cm  diameter 
f ield-of-view  Mach-Zehnder  interferometer.  The  shock  shapes  and  geometries, 
density  field  (isopycnics)  and  the  density  distribution  along  the  wedge 
surface  and  through  the  bow  shock  wave  in  detached  cases  were  ascertained 
from  the  corresponding  interferograms.  However,  from  the  shadowgrams,  only 
the  shock  shapes  and  geometries  could  be  ascertained.  The  experimental  data 
was  then  compared  with  the  numerical  analyses. 

2 . U  PSEUDO-STATIONARY  OBLIQUE-SHOCK-WAVE  REFLECTION 

2.1  Classification  of  Types  of  Reflection 


When  a  planar  incident  shock  wave  encounters  a  sharp  wedge  in  a  shock 
tube,  it  moves  with  constant  velocity  along  the  wedge  surface.  The  entire 
reflection  problem  can  be  considered  pseudo-stationary  in  a  frame  of 
reference  attached  to  the  confluence  point  of  the  shock  waves.  Hence, 
instead  of  three  independent  variables  x,  y  and  t,  the  phenomenon  becomes 
describable  in  terms  of  x/t  and  y/t,  and  the  flow  problem  is  self-similar 
[19],  The  variable  x  and  y  may  be  measured  relative  to  any  point  moving 
with  constant  velocity  with  respect  to  the  wedge  corner.  The  concepts  and 
definitions  correlated  to  regular  reflection,  single-Mach  reflection, 
complex-Mach  reflection,  and  double-Mach  reflection  are  stated  below  [13]. 
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The  regular  reflection  shock-wave  configuration  as  shown  in  Fig.  3a 
denotes  the  state  ahead  of  and  behind  the  incident  shock  wave  I  by  (0)  and 
(1),  respectively,  and  the  state  behind  the  reflected  shock  wave  H  by  (2). 
The  frame  of  reference  is  attached  to  the  two  shock  confluence  point  P  which 
is  moving  parallel  to  the  surface  of  the  wedge  with  angle  Qw  at  a  constant 
velocity  U0  =  Us  cosec  <t>0  or  a  Mach  number  MQ  =  M$  cosec  <j>0,  where  60  is  the 
incident  angle  between  the  flow  U0  and  the  incident  shock  wave  I  which  is 
given  as 

<J>0  =  90°  -  %  (2.6) 

When  the  flow  UQ  passes  through  the  incident  shock  wave  I,  it  is 
deflected  towards  the  wedge  surface  by  an  angle  Oj  from  its  original 
direction.  The  flow  has  a  new  value  U1  and  the  thermodynamic  properties  are 
changed.  This  supersonic  flow  U,  is  redeflected  by  the  reflected  shock  wave 
R  through  an  angle  02  =  <h.,  so  that  the  flow  U2  is  now  parallel  to  the  wedge 
surface,  thus  satisfying  the  required  boundary  condition.  With  the  frame  of 
reference  attached  to  the  reflection  point  P,  the  flow  configuration  becomes 
stationary  and  the  two-shock  theory  [4,6]  may  be  used  to  find  the  flow 
properties  of  the  regions  around  the  reflection  point  P.  The  method  of 
calculations  is  discussed  in  Appendix  A. 

If  the  flow  U2  behind  the  reflected  shock  R  is  sonic  or  subsonic,  the 
corner  signal  can  reach  the  point  P  and  affect  the  reflected  shock  R. 
Consequently,  R  can  be  curved  near  point  P  (Fig.  3b).  Since  the  velocity 
with  which  the  point  P  moves  along  the  wedge  surface  is  very  high  for  a 
small  incidence  wave  angle  <j>0,  the  flow  U2  is  generally  supersonic  relative 
to  point  P.  Therefore  the  corner  signal  cannot  reach  P  and  there  will  be  a 
region  behind  R  in  which  the  thermodynamic  state  of  the  gas  is  uniform  and 
a  portion  of  the  reflected  shock  R  is  straight  near  P  (Fig.  3c).  The 
reflected  shock  R  can  be  attached  to,  or  detached  from,  the  wedge  corner, 
depending  on  the  flow  Mach  number  M1  behind  the  incident  shock  wave  and  the 
wedge  angle  Qw  [11,  16]  (Fig.  3d  to  f) . 

2.1.2  Mach  Reflection  (MR) 

As  the  wedge  angle  o  is  decreased,  the  flow  behind  the  incident  shock 
cannot  negotiate  the  wedge  surface  through  just  the  reflected  shock  R.  An 
additional  shock  is  necessary  to  alter  the  flow  so  that  it  can  negotiate  the 
wedge  surface  again  by  turning  subsonical ly .  This  additional  shock 
protrudes  forward  into  state  (0)  at  the  lower  part  of  the  incident  shock 
wave  I  and  is  known  as  the  Mach  stem  M.  The  resulting  irregular  reflection 
configuration  is  called  Mach  reflection.  It  is  characterized  by  a 

confluence  of  the  three  shocks  and  a  fourth  discontinuity,  across  which 
entropy,  density,  temperature  and  flow  velocity  are  discontinuous,  but 
pressure  and  the  flow  direction  are  continuous,  and  is  referred  to  as  the 
contact  surface  or  slipstream  S  [33].  The  confluence  of  the  four 
discontinuities  is  called  the  first  triple  point  T. 

Three  types  of  reflection  can  occur  in  nonstationary  flows,  depending 
on  the  flow  Mach  number  behind  the  reflected  shock  wave  R,  namely, 
single-Mach,  com pi  ex-Mach ,  and  double-Mach  reflections.  One  major 

difference  between  these  Mach  reflections  and  regular  reflection  is  that 


they  all  produce  a  subsonic  flow  near  the  wedye  surface  behind  the  Mach 
stem  M,  while  in  RR  the  flow  is  usually  supersonic  near  the  wall  behind  the 
reflection  point  P.  Like  the  reyular  reflection,  the  reflected  shock  R  can 
be  attached  to  or  detached  from  the  wedge  corner. 

2.1.3  Sinyle-Mach  Reflection  (SMR) 

When  the  flow  U2  behind  the  reflected  shock  R  is  subsonic  relative  to 
the  triple  point  T  and  the  reflected  shock  R  is  curved  over  its  entire 
lenyth  up  to  the  triple  point  T,  this  is  identified  as  a  single-Mach 
reflection  (Fiy.  4).  The  frame  of  reference  is  attached  to  the  triple  point 
T  and  it  moves  along  a  straight  line  at  an  angle  x  (first  triple-point 
trajectory)  with  a  constant  velocity  U0  =  IL  cosec  $0,  where  <t>0  is  the 
incident  angle  between  the  incident  shock  wave  I  ana  the  oncoming  flow,  and 
is  given  as 

<t>0  =  90°  "  °‘w  (2.2) 

The  angle  y  is  defined  as  the  angle  between  the  first  triple-point 
trajectory  and  the  wedge  surface,  and  is  referred  to  as  the  effective  wedge 
angle.  They  are  related  by 

X  ■  %  *  X  (2.3) 

The  Mach  stem  M  lies  in  front  of  the  incident  shock  wave  1  and  is  normal  to 
the  wedge  surface,  but  not  necessarily  straight. 

The  flow  in  state  (0)  can  reach  the  region  bounded  by  the  shock  waves  R 
and  M  by  passing  through  two  shock  waves  I  and  R  if  it  is  above  T,  or  only 
one  shock  wave  M  if  it  is  below  T.  From  stability  consideration,  the  gas 
must  be  compressed  to  the  same  pressure,  and  must  move  in  the  same 
direction.  Consequently,  this  gives  rise  to  a  slipstream  S  dividing  the 
two  regions  of  equal  pressure  and  flow  velocities.  The  reflected  shock  wave 
R,  therefore,  needs  to  redeflect  the  flow  so  that  03  =  ©1  -  p>2. 

With  the  frame  of  reference  attached  to  the  triple  point  T,  the  flow 
configuration  becomes  stationary  and  the  three-shock  theory  [4,6]  may  be 
used  to  find  the  flow  properties  around  the  triple  point  T.  The  method  of 
calculations  is  discussed  in  Appendix  A. 

2.1.4  Complex-Mach  Reflection  (CMR) 

When  the  flow  U2  behind  the  reflected  shock  R  is  supersonic  relative  to 
triple  point  T,  a  seyment  of  the  reflected  shock  R  near  T  becomes  straight, 
and  the  flow  behind  it  is  uniform,  whereas  the  rest  of  R  has  continuous 
curvature  until  it  finally  terminates  at  the  surface.  The  resulting  Mach 
reflection  configuration  is  referred  to  as  complex-Mach  reflection  (Fiy.  5) 
which  is  similar  to  single-Mach  reflection  except  that  instead  of  a  curved 
reflected  shock  over  its  entire  length,  here  the  reflected  shock  wave  R 
develops  a  reversal  of  curvature  di scont inuously.  In  some  cases,  a  smooth 
reversal  of  curvature  takes  place.  The  reversal  of  curvature  happens 
because  of  the  existence  of  a  band  of  compression  waves  behind  the  reflected 
shock  wave  [10,11,21]  and  as  a  result,  a  kink  K  is  formed  in  the  reflected 
shock  wave  R.  The  method  of  calculations  used  for  SMR  applies  to  CMR. 


When  the  band  of  compression  waves  at  the  kink  K  of  a  complex-Mach 
reflection  converges  to  form  a  new  shock  wave  M'f  the  kink  becomes  the 

second  triple  point  T'  and  the  flow  Mach  number  MT  Dehind  the  reflected 

shock  wave  R,  in  a  frame  of  reference  attached  to  tne  second  triple  point 
T‘,  exceeds  unity.  Consequently,  two  systems  of  the  tnree-shock 

interactions  exist  in  the  flow,  hence  this  configuration  is  termed 

double-Mach  reflection  (Fig.  6).  In  the  second  tnree-shock  interaction 

system,  the  flows  behind  the  second  triple  point  T‘,  states  (4)  and  (6),  can 
be  obtained  from  state  (1)  by  passing  through  either  the  second  reflected 
shock  wave  R '  or  through  the  reflected  shock  wave  R  and  the  second  Mach  stem 
M'.  A  second  slipstream  S'  is  formed  at  the  second  triple  point  T ' ,  and  the 

second  triple-point  trajectory  angle  y'  is  defined  as  the  angle  between  the 

wedge  surface  and  the  line  extending  from  the  wedge  corner  to  the  triple 

point  T*.  The  method  of  calculations  used  for  SMR  applies  to  the  first 
triple  point  of  JMR.  Since  an  analogy  can  be  made  oetween  the  first  and 
second  triple  points,  the  three-shock  theory  may  also  be  used  to  find  the 

flow  properties  around  the  second  triple  point  T'  [11]. 

2.2  Transition  Criteria 

Witn  the  existence  of  the  four  different  types  of  nonstat i onary 
oblique-shock-wave  reflection,  it  oecomes  necessary  to  determine  the 
transition  criterion  for  one  type  of  reflection  to  terminate  and  another 
type  of  reflection  to  occur.  Establishment  of  tne  transition  criteria 
requires  the  understandi  ng  of  tne  physics  of  the  flow  fields,  and  it  is 
important  in  the  mapping  of  the  transition  boundaries. 

2.2.1  Transition  from  RR  to  MR 

There  are  three  most  quoted  termination  criteria  for  RR  in  the 

literature.  von  Neumann  [_4J  proposed  that  the  transition  from  RR  to  MK 
takes  place  when  the  wedge  angle  R.  is  lowered  to  an  extent  where  the  flow 
deflection  condition  +  ?2  =  d  (see  Sec.  2.1.1)  can  no  longer  be 

satisfied.  This  violation  occurs  when  the  flow  deflection  angle  ^  through 
the  incident  shock  wave  I  exceeds  the  maximum  flow  deflection  angle 
through  the  reflected  shock  R.  Tnis  criterion  is  referred  to  as  the 
detachment  criterion,  and  the  term  detachment  stems  from  steady  flows  where 
the  oblique  shock  wave  detaches  at  this  angle.  The  detachment  criterion  is 
sometimes  called  the  von  Neumann  criterion  and  has  the  analytical  form 


The  criterion  can  be  best  illustrated  by  using  the  pressure-deflection 
(P  -  o)  shock  polars  as  shown  in  Fig.  7,  where  the  I  and  R  polars  represent 
the  incident  and  reflected  shock  waves,  respectively.  Through  an  RR ,  the 
net  deflection  angle  is  zero,  hence  state  (2)  is  at  the  point  where  the  R 
polar  intersects  the  P / P 0  axis,  say  point  A  on  R.  (Fig.  7).  As  the  wedge 
angle  "w  decreases,  the  deflection  angle  \  through  the  incident  shock  wave 
I  increases,  and  the  R  polar  moves  away ’from  the  P/P  axis  until  it  is 
tangent  to  the  P / P Q  axis  (point  6).  decreasing  the  wedge  angale  "w  any 
further,  the  R  polar  no  longer  intersects  the  P/P0  axis  and  RR  is  not 
possible,  ano  MR  occurs.  Therefore,  the  detachment  criterion  corresponds  to 
the  ref  1 ected-shock  polar  R,  where  the  pressure  of  state  (2)  behind  the 


reflected  shock  wave  jumps  from  point  B  (RR  state)  to  point  C  (MR  state). 
Note  that  there  is  a  discontinuous  pressure  jump  during  transition  according 
to  the  detachment  criterion. 

An  alternative  criterion  was  proposed  by  Henderson  and  Lozzi  [10]  based 
on  the  assumption  that  the  system  should  remain  in  mechanical  equilibrium 
during  transition.  They  argued  that  a  system  which  develops  a  pressure 
discontinuity  during  transition  cannot  be  in  mechanical  equilibrium.  In 
order  to  enable  the  system  to  remain  in  mechanical  equilibrium,  the 

transition  between  RR  and  MR  should  occur  at  point  D  (Fig.  7)  where  the  R 
polar  intersects  the  I  polar  (MR  state)  at  the  same  point  where  it  meets  the 

P/P0  axis  (RR  state).  Consequently,  the  mechanical  equilibrium  criterion  is 

illustrated  by  the  refl ected-shock  polar  R  ,  and  analytically  expressed  as 

Oi  +  02  =  03  =  0  (2.5) 

where  ©3  is  the  flow  deflection  through  the  Mach  stem  (Fig.  4) 

Consider  figure  7  now,  according  to  the  mechanical  equilibrium 

criterion,  MR  will  take  place  for  all  reflected  shock  polar  R  above  point 
L),  since  the  transition  criterion  described  by  polar  has  been  exceeded. 
Yet,  according  to  the  detachment  criterion,  RR  will  take  place  for  all  shock 
polar  R  below  point  B,  since  the  transition  criterion  described  by  polar  R 
has  not  been  satistied.  Therefore,  RR  and  MR  are  theoretical ly  possible  for 
all  reflected  shock  polar  R  lying  between  polars  R^  and  Rm.  The 
dual -sol  ution  region  in  the  (M$  -  0^)  plane  is  shown  for  a  perfect  gas  with 
specific  heat  ratio  y  =  1.4  in  Fig.  8,  and  it  can  be  seen  that  the  area  of 
disagreemet  between  these  two  criteria  is  very  large.  However,  experimental 
results  showed  that  the  mechanical  equilibrium  criterion  made  agreement 
between  experiment  and  theory  worse  for  nonstationary  oblique-shock-wave 
reflection.  The  detachment  criterion  is  the  limit  for  the  two-shock  theory 
to  have  solutions.  The  mechanical  equilibrium  criterion  is  the  limit  for 
the  three-shock  theory  to  have  solutions  which  corresponds  to  *  =  0°. 

Hornung  and  Kychakoff  [12]  suggested  another  criterion  in  which  they 
argued  that  in  order  for  a  MR  to  form,  a  length  scale  must  be  available  at 
the  reflection  point;  that  is,  pressure  signals  must  be  communicated  to  the 
reflection  point.  They  argued  that  the  transition  between  RR  and  MR  takes 
place  at  the  sonic  deflection  angle  °2s  ’  the  deflection  angle  at  which  the 
flow  behind  tne  reflected  shock  R  is  just  sonic  relative  to  the  reflection 
point  P.  In  the  (P  -  0)  shock  polar  diagram  of  Fig.  7,  the  transition 
occurs  when  the  refl ected-shock  polar  intersects  the  P/ P 0  axis  at  the  sonic 
point  E.  Consequently,  this  criterion  is  referred  to  as  the  sonic  criterion 
and  is  expressed  by  the  relations 

Ol  +  02s  =  0  ,  M2p  =  1  (2.6) 

where  M2p  is  the  flow  Mach  number  in  state  (2)  in  a  frame  of  reference  fixed 
to  the  reflection  point  P.  The  sonic  criterion  is  illustrated  by  the 
refl ected-shock  polar  R$.  The  state  behind  the  reflected  shock  wave  jumps 
from  point  E  (RR  state)  to  point  F  (MR  state)  during  transition.  However, 
the  sonic  criterion  results  in  a  transition  boundary  which  is  too  close  to 
the  detachment  criterion  so  that  the  differences  are  experimentally 
unresol vabl e . 


Experimental  results  of  nonstationary  oblique  shock-wave  reflection 
show  that  the  detachment  criterion  is  the  correct  one,  although  it  has  taken 
thirty-five  years  to  confirm  this.  Experiments  also  show  that  MR  does  not 
occur  immediately  when  the  theoretical  RR  limit  is  exceeded.  This 
pers  istence  of  RR  in  the  MR  region  has  been  cal  lea  the  von  Neumann  paradox. 
Hornung  et  al  [23],  and  Shirouzu  and  Glass  [16]  discuss  this  paradox  in 
terms  of  boundary-1 ayer  displacement  effects,  and  Shirouzu  and  Glass 
suggests  that  this  can  explain  it.  However,  more  experimental  data  is 
presently  being  obtained  at  UTIAS  to  substantiate  these  findings  [43]. 

2.2.2  Transition  from  SMR  to  CMR 

In  SMR,  the  flow  behind  the  reflected  shock  wave  R  is  subsonic  in  a 

frame  of  reference  attached  to  the  triple  point  T,  and  it  can  negotiate  the 

wedge  surface  subsonical ly .  However,  when  the  incident  shock  wave  Mach 
number  is  increased,  or  in  some  cases  when  the  wedge  angle  is  raised,  the 
flow  Mach  number  behind  the  reflected  shock  wave  R  is  also  increased  and 
eventually  it  exceeds  unity.  Subsequently,  a  SMR  is  no  longer  sufficient, 
because  there  exists  again  a  supersonic  flow  directed  towards  the  wedge 
surface.  Theoretically,  this  supersonic  flow  has  to  be  made  either  parallel 
to  the  wedge  surface  or  subsonic  before  it  reaches  the  wedge  by  means  of  a 
compression  wave  or  a  shock  wave,  so  that  it  can  negotiate  the  wall 
subsonical ly.  The  flow  Mach  number  M2  behind  the  reflected  shock  wave  R 
appears  to  be  reduced  by  passing  through  a  compression  wave  and  a  CMR  is 
formed.  As  the  flow  Mach  number  M2  increases,  the  compression  wave  gets 

stronger  and  finally  this  compression  wave  becomes  a  shock  wave  and  CMR 

terminates.  Henderson  and  Lozzi  [10]  were  the  first  to  suggest  that  a  band 
of  compression  waves  must  exist  in  a  CMR  and  can  converge  to  a  shock  wave. 
Unfortunately,  they  put  forth  neither  analytical  nor  experimental 
substantiation  for  their  suggestion.  Ben-Dor  and  Glass  [13,14]  suggested 
that  the  transition  between  SMR  and  CMR  occurs  when  the  flow  behind  the 
reflected  shock  wave  R  is  sonic  with  respect  to  the  triple  point  T. 
Consequently,  the  criterion  for  the  termination  of  SMR  and  the  formation  of 
CMR  is 

M2T  =  1  (2.7) 

Recently,  Shirouzu  and  Glass  [16]  studied  the  problem  of  whether  a  flow 
is  in  equilibrium  or  frozen  from  a  more  fundamental  assessment  through  the 
appropriate  relaxation  length  as  related  by  the  angle  6  between  the 
incident  and  reflected  waves  (see  Fig.  9a)  rather  than  by  the  (M$  -  Ow) 
plots.  During  the  course  of  their  studies,  they  came  up  with  an  additional 
condition  for  the  exsistence  of  CMR.  They  argued  using  the  results  of  Law 
and  Glass  [16]  that  the  kink  K  moves  with  a  horizontal  flow  velocity  as  for 
state  (1).  This  means  that  the  flow  in  state  (1)  moves  parallel  to  the 
incident  shock  wave  I  and  downward  with  the  frame  of  reference  attached  to 
the  kink  K,  as  shown  in  Fig.  9a.  Since  the  flow  properties  upstream  on 
either  side  of  the  kink  K  are  the  same,  the  exsistence  of  a  band  of 
compression  wave  in  CMR  implies  that  P2  <  p!>  and  |  <t>^  -  90°  |  <  jdw  -  90°  | 
(or  |B|  <  | a | )  as  shown  in  Fig.  9b.  Of  the  four  possible  cases  illustrated 
in  Fig.  9c,  only  configurations  A  and  B  satisfy  the  condition  |p|  <  |a|,  and 
also  the  angle  6  is  larger  than  90°  in  both  cases.  Hence,  6  >  90°  is  a 
necessary  condition  for  the  exsistence  of  a  band  of  compression  waves  which 
is  required  for  the  formation  of  CMR.  This  is  only  a  necessary  condition 
because  configuration  C  also  has  6  >  90°,  although  it  does  not  satisfy  |B| 


<  | a | .  Consequently,  the  transition  from  SMR  to  CMR  takes  place  when  both 
the  former  criterion  Eq  .(2.7)  and  the  new  necessary  condition  are 
satisfied.  The  new  transition  criterion  from  SMR  to  CMR  can  be  expressed 
as 

M2T  >1  and  6  >  90°  (2.8) 

Shirouzu  and  Glass  [16]  showed  that  the  transition  boundaries  based  on  this 
new  transition  criterion  have  improved  the  agreement  with  experiments  in  Ar , 
N2,  air  and  C02. 

2.2.3  Transition  from  CMR  to  DMR 

L  w  and  Glass  [lb]  proposed  that  the  onset  of  the  second  triple  point 
T1  in  l)MR,  which  resulted  from  the  supersonic  flow  in  state  (2)  trying  to 
negotiate  a  compressive  wedge  when  striking  the  wall,  started  when  the  flow 
U2  in  a  laboratory  frame  of  reference  was  supersonic  with  respect  to  the 
wedge  surface.  Gvozdeva  et  al  [34]  suggested  that  DMR  is  formed  as  a  result 
of  the  excitation  of  the  internal  degrees  of  freedom,  which  leads  to  the 
increase  of  the  density  ratio  across  the  incident  shock  front.  Hence,  the 
flow  Mach  number  M2  becomes  supersonic  at  lower  incident  shock  wave  Mach 
number  and  smaller  wedge  angle  for  real  gases  than  for  perfect  gases. 
According  to  Semenov  et  al  [35],  DMR  is  formed  due  to  the  curling  of  the 
primary  slipstream  S  into  a  vortex.  In  Mach  reflection,  the  slipstream  S 
can  either  lie  along  the  wedge  surface  or  curl  into  a  large  vortex  near  the 
wedge  surface  and  the  Mach  stem  M.  However,  Bazhenova  et  al  [20]  indicated 
that  experiments  show  that  OMR  can  occur  when  slipstream  S  does  not  curl 
into  the  vortex  but  extends  over  the  wedge  surface  at  large  wedge  angles. 
Therefore,  the  curling  of  the  slipstream  S  does  not  establish  the  onset  of 
DMR  as  proposed  by  Semenov  et  al  [35].  Bazhenova  et  al  [20]  and  Gvozdeva 
et  al  [21]  used  piezo-gauges  to  measure  pressure  histories  on  the  wedge 
surface  and  found  that  their  results  did  not  agree  with  the  formation 
criterion  of  OMR  suggested  by  Law  and  Glass  [15].  They  considered  that  the 
flow  U2  behind  the  reflected  shock  K  becomes  supersonic  with  respect  to  the 
triple  point  T  as  the  necessary  but  not  sufficient  condition  for  the 
exsi stence  of  DMR.  They  argued  that  when  CMR  takes  place,  there  is  a  small 
secondary  rise  on  the  wall,  and  it  is  more  pronounced  for  DMR  with  the 

second  confluence. 

Ben-Dor  and  Glass  [13,14]  suggested  that  the  transition  from  CMR  to  DMR 
occurs  when  the  flow  U2  behind  the  reflected  shock  R  becomes  supersonic  in  a 
frame  of  reference  attached  to  the  kink  K  (fig.  5).  Their  transition 

criterion  is  expressed  analytically  as 

M2K  =  1  (2.9) 

Comparisons  between  the  transition  boundaries  based  on  this  criterion  and 
experiments  in  nitrogen  [13]  and  air  [26]  in  the  (M^  -  ',')  plane  show  good 
agreement  except  at  lower  M$  and  higher  near  the  RR  boundary  line.  In 
t he  case  of  argon  [14],  the  predicted  transition  line  has  to  be  1  owe r ed  to 
take  into  account  the  DMR  that  penetrated  into  the  CMR  region.  However,  for 
C02  [16, 3U],  the  calculated  boundary  line  fails  to  include  many  CMR  points 

that  lie  in  the  DMR  region,  and  the  boundary  line  must  be  shifted  upward. 

These  discrepancies  can  be  considered  as  inappropriate  for  a  transition 
criterion.  The  criterion  M  ^  -  1  is  only  one  of  the  necessary  conditions 


and  the  exsistence  of  other  necessary  conditions  may  improve  the  agreement. 
3.0  EXPERIMENTAL  INVESTIGATION 


3.1  Experimental  Instrumentation 


The  UTIAS  10  cm  x  18  cm  Hypervelocity  Shock  Tube  [36]  was  used  to 
perform  the  experiments  for  the  present  work.  Bristow  [37]  reported  on 
further  modifications  and  improvements.  An  updated  and  detailed 

documentation  of  the  shock  tube  and  other  associated  facilities  as  well  as 
their  maintenance  and  operation  procedures  can  be  found  in  Ben-Dor  and 
Whitten  L38J.  Therefore,  only  a  brief  discussion  of  the  facilities  relevant 
to  the  present  study  will  be  given. 

Cold-gas  runs  were  used  to  generate  the  range  of  incident  shock  wave 
Mach  numbers  required  for  the  present  study.  The  diaphragm  used  for  these 
cold-gas  driven  shock  runs  consisted  of  several  different  thickness  layers 
of  Mylar  polyester  film.  Initially,  a  diaphragm  of  the  calculated  thickness 
to  give  the  required  pressure  ratio  P30  across  the  diaphragm  was  inserted  to 
separate  the  driver  from  the  channel .  The  channel  was  then  evacuated  down 
to  a  vacuum  of  5  to  10  millitorr  (0.66  to  1.33  Pa)  by  the  mechanical,  roots 
and  diffusion  pumps,  which  were  connected  in  series,  and  the  vacuum  level 
was  measured  by  a  Pi  rani  vacuum  gauge  (Type  GP14U).  Concurrently,  the 
driver  section  was  evacuated  by  a  vacuum  pump. 

The  test  gas  used  for  the  present  work  was  commercial  grade  SF,  with  a 
purity  rating  exceeding  99.8%.  Relatively  high  initial  pressures,  4  to  267 
torr  (0.53  to  35.60  kPa),  were  employed  in  the  experiments  to  meet  the 
triggering  threshold  of  the  shock  wave  velocity  measuring  pressure 

transducers,  to  intensify  the  density  contours  for  the  case  of 

infinite-fringe  interferometry ,  as  well  as  to  preserve  the  reproducibility 

of  shocks.  Carbon  dioxide  and  helium  were  employed  as  the  driver  gases  to 

generate  incident  shock  wave  Mach  numbers  M$  in  the  range  1.25  <  M$  <  8.0  in 


The  initial  pressure  PQ  in  the  channel  section  was  monitored  by  Wallace 
&  Tiernan  type  FA  160  dial  pressure  gauges  (0  -  50  torr,  0  -  200  torr,  0  - 
40U  torr,  and  400  -  800  torr).  The  accuracy  of  measurement  using  these 
gauges  was  listed  below: 


Gauge 

Overal 1 

Error 

Absolute 

Error 

w 

& 

T 

0  - 

50 

±0.1 

torr 

(13  Pa) 

0.2 

torr 

(26  Pa) 

w 

& 

T 

0  - 

200 

±0.5 

torr 

(66  Pa) 

1.0 

torr 

(133  Pa) 

w 

& 

T 

0  - 

400 

+  1.0 

torr 

(133  Pa) 

2.0 

torr 

(266  Pa) 

w 

A 

T 

400  - 

800 

+  1.0 

torr 

(133  Pa) 

2.0 

torr 

(266  Pa) 

The  pressure  of  the  driver  gas  P3  admitted  into  the  driver  section  was 
monitored  by  Heise  type  H1/507C  dual  pressure  gauge  (0  -  500  psi), 

calibrated  in  0.5  psi  (3.45  kPa)  intervals.  The  bursting  driver  pressure  P3 
was  interpolated  to  the  nearest  0.25  psi  (1.72  kPa)  with  an  estimated 
overall  error  of  ±0.25  psi  (1.72  kPa),  and  an  absolute  error  of  0.5  psi 
(3.45  kPa). 
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The  leak  rate  of  room  air  into  the  shock  tube  was  found  to  be  6.1  torr 
(U.81  kPa)  in  a  typical  day.  Since  most  of  the  experiments  were  conducted 
at  relatively  high  initial  pressures,  P  >  4  torr  (0.53  kPa),  and  the  period 
between  the  admission  of  the  test  gas  to  the  bursting  of  the  diaphragm  was 
about  15  minutes.  Ihe  worst  possible  increase  of  the  initial  pressure  due 
to  leakage  would  be  0.08  torr  (10.6/  Pa),  or  no  more  than  2%  of  the  measured 
initial  pressure. 


The  initial  temperature  of  the  test  gas  Tn  was  measured  by  a  mercury 
bulb  thermometer  calibrated  in  0.1°C  interval  s  °which  was  embedded  in  the 
midsection  of  the  channel.  The  temperature  of  the  test  gas  was  assumed  to 
reach  a  thermal  equilibrium  state  in  the  shock  tube  within  5  minutes  after 
the  admission  of  the  test  gas.  The  initial  temperature  was  interpolated  to 
the  nearest  0.0buC  with  an  estimated  overall  error  of  ±0.05°C,  and  an 
absolute  error  of  0.1UC. 


The  shock  speed  was  calculated  from  the  measured  traverse  time  interval 
between  several  stations  in  the  channel  (Fig.  10).  Three  Hewlett-Packard 
(Type  3734A)  and  two  Racal  (Type  SA.45)  digital  counters  were  connected  to  a 
common,  external,  1  MHz  oscillator  providing  all  counters  a  uniform  time 
base  with  a  one  microsecond  resolution.  When  the  shock  wave  arrived  at 
station  D  (5./  m  upstream  of  the  test  section),  all  five  counters  were 
simultaneously  triggered.  tach  counter  was  then  stopped  in  succession  as 
the  shock  wave  travelled  past  subsequent  pressure  transducers  at  stations  F, 
G,  H,  I  and  J.  Atlantic  Research  LD  -  25  piezoelectric  pressure  transducer 
flush  mounted  in  the  shock  tube  were  used  to  detect  the  arrival  of  the  shock 
wave.  The  absolute  error  in  calculating  the  incident  shock  wave  Mach  number 
is  outlined  by  Ben-Dor  and  Whitten  L 38 J .  For  SF6,  the  absolute  error  is 
0.01  at  M$  =  1.25  while  the  absolute  error  is  0.11  at  Ms  =  8.0. 


To  study  the  reflection  process,  a  23-cm  dia  field  of  view  Mach-Zehnder 
interferometer  was  used.  Details  of  the  design  and  operation  of  this 
i nterferometer  were  given  by  Hall  L 3 9 J .  A  detailed  description  of  the 
spatial  resolution  and  alignment  of  the  interferometer  can  be  found  in 
Ben-Dor  and  Whitten  L38J.  The  adjustments  for  the  infinite-fringe 
interferograms  were  very  difficult  in  order  to  obtain  just  one  fringe  over 
the  entire  no-flow  field  of  view,  fortunately,  only  a  few  residual  fringes 
remained  far  from  the  wedge  surface,  and  their  effects  on  the  quality  of 
the  infinite-fringe  interferograms  were  negligible.  The  shapes  of  the 
density  contours  were  initially  set  almost  identical  before  each  run,  and 
all  the  photographs  were  taken  under  that  condition.  The  infinite-fringe 
interferometric  technique  not  only  provides  excellent  isopycnic  contours,  it 
also  is  a  very  sensitive  method  of  measuring  density  values.  The 
shadowgrams  taken  in  the  present  study  used  the  same  equipment  with  the 
light  path  of  the  lower  arm  blocked  by  a  cardboard  which  was  placed  in  front 
of  a  mirror. 


Ihe  interferometric  light  source  used  was  a  giant-pulse  ruby  laser 
consisting  of  a  1RG  model  104A  laser  head  and  power  supply,  and  an 
integrated  TRG  model  2113-1  harmonic  generator  and  TRG  Pockels  Q-switch.  It 
was  found  exper imental  ly  that  the  Pockels  cell  shutter  opened  at  a  time  of 
9U0  microseconds  into  the  flash  lamp  pumping  cycle  gave  an  output  pulse 
width  of  3U  nanoseconds,  which  was  sufficiently  fast  to  freeze  the  shock 
motion  during  the  recording  of  the  interferogram.  The  pulsed  laser  provided 


a  monochromatic  light  at  a  wavelength  of  6943.0  A  (first  harmonic) 
3471.5  A  (second  harmonic). 


anc 


Both  the  interferograins  and  shadowgrams  were  recorded  on  Kodak  Koyal-X 
pan  (1250  ASA)  10  cm  x  13  cm  sheet  films.  The  film  was  developed  in 
dilution  A  of  HC -110  developer  for  9  minutes  for  the  first  harmonic 
i nterferograms  and  shadowgrams,  and  12  minutes  for  the  second  harmonic 
interferograins .  The  control  of  the  laser  light  source  operation  is  shown 
in  Fig.  10.  The  laser  flashlamp  capacitor  bank  was  triggered  by  the  shock 
arrival  at  station  F,  correctly  delayed  by  a  Tektronic  type  555  oscilloscope 
to  get  900  microseconds  of  energy  pumping  by  the  time  the  Pockels  cell 
shutter  opened.  The  Pockels  cell  0-switch  was  triggered  from  station  1,  and 
delayed  suitably  in  a  delay  unit  to  take  the  interferoy ram  of  the  shock  wave 
at  a  desired  location  with  respect  to  station  1. 

In  the  present  work,  wedges  with  angle  of  4°,  10°,  2U°,  30°,  37°, 
40°,  42°,  45°  and  47°  were  used,  with  a  tolerance  of  +10  minutes  of  arc. 
The  wedges  were  constructed  of  mild  steel,  bolted  to  the  bottom  of  and 

mounted  flush  with  the  inside  shock  tube  walls  at  the  test  section  (Fig. 
11a).  This  arrangement  was  adopted  due  to  the  simplicity  in  desiyn  and 

stability  during  impact  of  a  strong  shock  wave.  However,  it  has  a  drawback 
of  producing  a  boundary  layer  interactions  at  the  wedge  corner  which  is  not 
of  important  in  the  present  work.  Gvozdeva  et  al  L4U]  used  the  elevated 
wedge  model  in  Fig.  lib  to  eliminate  the  boundary-1 ayer  interaction  at  the 
corner,  but  a  more  serious  yasdynainic  problem  arose.  In  the  case  of  a 

detached  bow  shock  at  the  wedge  corner,  the  flow  field  beneath  the  wedge 
will  influence  that  above,  thus  disturbing  the  region  of  interest.  A 

symmetric  wedge  (Fig.  He)  would  be  freed  of  the  above  problems,  but  it 
occupies  a  larger  volume,  consequently  generating  greater  disturbances  that 
could  choke  the  flow  and  even  raise  the  temperature  high  enough  to  cause 
damaye  to  the  window  surfaces.  An  asymmetrical  wedye  (Fig.  lid)  with  two 
wedye  angles  owi  and  0W2  has  the  merit  of  producing  two  sets  of  results  for 

each  run,  however,  it  has  all  the  d  i  sadvantages  of  both  models  b  and  c. 

3 . 2  Evaluation  of  Interferograins 

Usiny  inf i n i te- fringe  interferometr ic  technique,  very  small  density 
changes  can  be  recorded  on  an  i nterferoyram .  The  density  difference  Ap 

between  any  two  consecutive  isopycnics  is  determined  by  the  relation  [3D] 


where  \  is  the  wavelength  of  the  light  source,  K  is  the  Gl  adstone-Da  1  e 
constant,  and  L  is  the  depth  of  the  test  section  of  the  shock  tube.  The 
values  of  K  for  several  gases  listed  in  reference  [26]  is  reproduced  in 
Table  1. 

By  using  a  light  source  of  shorter  wavelength,  the  density  contours 
become  more  numerous,  and  a  more  sensitive  density  distribution  is  achieved. 
Compare  Fig.  12a  with  Fig.  12b,  which  were  obtained  simultaneously  in  the 
same  experiment.  The  wavelenyth  of  the  1 iyht  source  fur  Fig.  12a  is  6943  A 
(first  harmonic)  while  the  wavelength  of  the  light  source  for  Fig.  12b  is 
3471.5  A  (second  harmonic).  Although  there  are  twice  as  many  isopycnics 


v:s:ble  with  the  second  harmonic  than  with  the  first  harmonic,  this 

technique  suffers  from  the  limitations  of  the  optics  and  photographic 
emulsions  used,  hence  large  losses  in  -intensity  and  improper  exposure  of 
the  interf erograms .  An  alternative  to  ach4eve  more  sensitive  measurements 
of  the  density  contours  is  to  'increase  the  initial  pressure  of  the  test  gas, 
for  a  fixed  Mach  number.  There  will  be  larger  relative  pressure  changes, 

and  larger  density  changes,  thus  g’ving  r-ise  to  more  detailed  and  accurate 

isopycnics.  Deschambault  [26]  reported  a  case  where  a  CMR  observed  at  lower 
initial  pressure  needed  to  be  reclass4 f ied  as  a  DMK  when  the  same  experiment 
was  repeated  with  the  initial  pressure  threefold  h4gher.  Th's  method  also 
has  its  practical  limits,  as  the  incident  Mach  number  is  governed  by  the 

required  driver  pressure  and  burning  of  the  window  may  occur  at  h4gh  Mach 
numbers  and  h4gh  wedge  angles.  Experimentalists  should  be  aware  of  the 
limitations  of  the  equipment  and  how  the  initial  conditions  can  change  the 
state  of  a  real  gas  and  affect  the  classification  of  a  reflection. 
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ve  is  regarded  as  stable  if  the  4 nterf erogram 
unstable  if  it  is  non-planar.  In  the  present 
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Since  it  is  not  possible  to  analyze  and  discuss  each  i  nterf  erogram 
obtained  in  the  present  study  in  detail,  only  several  interesting  cases  will 
be  presented.  Table  2  is  a  record  of  all  the  experiments  performed 
including  the  initial  conditions  and  important  measurements  of  each 
experiment.  In  the  following  sections,  examples  of  attached  and  detached 
RR,  SMR,  CMR  and  DMR  are  presented,  including  their  interferometri c  data. 


3.2.1  Regular  Reflection 


Case  1 


=  1.82,  0W=  47°, 


PQ  =  149  torn  (19.87  kPa),  TQ  =  23.1°C 


An  4 nf i ni te-f ri nge  i nterferogram  of  a  detached  RR  is  shown  in  Fig.  13. 
In  the  same  figure,  the  evaluation  of  the  density  field  is  also  presented. 
Note  how  clearly  the  isopycnics  appear  in  the  i nterf erogram .  The  density 
contours  are  enhanced  due  to  the  relatively  high  initial  pressure  of  the 
test  gas  used  for  the  experiment.  There  exists  a  small  uniform  superson4c 
region  4mmediately  behind  the  reflection  point  and  the  reflected  shock  wave 
is  stra4ght  just  in  this  region.  As  one  moves  away  from  the  reflection 
po4nt  toward  the  wedge  corner,  the  flow  decreases  in  strength  and  becomes 
subsonic.  Since  the  corner  s4gnal$  are  able  to  influence  the  flow  up  to  the 
reg4on  where  it  just  becomes  superson4c,  the  reflected  shock  wave  in  the 
subsonic  region  becomes  curved  throughout  the  rest  of  its  length  and 
terminates  on  the  bottom  wall  of  the  shock  tube  (detached). 


In  Fig.  14,  the  wall-density  distribution  p/pQ  plotted  against  the 
normalized  distance  L  is  shown  for  this  example.  For  regular  reflection,  L 
is  defined  to  be  the  distance  from  the  reflection  point  P  to  the  wedge 
corner  (see  Fig.  16a).  For  Mach  reflection,  L  is  defined  as  the  distance 


from  the  foot  of  the  Mach  stem  M  to  the  wedge  corner  (Fig.  15b).  The  points 

on  the  wall-density  distribution  plot  were  measured  from  the 

i  nterferograms .  The  wall-density  plot  also  indicated  a  small  constant 
region  behind  the  reflection  point  before  the  density  drops  continuously. 
Near  the  corner  of  the  wedge  surface,  the  density  increases  slightly  and 

then  drops  sharply  but  continuously  behind  the  corner.  The  sharp  decrease 
in  density  is  seen  in  the  interferogram  as  a  small  region  at  the  corner 
where  the  isopycnics  converge,  which  may  be  caused  by  the  interaction  of  the 
shock  wave  with  the  boundary  layer.  If  the  interaction  were  stronger,  a 
lambda-shock  configuration  as  shown  later  in  case  4  would  occur.  Finally, 
after  another  drop  in  density  due  to  the  reflected  shock,  the  density  in 
state  (1)  is  reached. 

Case  2  :  Ms  =  4.57,  Ow  =  40°,  PQ  =  7.5  torr  (1  kPa),  TQ  =  24.2°C 

An  example  of  an  attached  RR  is  shown  in  the  infinite-fringe 
i  nterf erogram  of  Fig.  16.  The  corresponding  wall-density  distribution  plot 
is  given  in  Fig.  17.  The  incident  shock  wave  Mach  number  for  this  case  is 
approximate ly  2.5  times  higher  than  that  in  case  1.  The  flow  in  region  (1) 

is  strong,  thus  forcing  the  reflected  shock  wave  close  to  the  wedge  surface. 
The  region  behind  state  (2)  is  uniform,  hence  the  flow  behind  the  reflected 
shock  is  subsonic  and  there  is  communication  of  the  corner  signal 
everywhere.  Note  that  in  the  wall-density  distribution  plot,  the  density 
near  the  reflection  point  is  about  the  same  as  the  density  at  the  corner, 
and  a  density  drop  i n  between  of  less  than  10%.  A  weak  disturbance  ahead 
of  the  reflected  shock  wave  which  had  been  observed  previously  by  Ando  [24] 
in  CO.  and  Deschambault  [26]  in  air  is  also  seen.  This  phenomena  occurs  at 
high  Mach  number  and  Deschambault  [26]  tried  to  explain  that  the  shock 
bifurcation  seen  at  the  wedge  corner  continues  along  the  window  boundary 
layer,  gradually  diminishing  in  strength  and  interaction. 

3.2.2  Single-Mach  Reflection 

Case  3  :  Ms  =  1.54,  ow  =  42°,  PQ  =  138  torr  (18.40  kPa),  TQ  =  23.9°C 

A  detached  SMR  i nf i nte-fri nge  interferogram  with  the  evaluation  of  the 
density  field  is  shown  in  Fig.  18.  The  well  defined  and  numerous  isopycnic 
lines  in  the  i nterferogram  are  the  result  of  using  relatively  high  initial 
pressure.  No  uniform  region  is  observed  behind  the  triple  point  and  the 
reflected  shock  since  the  flow  is  subsonic  in  state  (2).  The  corner  signal 
can  propagate  upstream  to  the  triple  point  and  influence  the  reflected  shock 
wave.  As  a  result,  the  reflected  shock  shape  is  curved  throughout  its 
entire  length  from  the  triple  point  to  the  bottom  of  the  shock  tube,  where 
it  terminates  as  a  detached  shock.  The  straight  slipstream  clearly  emanates 
from  the  triple  point  and  merges  with  the  boundary  layer  on  the  wedge 
surface.  The  wall-density  distribution  is  plotted  in  Fig.  19.  The  density 
decreases  gradually  behind  the  Mach  stem  and  then  increases  slightly  at  the 
corner.  As  one  moves  a!ong  the  bottom  wall  of  the  shock  wave  away  from  the 
wedge  corner,  the  density  is  seen  to  be  dropping  again,  and  the  sharp 
discontinuous  drop  in  density  across  the  reflected  shock  to  state  (1)  is 
clearly  shown. 
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Case  4  :  Ms  =  2.11,  0W  =  20°,  PQ  =  29.4  torr  (3.92  kPa),  TQ  =  23.8°C 

In  Fig.  20,  an  example  of  an  attached  SMR  is  shown.  Compare  the 
interferogram  of  this  case  to  the  previous  case  (Fig.  18).  The  isopycnics 
shapes  are  still  clearly  shown  in  the  interferogram,  but  note  that  in  this 
case  there  are  fewer  isopycnics  recorded.  This  is  due  to  the  relatively 
lower  initial  pressure  used  and  a  smaller  compressive  wedge  angle,  although 
the  incident  Mach  number  is  slightly  higher  in  the  present  case.  Since  the 
triple  point  trajectory  angle  x  is  larger  for  smaller  wedges  and  the 
i  ..terferogram  was  taken  when  the  shock  system  almost  reached  the  crest  of 
the  wedge,  the  Mach  stem  is  much  longer  and  clearer  than  in  the  previous 
case.  The  slipstream  is  seen  to  be  straight  from  the  triple  point  and  then 
it  curls  up  into  a  vortex  structure  near  the  wedge  surface.  Also  seen 

clearly  in  the  interferogram  is  the  shock-wave  boundary-layer  interaction  at 
the  corner.  The  reflected  shock  wave  is  curved  and  attached  to  the  wedge 
surface  at  the  corner.  Another  shock  which  is  running  ahead  of  the 
reflected  shock  at  the  wedge  corner  and  this  system  is  referred  to  as  a 
lambda-shock  configuration,  which  will  not  be  studied  in  this  work.  The 

portion  of  the  reflected  shock  immediately  behind  the  triple  point  is  fairly 
straight  before  it  curves  downward.  It  can  be  expected  that  this 

experiment  is  close  to  the  transition  boundary  from  SMR  to  CMR  since  the 

reflected  shock  is  straight  at  the  triple  point  and  the  closeness  of  the 
angle  6  to  90°. 
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The  wall-density  distribution  is  presented  in  Fig.  21.  The  density  has 
an  increase  as  it  crosses  the  slipstream  from  state  (3)  to  state  (2),  and  it 
remains  steady  over  a  large  region.  As  one  approaches  the  wedge  corner,  the 
density  increases  rapidly  and  drops  suddenly  when  crossing  the  stem  of  the 
reflected  shock.  Then  the  density  increases  sharply  again  before  it  drops 
discontinuously  back  to  the  value  of  state  (1)  as  one  crosses  the  outer  stem 
of  the  lambda-shock  configuration. 
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3.2.3  Complex-Mach  Reflection 

Case  5  :  Ms  =  1.72,  Ow  =  40°,  PQ  =  226  torn  (30.13  kPa),  TQ  =  24.9°C 

An  excellent  infinite-fringe  i nterferograin  of  a  detached  CMR  is  shown 
in  Fig.  22.  The  clear  and  enhanced  isopycnics  are  the  result  of  a  high 
initial  pressure  and  a  large  wedge  angle.  Note  the  shallow  kink  is  situated 
at  the  intersection  of  the  straight  reflected  shock  and  the  curved  reflected 
shock.  The  flow  behind  the  reflected  shock  between  the  triple  point  and 
the  kink  is  supersonic  relative  to  the  triple  point,  causing  the  reflected 
shock  to  be  straight  in  that  region.  Note  the  spreading  of  the  straight 
slipstream  and  its  termination  at  the  wedge  surface  without  curling  up.  The 
boundary-layer  behind  the  slipstream  can  also  be  seen  in  the  i nterferograin . 
The  isopycnics  near  the  wall  on  the  upstream  side  of  the  slipstream  slant 
down  to  the  left,  whereas  those  on  the  downstream  s’de  slant  down  to  the 
right.  The  boundary  layer  profile  changes  shape  between  the  isopycnics  ' f ’ 
and  ' g '  because  the  highest  density  occurs  right  behind  the  slipstream.  The 
wall-density  distribution  is  plotted  in  Fig.  23.  The  density  jumps  behind 
the  Mach  stem  and  rises  up  gradually  till  one  passes  the  slipstream  and  then 
it  decreases  to  a  minimum  slowly.  Close  to  the  corner  on  the  wedge  surface, 
the  density  rises  to  a  maximum  and  it  drops  at  the  vortex.  Finally,  the 
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density  has  a  sharp  discontinuous  drop  as  one  goes  through  the  reflected 
shock  back  to  state  (1). 

Case  6  :  M$  =  3.14,  Ow  =  10°,  PQ  =  7.95  torr  (1.06  kPa),  TQ  =  24.3°C 

An  example  of  an  attached  CMR  along  with  the  evaluation  of  the  density 
field  is  presented  in  Fig.  24a.  The  number  of  isopycnics  in  the  first 
harmonic  interferogram  is  dramatically  reduced  because  of  the  relatively  low 
initial  pressure  used.  Therefore,  the  second  harmonic  interferogram  is 
also  being  analysed  in  Fig.  24b,  due  to  the  fact  that  it  provides  twice  as 
many  isopycnics  as  the  first  harmonic.  There  is  a  smooth  reversal  of 
curvature  of  the  reflected  shock  wave  and  it  is  difficult  to  locate  exactly 
the  position  of  the  kink.  The  slipstream  is  straight  emanating  from  the 
triple  point  and  is  pushed  forward  near  the  wedge  surface  toward  the  Mach 

stem  by  the  higher  pressure  behind  the  curved  reflected  shock.  The  Mach 

stem  is  seen  slightly  curved  in  this  picture.  The  bow  shock  wave 
terminates  at  the  corner  straight  and  attached,  and  no  lambda  shock  occurs. 
The  wall-density  distribution  is  shown  in  Fig.  25.  The  density  jumps 
behind  the  Mach  stem  and  has  a  slight  rise  across  the  slipstream  before  it 
increases  gradually  to  a  fairly  constant  density  region  near  the  wedge 
corner.  The  density  then  increases  at  the  corner  and  drops  sharply  as  one 

goes  through  the  reflected  bow  shock  back  to  state  (1). 

3.2.4  Double-Mach  Reflection 

Case  7  :  Ms  =  2.18,  Ow  =  37°,  PQ  =  35.2  torr  (4.69  kPa),  TQ  =  23.7°C 

An  infinite-fringe  interferogram  for  the  case  of  an  attached  DMR  is 
shown  in  Fig.  26.  The  sharp  kink  now  becomes  the  second  triple  point.  The 
band  of  compression  waves,  as  indicated  by  the  enhanced  isopycnics  emanating 
from  the  second  triple  point  to  the  wedge  surface  interacting  with  the 
slipstream,  converges  and  merges  with  the  curved  reflected  shock  forming  a 
minute  second  Mach  shock.  The  slipstream  develops  instabilities,  forming 
small  ripples  before  it  rolls  up  at  the  wedge  surface.  Since  the  second 
Mach  shock  is  weak  and  not  yet  fully  developed,  the  second  slipstream  is 
not  recorded  in  the  i nterf erogram .  At  the  wedge  corner,  the  intense 
lambda-shock  interaction  is  prominent.  In  Fig.  27,  the  corresponding 

wall-density  distribution  is  presented.  After  the  increase  of  density 
across  the  Mach  stem,  the  density  remains  steady  in  state  (3)  before  it 
jumps  rapidly  to  a  maximum  density  across  the  slipstream.  Then  it  decreases 
gradually  and  remains  constant  till  one  moves  further  down  closer  to  the 
wedge  corner  where  the  density  increases  almost  back  to  the  maximum  value 
again.  The  strong  shock-wave  boundary-layer-interaction  effect  at  the 
corner  is  also  clearly  indicated  in  the  wall-density  profile.  As  one 
crosses  the  two  shocks  of  the  lambda  shocks  on  the  wall,  the  density 

experiences  a  sharp  decrease  reaching  a  minimum  at  the  wedge  corner,  and 
then  the  density  rises  slightly  and  falls  sharply  back  to  state  (1). 

Case  8  :  M$  -  2.59,  ow  =  30°,  PQ  =  11.2  torr  (1.49  kPa),  TQ  =  24.2°C 

A  good  example  of  attached  DMR  with  the  evaluation  of  the  density  field 
is  shown  in  Fig.  28.  Although  the  experiment  was  performed  with  relatively 
low  initial  pressure,  there  is  sufficient  information  in  the  interferogram 
to  study  the  density  field.  Two  fairly  uniform  regions  are  seen  behind  the 

straight  reflected  shock  and  behind  the  Mach  stem.  Ti.o  straight  slipstream 


gradually  develops  turbulence  and  is  pushed  forward  to  curl  up  toward  the 
Mach  stem.  A  few  isopycnics  behind  the  curled  up  slipstream  can  be  seen 
clearly  converging  to  form  a  short  second  Mach  shock  which  merges  with  the 
curved  reflected  shock  at  the  second  triple  point.  The  second  slipstream  is 
again  too  weak  to  be  recorded  in  the  interferogram.  Note  how  sharp  and 
detailed  the  isopycnics  are  recorded  of  the  lambda-shock  interaction  at  the 
wedge  corner.  The  wall-density  distribution  for  this  case  appears  in  Fig. 
29.  The  effect  of  the  rolling  up  of  the  slipstream  and  the  compression  wave 
is  illustrated  in  the  wall-density  profile  as  a  sharp  density  jump  after  the 
initial  density  increase  across  the  Mach  stem.  The  effect  of  the 
lambda- shock  interaction  is  more  dramatic  in  this  c=»se.  As  one  moves  along 
the  wedge  surface  down  to  the  corner,  the  density  rises  to  a  maximum  and 
then  drops  sharply  as  one  crosses  the  inner  lambda  shock.  The  density  then 
rises  again  before  it  drops  sharply  back  to  state  (1)  as  one  crosses  the 
outer  lambda  shock. 

3.2.5  Mach  Reflections  in  Ar,  air,  C09  and  SF^ 

Extensive  research  of  shock  wave  reflection  in  various  types  of  gases, 
such  as  monatomic  argon,  diatomic  02,  N2  and  air,  triatomic  C02  and 

polyatomic  SF(  ,  have  been  carried  on  at  UTIAS  for  more  than  a  decade.'  Lee 

and  Glass  L2(>J  did  some  numerical  work  in  showing  the  effects  of  the 

specific  heat  ratio  y  on  the  transition  boundaries  and  the  shock-wave 
configuration.  Their  calculations  show  that  the  transition  boundary  shifts 
downward  significantly  with  decreasing  y,  and  that  a  decrease  in  y  lowers 
the  value  of  the  reflection  angle  u*  .  It  became  necessary  to  compare  some 
experiments  done  in  gases  with  different  y.  Four  experimental  results  of  MR 
in  four  different  types  of  gases  with  the  same  wedge  angle  of  GL  =  20°  and  a 
Mach  number  of  M$  =  3.4  ±  0.2  are  gathered  for  a  comparison  of  their  types 

of  reflection,  shock  shapes,  wall-density  distribution  profiles,  the 

differences  in  the  lambda- shock  interaction,  and  the  shape  of  the  slipstream 
and  its  effect  on  the  Mach  stem.  The  initial  conditions  of  the  four  runs 
are  listed  below. 

Case  9: 

Monatomic  Ar,  M  =3.20,  GL=20.0°,  Po=100.0  torr  (13.3  kPa)  ,  T  =24.7°C  [32] 

Diatomic  air,  M  =3.60,  0=20.0°,  P0=  60.0  torr  (8.0  kPa) ,  T0=25.1°C  [26] 

Triatomic  C02,  Ml=3 .43 ,  0  =20.3°,  P0=  30.0  torr  (4.0  kPa),  T=20.8°C  [24] 

Polyatomic  SFg,  Ms=3.51,  ^=20.0°,  PQ=  4.6  torr  (0.6  kPa) ,  TQ=24.0°C 

The  first  and  second  harmonic  infinite-fringe  interferograms  for  the 
run  in  Ar  are  shown  in  Figs.  30(a)  and  30(b),  respectively.  The  type  of 
reflection  is  a  detached  SMR.  Because  the  reflected  bow  shock  is  curved 
throughout  its  entire  length,  the  corner  signals  are  able  to  propagate 
upstream  all  the  way  to  the  triple  point.  No  significant  shock-wave 
boundary-layer  interaction  near  the  corner  is  observed.  The  slipstream  is 
seen  to  extend  from  the  triple  point  and  is  laminar.  Near  the  wedge 
surface,  the  slipstream  is  pushed  ahead  slightly  and  curls  upstream  because 
of  the  high  pressure  behind  the  bow  shock.  Since  the  slipstream  terminates 
at  the  wall  well  Dehind  the  Mach  stem,  it  does  not  disturb  the  flow  in 
region  (3).  Hence,  the  Mach  stem  is  straight  and  perpendicular  to  the  wall. 
The  wall-density  distribution  is  plotted  in  Fig.  31(a)  for  this  run.  There 
is  a  small  density  jump  across  the  Mach  stem  and  the  density  behind  the  Mach 
stem  is  fairly  constant  along  the  wall,  except  for  the  slight  increase  and 
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then  gradual  fall  of  the  density  before  crossing  the  bow  shock  near  the 
corner  into  state  (1). 

The  second  run  is  an  attached  CMR  in  air  as  shown  in  Fig.  3U(c).  The 
reflected  shock  near  the  triple  point  is  fairly  straight  and  then  it  curves 
up  with  the  bow  shock  terminating  at  the  wall.  The  kink  is  not  well  defined 

in  this  run.  The  bow  shock  wave  in  the  air  run  is  not  as  bulged  as  in  the 

Ar  run  and  a  shock-wave  boundary-layer  interaction  at  the  corner  is 
beginning  to  occur.  The  slipstream  is  still  laminar  but  the  curling 

upstream  is  more  apparent.  The  gap  between  the  Mach  stem  and  the  slipstream 
has  decreased.  The  Mach  stem  remains  straight  and  perpendicular  to  the 
wedge  surface.  The  wall-density  distribution  is  plotted  in  Fig.  31(b).  The 
relative  density  jumps  across  the  Mach  stem  and  across  the  slipstream  are 
higher  in  the  air  run  than  in  the  Ar  run.  The  fairly  constant  density 

region  near  the  center  diminishes  and  near  the  corner,  the  density  increases 
swiftly  to  a  maximum  before  it  drops  rapidly  to  state  (1)  across  the  bow 
shock  with  a  shock-wave  boundary-layer  interaction. 

The  third  run  is  an  attached  CMR  in  C02  with  the  finite- fringe 
interferogram  shown  in  Fig.  30(d).  The  vibrational  relaxation  behind  the 
incident  shock  is  clearly  seen  in  the  picture  and  the  relaxation  length  is 
comparatively  much  shorter  than  the  straight  portion  of  the  reflected  shock 
near  the  triple  point  and  the  flow  is  in  equilibrium  behind  the  relaxation 
zone.  The  bow  shock  has  moved  closer  to  the  wedge.  The  kink  is  more 
distinct  here  than  in  the  air  run.  The  shock-wave  boundary-layer 
interaction  is  strong  and  the  lambda-shock  configuration  is. more  well 
defined.  The  slipstream  curls  up  at  the  wedge  surface  in  a  vortical  flow, 
influencing  the  flow  in  region  (3)  pushing  the  Mach  stem  further  upstream, 
tnus  the  Mach  stein  becomes  curved.  Because  the  interferogram  was  not  taken 
in  the  inf  in  i  te-fr  inge  mode,  the  wall-density  profile  cannot  be  obtained 
directly.  Analysis  would  require  the  reduction  of  the  finite-fringe 
i nterterogram  which  was  not  done  for  the  present  study. 

Tne  tint  and  second  harmonic  interferograms  for  the  fourth  run  of  an 
nranifd  i)MR  in  $F6  are  shown  in  Figs.  30(e)  and  30(f),  respectively.  The 
r,. fleeted  shock  is  seen  closer  to  the  wedge  surface  when  compared  to  the  C02 

re  w;:n  a  visible  uniform  region  (2)  behind  it.  There  is  a  very  strong  ana 

:  ♦  jnn  r  1  ow  region  close  to  the  corner  behind  the  bow  shock  and  the 

' nt er dt t i un  of  the  boundary-1 ayer  with  the  straight  attached  bow  shock  is 
reiNdrudhle  here.  The  outer-most  shock  of  the  lambda-shock  system  slants  out 
further  and  lies  closer  to  the  wall  here  than  in  the  previous  runs  of  air 
and  CO,.  The  bow  shock  wave  also  lies  closest  to  the  wedge  among  all  four 

runs.  ‘The  curling  up  effect  of  the  laminar  slipstream  is  the  greatest  here 
also.  The  strong  vortical  flow  of  the  slipstream  is  situated  very  close  to 

tne  Macn  stem.  It  pushes  the  Mach  stem  further  upstream  with  a  bulge 

clearly  seen  in  the  i  nter  ferograin .  Figure  31(c)  illustrates  the 
wall-density  distribution  of  this  run.  The  first  jump  of  density  across  the 
Mach  stem  and  the  second  jump  across  the  slipstream  are  much  larger  than  in 
the  air  run.  Then  it  increases  gradually  as  one  moves  toward  the  wedge 
corner  reaching  a  uniform  region  at  the  corner,  with  the  maximum  density 

occuring  there.  Finally,  there  is  a  very  sharp  jump  across  the  lambda 

shocks  back  to  state  (l). 

Some  statements  can  be  made  from  the  observations  taken  of  the  four 

runs  in  the  various  gases  with  the  same  wedge  angle  and  approximately  the 


same  Mach  number.  Firstly,  the  statement  of  Lee  and  Glass  [25]  that  the 
transition  boundary  lines  shift  downwards  significantly  and  u>'  decreases 
with  decreasing  y  are  supported  by  the  comparison.  The  flows  in  Ar  and  air 
are  both  frozen,  the  flow  in  C02  is  almost  in  equilibrium,  whereas  the  flow 
in  SF6  is  in  equilibrium  under  the  run  conditions  given  above.  For  frozen 
Ar,  y  =  1.667;  frozen  air,  y  =  1.400;  nearly  equilibrium  C02,  1.143  <  y  < 
1.290  and  equilibrium  SF6,  1.056  <  y  <  1.093.  Ar  has  the  largest  y  and  the 
type  of  reflection  here  is  SMR.  Air  has  a  smaller  y  than  Ar  and  the  type  of 
reflection  becomes  a  CMR.  This  implies  that  the  transition  line  of  SMR  <-> 
CMR  has  shifted  down  in  the  transition  map  with  decreasing  y  as  shown  in 
Fig.  31(d).  If  C02  is  used,  y  is  further  lowered.  However,  the  type  of 
reflection  remains  as  CMR  but  is  a  stronger  one  than  in  the  air  run.  In 
SF6,  y  is  the  smallest  and  the  type  of  reflection  is  a  DMR.  This 
corresponds  to  the  CMR<->L)MR  transition  boundary  is  lowered  as  y  is  decreased 
as  shown  in  Fig.  31(e). 

Secondly,  the  claim  of  Lee  and  Glass  [25]  that  w1  decreases  with 
decreasing  y  is  also  substantiated  here.  Listed  below  are  the  measured 
physical  quantities  for  the  four  runs. 

Gas  _6  _x  (angles  are  in  degrees) 

Ar  83.4  40.9  14.3 

Air  99.5  23.1  12.6 

C02  110.8  10.5  11.0 

SF6  127.0  -  5.0  10.4 

The  angle  6  behaves  just  the  opposite  to  u>‘ .  Therefore  lowering  y  will 
increase  6  but  decrease  w' .  The  angle  x  is  also  found  to  be  a  decreasing 
function  of  y.  The  behaviour  of  these  measured  angles  in  different  y  gases 
is  similar  to  increasing  at  the  same  Mach  number  for  a  given  y  gas. 

Thirdly,  the  wall -density  profiles  show  more  drastic  changes  as  y  is 
lowered.  The  density  jumps  across  the  Mach  stem,  slipstream,  bow  shock  or 
lambda  shocks  at  the  corner,  and  the  maximum  density  ratio  attained  behind 
the  slipstream  and  at  the  corner  are  all  increased  with  decreased  y.  Sinde 
pressure  is  directly  related  to  density,  the  pressure  profiles  would  behave 
in  a  similar  manner.  Therefore  the  lower  the  y  of  the  gas  is,  the  stronger 
the  pressure  behind  the  bow  shock  will  be,  and  this  is  evidenced  by  the 
shape  of  the  bow  shock  and  where  it  terminates  at  the  wall. 

Fourthly,  the  strength  of  the  shock-wave  boundary-layer  interaction  is 
observed  to  increase  with  decreasing  y.  In  the  Ar  run,  only  the  boundary 
layer  is  observed  near  the  corner.  However,  definite  shock-wave 
boundary-layer  interactions  at  the  corner  with  the  attached  CMR  is  observed 
in  the  air  run.  The  interaction  turns  into  the  lambda-shock  configuration 
in  both  the  CU2  and  SF6  runs,  with  the  SFfe  run  being  the  stronger  one. 

Lastly,  the  curling  effect  of  the  slipstream  and  the  curving  of  the 
Mach  stem  are  enhanced  with  increasing  y.  The  laminar  slipstream  curls  up 
slightly  in  the  Ar  and  air  runs.  It  curls  into  a  vortical  flow  in  the  C0^ 
run  and  is  more  pronounced  in  the  SF6  run.  The  distance  between  the  Mach 
stem  and  the  slipstream  decreases  with  lowering  y.  The  Mach  stems  are  found 
straight  for  both  the  Ar  and  air  runs.  The  slipstream  is  getting  closer  to 


the  Mach  stem  in  the  C02  and  SF6  runs  and  as  a  result  the  Mach  stem  is 
pushed  out  and  curved  with  the  SFg  run  being  the  more  drastic  example. 

4.0  RESULTS  AND  DISCUSSIONS 


For  the  present  experimental  investigation  of  pseudo-stationary  oblique 
shock-wave  reflections  in  SF.,  over  ISO  experiments  were  performed  in  the 
UT1AS  Hypervelocity  Shock  Tube.  A  listing  of  the  initial  conditions  for 
each  experiment  performed  in  this  study  can  be  found  in  Table  2.  Since 
several  postulates  were  raised  on  transition  criteria  as  described  in 
Sec.  2.2,  it  is  necessary  to  compare  them  with  experimental  results  and 
justify  their  validity.  The  transition  boundary  maps  in  the  (M  -  Q  )  and 
the  (Ms  -  G^J  planes  for  the  case  of  frozen  and  vibrational  equilibrium  SF6 
were  computed  according  to  the  criteria  proposed  by  Ben-Dor  and  Glass 
[13]  and  Shirouzu  and  Glass  [16].  To  verify  the  contribution  of  real-gas 
effects  and  the  validity  of  the  transition  criteria,  experiments  were  done 
well  inside  the  transition  region.  Comparisons  of  some  physical  phenomena 
of  shock-wave  reflection  between  experimental  and  numerical  results  will  be 
presented  and  discussed. 

4.1  Transition  Boundaries 


The  computation  of  the  transition  boundaries  in  this  study  uses  the 
equations  r r  state  for  frozen  SFfe  (v  =  4/3)  and  vibrational  equilibrium  SF^. 
The  thermodynamic  states  of  SFb  used  in  the  calculation  are  considered  in 
more  detail  in  Appendix  B.  In  analyses  of  gasdynamic  phenomena  including 
shock  waves,  relaxation  effects  are  very  important.  At  room  conditions,  SF, 
is  already  partially  vibrat ional ly  excited.  By  examining  the  vibrational 
relaxation  lengths  behind  the  shock  waves,  analysis  shows  that  the  flow 
states  behind  the  shock  fronts  under  the  present  experimental  initial 
conditions  are  in  vibrational  equilibrium.  The  vibrational  relaxation 
effect  in  SF6  is  discussed  in  Appendix  C.  It  will  be  seen  that  relaxation 
lengths  play  an  important  role  in  the  interpretation  of  the  transition 
boundaries  and  the  regions  they  enclose. 

The  experimental  results  are  classified  into  the  four  possible  types  of 
reflection  and  are  plotted  in  their  transition  maps  for  comparison.  The 
comparisons  are  twofold:  first,  to  determine  whether  the  experimental 

results  agree  with  the  frozen-yas  or  the  equilibrium-gas  analysis;  second, 
to  verify  tne  validity  of  the  transition  criteria.  The  transition 
boundaries  for  frozen  SF6  in  the  (M  -  0  )  and  (M  -  (V)  planes  are  shown  in 
Figs.  32(a)  and  32(b),  respectively.  In  these  transition  boundary  plots, 
the  transition  between  SMR  and  CMR  is  represented  by  two  lines.  The  dashed 
line  shows  the  former  criterion  M2j  =  1  while  the  solid  line  shows  the  new 
additional  criterion  6  =  yu°  (see^Sec.  2.2.2).  Transition  from  SMR  to  CMR 
will  occur  when  both  the  former  and  the  new  additional  criteria  are 
satisfied.  Note  that  the  two  SMR«-»CMR  lines  cross  over.  The  new  transition 
criterion  to  the  left  of  this  cross-over  point  is  the  dashed  line;  whereas 
to  the  riyht  of  the  cross-over  point,  the  solid  line  is  the  new  transition 
criterion.  In  (Ms  -  n^)  plane,  there  is  an  additional  boundary  line  in  the 
SMR  region.  It  corresponds  to  the  state  of  Mj  =  1.0,  below  which  no 
reflection  NR  can  occur.  It  can  be  seen  in  both  fiyures  that  there  is  no 
agreement  between  the  analytical  and  experimental  results.  However,  there 
is  definite  agreement  with  the  equilibrium  SF^  model  as  shown  in  Figs.  33(a) 
and  33(b). 
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(3  runs)  for  9  =  42°,  and  =  b.48  (1  run)  and  7. do  (2  runs)  for 
Note  that  tne  RR<-41R  boundary  is  predicted  by  the  detachment 
and  is  the  limit  for  the  two-shock  theory  to  have  solutions. 
Therefore,  any  RR  which  exists  beyond  this  boundary  limited  by  the 
two-shock  theory  is  a  persistence  of  RR  or  is  known  as  the  von  Neumann 
paradox.  Furthermore,  the  mechanical  equilibrium  criterion  is  the  limit  for 
the  three-shock  theory  to  have  solutions,  thus  any  MR  lying  above  the  RR <  *MR 
boundary  into  RR  domain,  on  the  contrary,  is  not  a  persistence  of  MR  because 
it  has  not  exceeded  the  limit  of  the  three- shock  theory  for  a  MR  to 
exist.  However,  the  three  DMR  points  lying  in  the  RR  domain  are  in  the 
transition  region  and  very  close  to  the  RR<-*MR  boundary  line.  Therefore, 

they  become  very  sensitive  to  any  conditions  affecting  their  type  of 
reflection.  It  can  be  seen  that  the  termination  of  RR  line,  according  to 
the  detachment  criterion,  is  in  general  very  good  for  engineering 
applications,  however,  improvement  is  still  necessary. 


In  the  (M$  -  9^ )  plane,  the  results  for  do  not  require  the 
measurement  of  y.  Transition  boundaries  in  the  MR  region  are  more  accurate 
in  the  (M,  -  0^)  plane  than  in  the  (M$  -  9^)  plane.  In  Fig.  33(b),  there 
are  two  SMR  points  lying  in  the  CMR  according  to  the  former  criterion 
M2j  =  1.  But  according  to  the  new  SMR<->CMR  criterion,  there  is  only  one  CMR 
data  point  lying  below  the  transition  line,  at  M$  =  3.47  and  =  19.5°. 
Therefore  the  new  criterion  for  SMR<-->CMR  transition  agrees  best  with 
experiment.  There  are  also  some  di screpanci es  in  the  transition  boundaries 
between  CMR  and  DMR.  At  low  Mach  numbers  and  high  wedge  angles,  <  1.95 
and  O'  >  42°,  there  are  several  DMR  points  lying  beyond  the  CMR*--4)MR 
transition  line  (M2K  =  1.0)  into  the  CMR  region.  On  the  other  hand,  there 
is  one  CMR  point  just  lying  above  the  transition  line  in  the  DMR  domain,  at 
M$  =  5.5  and  O'  =  17.8°.  Generally,  except  for  low  Mach  numbers  and  high 
wedge  angles,  tne  CMR  ♦-♦DMR  transition  criterion  is  reasonable,  though  not 
precise. 


In  some  cases,  the  disagreement  between  experiments  and  the  transition 
boundary  of  CMR  and  DMR  are  due  to  experimental  impreciseness  through  the 
use  of  interferograms.  An  i nterferogram  responds  to  changes  in  density 
whereas  a  shadowgram  responds  to  the  second  derivative  of  the  density 
changes.  Thereby,  a  shock-wave-reflection  configuration  is  much  easier  to 
interpret  in  a  shadowgram  than  in  an  interferogram ,  although  it  is  true  that 
an  interferogram  provides  more  quantitative  information.  Consequently,  it 
is  advisable  to  employ  shadowgram  near  the  transition  region  of  CMR  and  DMR 
to  give  more  precise  conf igurations  for  the  classification  of  the  type  of 
reflection. 


As  mentioned  earlier,  the  transition  of  RR<->MR  should  be  compared  in 
the  (Ms-r>w)  plane,  whereas  the  SMR^+CMR  and  CMR  ♦-♦DMR  transitions  should  be 
compared  in  the  (Ms  -  9^)  plane.  If  an  accurate  prediction  of  y  can  be 
found  in  the  future,  then  the  (M  -  9W)  plane  could  be  used  throughout  for 
comparison  of  all  transition  boundaries.  Note  that  the  (M$  -  0^)  plane  has 
another  di sadvantage ,  beside  the  multi-valued  portion  near  the  RR<~>MR 


boundary,  that  the  transition  boundary  lines  being  too  close  together, 
expecially  at  high  Mach  numbers,  and  it  becomes  difficult  to  determine  if 
experiment  and  analysis  of  the  various  domains  are  in  good  agreement.  In 
Fig.  33(b),  former  SMR<->CMR  transition  line  approaches  the  NR  line  for 
>  b.  This  means  that  the  flow  Mach  number  M1  behind  the  incident  shock 
wave  diminishes  to  nearly  unity  and  the  reflected  shock  occurs  as  a  very 
weak  shock  wave.  Also  in  the  same  figure,  it  can  be  seen  that  the  CMR«->DMR 
and  the  new  additional  SMR<->CMR  transition  lines  approach  the  NR  and  level 
off  at  >  8,  thus  narrowing  the  domains  of  SMR  and  CMR,  and  any  comparison 
in  this  range  would  be  difficult. 

It  is  evident  from  Fig.  33(a)  to  33(b)  that  real-gas  effects  shift 
every  boundary  line  downwards  in  both  the  (Ms  -  m  )  and  the  (M$  -  0^) 
planes.  The  equilibrium  boundaries  begin  to  drop  from  their  frozen-gas 
values  even  at  very  low  Mach  numbers,  due  to  vibrational  excitation.  The 
divergence  between  the  frozen  and  equilibrium  boundaries  becomes  greater 
with  increasing  Mach  number.  At  higher  Mach  numbers,  M$  >  5,  the  frozen 
RR+-+MR  boundary  line  levels  off  so  that  it  is  almost  independent  of  the 
incident  shock  Mach  number.  However,  the  equilibrium  RR<-+MR  transition  line 
strongly  depends  on  the  incident  shock  Mach  number  due  to  real-gas  effect. 
Both  the  SMR«-»CMR  and  the  CMR<-»DMR  transition  boundaries  drop  rapidly  with 
increasing  Mach  number  in  the  case  of  equilibrium  SF6.  Thus,  for 
equilibrium  SF6,  the  domains  of  SMR  and  CMR  are  diminished  while  the  domains 
of  RR  and  DMR  are  enlarged.  The  behaviour  of  the  boundary  lines  due  to 
real-gas  effects  is  analogous  to  the  results  for  gases  with  lower  values  of 
specific  heat  ratio  y  as  reported  by  Lee  and  Glass  [35]. 


Persistence  of  Regular  Reflection 


Several  previous  experimental  studies  (Smith  [7],  Bleakney  and  Taub 
[6],  White  [8],  Kawamura  and  Saito  [41],  Henderson  and  Lozzi  L 10] )  have 
observed  that  RR  exists  beyond  the  detachment  criterion.  This  is  referred 
to  as  the  persistence  of  RR  or  von  Neumann  paradox  and  the  phenomenon  has 
not  been  clearly  understood  or  explained.  Recently,  Shirouzu  and  Glass  [16] 
analysed  the  data  of  Ando  and  Glass  [30]  and  Matsuo  et  al  [42]  in  C02  and 
noted  the  same  trend.  Experimental  results  of  Deschambault  [26]  in  air  show 
that  RR  persists  not  only  at  low  incident  shock  Mach  numbers,  but  also  at 
extremely  high  Mach  numbers. 


Consider  a  transition  map  in  the  (M$  -  o^)  plane  (see  Fig.  8).  At  the 
RR  boundary  of  the  MR  domain,  y  has  a  finite  value  not  equal  to  zero.  Thus 
0*  =  Ow  +  x  has  a  greater  value  than  ow  at  the  boundary.  On  the  contrary, 
0^  is  equivalent  to  n  in  the  RR  domain.  Therefore,  the  RR<-+MR  transition 
boundary  corresponds  to  two  different  lines,  ie.  n  and  ow  +  y,  in  the 
(Ms  -  0^)  plane.  Consequently,  a  point  lying  between  these  two  lines  in  the 
(M$  -  f)^)  plane  corresponds  to  two  physical  points,  one  in  RR  and  the  other 
in  MR.  In  order  to  avoid  this  multi-valued  portion  near  the  RR*->MR 
boundary,  the  transition  between  RR  and  MR  is  best  studied  in  the  (M$  -  Ow) 
plane.  The  persistence  of  RR  has  also  been  found  in  the  present  study  at 
f>w  =  45°  and  low  incident  shock  wave  Mach  number  Ms  =  1.52  ,  1.66,  2.U7  and 
2.11  in  Fig.  33(a). 


Henderson  and  Lozzi  [10]  tried  to  explain  the  anomaly  as  the  apparent 
persistence  of  the  RR  system  is  really  due  to  a  OMR  system  that  has  not 
developed  sufficiently  to  be  observable  as  such.  However,  Ando  [24] 


disagreed  and  remarked  that  the  resolution  of  the  interferoyrams  (0.2  mm) 
would  enable  the  angle  x  of  the  range  encountered  in  his  experiments  (1.5° 
to  2.5°)  to  be  clearly  seen.  Deschambault  [26]  showed  that  very  small 
triple-point  trajectory  angles  of  less  than  1°  (far  smaller  than  those  of 
the  nascent  DMR)  were  observed.  In  the  present  work,  the  values  of  x  less 
than  1°  were  also  observed.  For  example,  the  infinite-fringe  interferogram 
in  Fig.  34  shows  a  double-Mach  reflection  with  x  =  0.7°  and  part  of  the 
first  slipstream  can  be  seen  clearly  as  a  white  line  along  the  wedge 
surface.  Deschambault  [26]  indicated  that  persistence  in  some  cases,  may 
only  be  a  limitation  of  the  resolution  of  the  experimental  apparatus.  He 
reported  one  case  where  the  reflection  was  classified  as  RR  when  the 
incident  shock  wave  was  only  60%  up  the  wedge,  and  a  run  under  similar 
initial  conditions  with  the  incident  shock  wave  80%  past  the  wedge  corner, 
the  reflection  was  a  DMR.  It  becomes  desirable  to  catch  the  incident  shock 
wave  as  far  up  the  wedge  as  possible  to  observe  fine  details  of  the  flow 
field  more  clearly. 

Hornung  et  al  [23]  conducted  experiments  in  dissociating  gases  and  both 
the  pseudo- steady  and  the  steady  flow  experiments  yielded  results  which 
indicate  that  the  equilibrium  transition  conditions  apply.  They  showed  that 
the  finite  relaxation  length  was  not  sufficient  enough  to  resolve  the 
discrepancy  between  the  analytical  and  experimental  transition  anyle  from  RR 
to  MR.  They  then  suggested  that  viscous  boundary  layer  effect  could  be  used 
to  resolve  the  persistence  of  RR.  Analytically,  the  boundary  layer  imposes 
a  negative  displacement  thickness  and  causes  the  displaced  wedye  to  have  a 
larger  inclination,  although  the  wedge  angle  in  front  of  the  incident  shock 
wave  is  still  equal  to  the  actual  wedge  angle  ew  (see  Fig.  35).  With  the 
displaced  wedge  being  further  away,  the  deflection  through  the  reflected 
shock  necessary  to  direct  the  flow  parallel  to  the  displaced  wall  is 

therefore  less  than  that  required  in  the  inviscid  flow.  The  displacement 
anyle  0d  thus  formed  enables  RR  to  persist  to  lower  wedge  angles.  In 
response  to  the  displacement  effect,  the  angle  made  by  the  reflected  shock 
with  the  wedge,  w'  should  decrease.  Hornung  et  al  [23]  observed  the 

opposite  behaviour  of  w'  in  their  results.  However,  Shirouzu  and  Glass  [16] 
analysed  the  previous  experiments  in  C02  and  air,  and  obtained  results  which 
support  the  viscous  boundary-layer  theory.  More  experiments  are  required 
for  further  analysis  of  this  problem,  not  only  at  the  transition  boundary 
but  around  the  transition  reyion  also  to  provide  information  on  the 

behaviour  of  some  reflection  quantities  such  as  u>'  in  RR  and  x  1n  MR.  At 

present,  an  investigation  of  the  problem  of  persistence  of  RR  due  to  the 
viscous  boundary-layer  effect  in  air  is  being  done  (Wheeler  [43]). 


Transition  Boundary  Between  CMR  and  DMR 


In  the  calculation  of  the  transition  boundary  between  CMR  and  DMR,  the 
location  of  the  kink  K  is  necessary  in  order  to  obtain  the  flow  velocity 
behind  the  reflected  shock  wave  R  with  respect  to  the  kink.  An  assumption 
was  made  by  Law  and  Glass  [15]  that  the  kink  K  has  the  same  horizontal 
velocity  as  the  flow  velocity  behind  the  incident  shock  wave  U^.  As  a 
result,  the  ratio  ?*/L*  (see  Fig.  9a)  is  equal  to  the  density  ratio  across 
the  incident  shock  wave  irrespective  of  the  wedge  anyle,  that  is 


0 


Bazhenova  et  al  L 2U J  investigated  the  ratio  f*/L*  experimentally  and  found 
that  this  ratio  for  a  given  Mach  number  decreased  and  approached  zero  as  the 
wedye  anyle  increased.  Since  the  ratio  f*/L*  vanished  at  the  critical  wedye 
anyle  correspond!  ny  to  the  RR+-»MR  transition,  as  a  result,  it  can  be 
concluded  on  physical  yrounds  that  the  transition  boundaries  of  $MR<->CMk  and 
CMR<-+DMR  should  merye  at  the  RR  boundary  [16].  The  trend  of  f*/L*  is 
plotted  ayainst  wedye  anyle  0  for  three  Mach  numbers  in  Fiy.  36  usiny  the 
present  experimental  results.  Note  ttiat  the  ratio  ?*/ L*  does  approach  zero 
as  the  wedye  anyle  approaches  the  critical  value  at  the  RR  boundary. 

Some  interestiny  points  are  found  in  the  behaviour  of  the  ratio  J*/L*. 
Dashed  lines  are  used  to  link  up  data  points  correspond i ny  to  the  same  Mach 
number  ranye.  All  three  curves  are  seen  to  behave  in  a  similar  fashion.  At 
lower  wedye  anyles,  the  ratio  ?*/L*  increases  with  o  and  reaches  a  maximum 
then  decreases  with  further  increase  in  o  and  eventually  drops  sharply  to 
zero  at  very  high  ow.  It  is  seen  that  as  M$  increases,  the  curves  shift  to 
the  left  and  the  relative  maximum  value  becomes  smaller.  Refer  to  Fiy. 
33(a),  with  the  incident  shock  wave  Mach  number  fixed  at  M  =  4,  all  four 
different  types  of  reflection  will  occur  in  the  sequence  or  SMR,  CMR,  DMR 
and  RR  as  the  wedye  angle  is  raised.  With  SMR  and  RR,  Z*  and  j?*/L*  are  both 
zero.  If  an  experimental  transition  boundary  line  is  drawn  between  SMR  and 
CMR  in  Fiy.  33(a),  SMR  will  just  terminate  at  M$  =  2.1  and  ow  =  20°,  M$  * 
3.0  and  Ow=  8°,  and  M$  =  4.0  and  ow  =  2°.  Since  the  transition  from  SMR  to 
CMR  is  smooth,  the  value  of  ?*/L*  will  beyin  at  zero  and  increase  gradually 
to  the  value  corresponding  to  a  CMR.  If  one  extrapolates  the  curves  of  Ms  = 
3.0  and  M$  =  4.04  at  the  lower  wedge  anyle  side,  the  ratio  Z*/l*  will  be 
zero  at  nw  =  7°  and  1.6°,  respectively,  which  agree  well  with  the 
termination  conditions  of  SMR  on  the  transition  map  found  above.  Consider 
the  curve  of  Ms  =  4,04,  the  ratio  ?*/L*  is  zero  as  SMR  terminates  at  Pw  = 
1.5°.  The  value  of  the  ratio  increases  as  the  wedge  anyle  is  increased,  and 
reaches  a  maximum  at  O  ~  20°.  Then  it  decreases  and  vanishes  at  Ow  =  40° 
which  corresponds  to  where  the  termination  of  DMR  and  RR  occurs.  In 
general ,  the  pressure  in  states  (4)  and  (5)  behind  the  reflected  shock  R'  is 
laryer  than  that  in  states  (2)  and  (3)  behind  the  reflected  shock  R  (see 
Fig.  6),  as  evident  from  the  pressure  histories  taken  by  Bertrand  L44], 
Merritt  [45]  and  Deschambault  [2b].  With  the  increase  in  wedge  angle,  the 
reflected  shock  R'  becomes  more  and  more  normal  to  the  induced  flow  behind 
the  incident  shock  I.  The  pressure  in  states  (4)  and  (5)  become  higher  and 
the  second  triple  point  T‘  is  pushed  further  upstream  and  hence  the  ratio 
J!*/L*  decreases. 

Several  experiments  were  conducted  at  hiyh  wedye  anyle  near  the  RR 
boundary  in  the  CMR  region.  Refer  to  Fiy.  33(a)  and  33(b),  at  Ow  =  45°, 
M$  =  1.81  and  1.91,  DMR  is  observed  in  the  proposed  CMR  domain.  Two  DMR  are 
found  penetrating  into  the  CMR  regime  at  M$  =  1.94,  nw  =  4U°  and  M$  =  1.93, 
Ow  =  42°.  Flence  the  present  results  support  the  postulation  that  the 
transition  lines  meet  at  a  point  on  the  RR  boundary.  Flowever,  the  MR 
transition  boundaries  computed  usiny  the  criteria  suggested  by  Ben-Dor  and 
Glass  [13]  terminate  at  the  RR  boundary  at  two  separate  points.  Further 
analysis  for  the  prediction  of  the  location  of  the  kink  K  at  hiyher  wedye 
anyle  is  required  to  modify  the  transition  line  between  CMR  and  DMR  as  shown 
approx imately  usiny  the  experimental  data  by  the  dot-dash  line  extendiny  to 
the  triple  point  P  on  the  RR<~>MR  boundary  lines. 


4.4  Comparison  of  Analytical  and  Experimental  Results  of  Physical 

Quantities  in  Mach  Reflection 

Analytical  results  based  on  the  frozen-gas  model  is  distinct  from  the 
equilibrium-gas  model.  By  comparing  the  experimental  results  with  the 
frozen  or  equi 1 i bri um-gas  model  for  certain  flow  properties,  the  validity  of 
each  model  can  be  justified.  From  a  practical  point  of  view,  the  chosen 
properties  should  be  measured  accurately  and  easily,  and  it  is  important 
that  they  are  fundamental  quantities  which  are  not  affected  by  other 
assumptions  used  in  the  computations.  The  angle  6  between  the  incident  and 
reflected  shock  waves,  and  the  angle  u>‘  between  the  reflected  shock  wave  and 
the  wall  or  the  triple-point  trajectory  path  are  chosen  as  the  fundamental 
quantities  to  serve  this  purpose  (see  Fig.  37). 


In  the  present  comparison,  the  experimental  values  are  compared  with 
the  analytical  results  for  the  same  instead  of  o  used  in  previous 
comparisons  [24,25],  In  the  three-shock  theory,  the  effective  wedge  angle 
0^  plays  an  important  role.  With  a  given  value  of  at  a  fixed  M$,  the 
directions  of  the  reflected  shock  wave  and  the  Mach  stem  are  determined 
together  with  the  flow  variables.  In  this  calculation,  the  Mach  stem  at  the 
triple  point  is  assumed  perpendi cul ar  to  the  wedge  surface  and  a  value  of  x 
is  predicted.  The  effective  wedge  angle  0^  is  related  to  x  by  0^  =  CL  +  x. 
This  means  that  0^  is  a  function  of  \  and  the  value  of  x  is  affected  by  the 
orientation  of  the  Mach  stem.  The  assumption  of  a  perpendicular  Mach  stem 
to  the  wedge  surface  is  satisfied  precisely  at  the  foot  of  the  Mach  stem  to 
meet  the  necessary  boundary  condition  that  flows  on  both  sides  of  the  Mach 
stem  must  be  parallel  to  the  wedge  surface  with  an  inviscid  flow.  It  is 
also  valid  at  the  triple  point  if  the  Mach  stem  is  perfectly  straight. 
However,  experimental  results  show  that  the  Mach  stem  generally  is  curved. 
Shirouzu  and  Glass  L 16 J  tried  to  take  into  account  the  effect  of  the  curved 
Mach  stem  in  the  prediction  of  x  for  several  gases,  but  the  discrepancies 
are  not  minimized.  It  is  important  to  recognize  that  comparison  of 
experimental  results  and  analytical  results  which  include  the  assumption  of 
a  straight  Mach  stem  are  no  longer  accurate  if  they  are  based  on  the  same 
value  of  0W .  Therefore,  the  experimental  data  should  be  compared  with 
analysis  based  on  the  same  value  of  0^. 

Besides  the  comparison  of  the  fundamental  angles  6  and  to',  the 
experimental  first  triple-point  trajectory  angle  x  and  second  triple-point 
trajectory  angle  x'  are  also  compared  with  analysis.  In  the  following,  the 
angles,  6,  w'and  x' ,  are  compared  twice  with  and  QT  as  the  parameter, 
while  the  angle  x  is  compared  with  ^  as  the  parameter  for  both  the  frozen 
and  equilibrium-gas  models  in  the  range  of  1  <  M$  <  5.5. 

It  is  important  to  note  that  the  physical  quantities  obtained  in  the 
present  work  are  measured  from  photographs.  Consequently,  there  are  errors 
due  to  measurement  and  magnification  of  the  prints.  Errors  involved  in 
measurements  are  ±1.5“  for  6  and  w' ,  and  ±1°  for  x  and  x1  •  A  line  from  the 
reflection  point  P  or  the  first  triple  point  T,  which  is  tangent  to  the 
reflected  shock  R,  needs  to  be  drawn  before  the  angles  6  and  w1  can  be 
measured  (see  Fig.  37).  Unless  region  (2)  is  supersonic  relative  to  the 
reflection  point  P  or  the  first  triple  point  T,  or  region  (2)  is  fairly 
uniform,  the  reflected  shock  is  straight,  otherwise  the  reflected  shock  R 
curves  continuously  till  it  meets  the  incident  shock  wave  I.  This  makes  it 
difficult  to  locate  the  tangent  line  to  R,  thus  causing  the  error  of 


measurement  in  ft  and  «  niyher.  When  o  is  used  as  the  parameter  for 
comparison,  the  numerical  results  are  based  on  the  measured  o'  and  are 
iterated  within  0.1°  of  those  values.  * 

Table  3  presents  the  measured  experimental  data  of  ft,  m' ,  y  an(J 
( =  0  +  x '  ) .  The  frozen  and  equilibrium-gas  analytical  results  at  fixed 

0^  are  also  given.  Note  that  the  calculated  analytical  wedge  angle  o  ]  is 
different  from  the  experimental  r>w  t  and  it  can  be  obtained  by  riWCd]=  0^  - 
Xc a  1  *  where  Xca ]  is  the  calculated  analytical  value  of  y.  For  example,  in 
experiment  no.  118,  both  the  experimental  and  analytical  results  have  the 
same  o^  =  45.7°.  The  experimental  value  of  Ow  =  4b°,  whereas  the 
equilibrium  analysis  has  o  i  =  45.7°  -  1.7°  =  44°,  and  the  frozen  analysis 
has  ciwcal  =  45.7°  -  2.9°  =  42.8°. 


A  numerical  comparison  of  the  physical  quantities  in  MR  between 
experimental  and  analytical  results  at  fixed  Ow  is  tabulated  in  Table  4. 
Note  that  the  analytical  values  of  y  are  different  from  the  experimental 
ones,  and  so  are  the  values  of  o'.  For  example,  in  experiment  no.  118,  both 
the  exper imental  and  analytical  data  have  o^  =  45°.  The  experimental 
X  =  0.7°,  o^  =  4b. 7°,  whereas  the  equilibrium  analysis  has  x  =  1.5°, 
0^  =  46. b°  and  the  frozen  analysis  has  y  =  2.2°  and  o^=  47.2°. 

4.4.1  Angle  Between  Incident  and  Reflected  Shock  Waves  6 

The  comparison  of  ft  is  essentially  equivalent  to  the  comparison  of  to' 
if  an  experimental  result  is  compared  with  a  calculated  result  which  has  the 
same  value  of  o^  since  they  are  related  by 


ft  =  90°  +  0^  -  u' 


(4.2) 


as  shown  in  Fig.  37.  Nevertheless,  the  comparison  of  ft  has  some  advantages 
over  u)1  ,  for  ft  is  a  monotonic  function  of  o  and  o^  at  fixed  Ms,  whereas  m' 
has  a  minimum  point  at  a  certain  ow  and  o^  for  some  part  of  the  Mach  number 
range  [25].  Moreover,  ft  is  used  as  a  new  criterion  for  the  transition 
between  SMR  and  CMR ,  it  is  necessary  to  check  the  agreement  between  measured 
and  calculated  results  of  6. 

Figure  38  shows  the  experimental  results  for  o'  plotted  in  the  (M$  -  ') 
plane.  The  number  beside  each  experimental  point  indicates  the  measured 
value  of  0^ .  The  line  that  extends  above  or  below  a  data  point  shows  the 
discrepancy  of  the  analytical  value  from  the  experimental  result.  The 
calculated  results  are  shown  in  this  figure  for  a  vibrational  equilibrium 
gas  for  each  fixed  o^  in  solid  lines  and  are  used  only  as  a  guide  to 
illustrate  the  behaviour  of  ft  in  relation  to  M^.  It  can  be  seen  that  the 
maximum  discrepancy  found  is  less  than  2.b°.  Thus,  the  agreement  between 
the  experimental  results  and  equilibrium-gas  calculations  is  excellent.  The 
same  good  agreement  is  expected  in  the  comparison  of  the  angle  m'  since  both 
ft  and  oj*  are  fundamental  quantities. 


The  plot  in  Figure  39  uses  ow  instead  of  o'  for  comparison.  The 
numerical  results  are  shown  for  a  vibrational  equilibrium  gas  for  each  fixed 
0W  as  a  solid  line.  The  experimental  results  at  a  fixed  wedge  angle 
supposedly  should  fall  on  top  of  the  corresponding  ow  analytical  line.  The 
agreement  is  not  as  yood  as  obtained  in  Fig.  38.  In  most  cases,  the 
experimental  ft  is  larger  than  the  calculated  result  at  lower  o  At 


high  the  numerical  lines  are  closer  to  one  another  and  dS/d^  is  small. 
Therefore  any  discrepany  at  higher  wedge  angles  is  not  easily  seen  in  the 
figure.  The  disagreement  shown  in  this  plot  illustrates  the  additional 
errors  caused  by  x  when  using  rather  than  r^. 

Figure  40  compares  experimental  results  with  the  frozen-gas  analysis  in 
o'.  The  values  of  corresponding  to  each  data  point  are  identical  to 
tnose  numbered  in  Fig.  38  and  are  not  repeated  in  this  plot.  It  is  seen 
that  the  numerical  lines  are  all  shifted  lower  to  the  right  from  the 
vibrational  equilibrium  case.  Thus,  for  a  given  Mach  number  and  OA,  6  will 
have  a  smaller  value  with  frozen-gas  analysis  than  with  vibrational 
equilibrium  analysis.  For  example,  at  M$  =  5.0  and  0^  =  40°,  6  =  150°  with 
equilibrium-gas  analysis,  while  6  =  126.7°  according  to  frozen-gas 
analysis.  It  is  important  to  realize  that  because  of  such  significant 
difference  between  the  frozen  and  equilibrium-gas  calculations,  it  becomes 
possible  to  conclude  whether  the  frozen  or  equilibrium-gas  model  agrees  best 
with  experiments.  A  similar  trend  of  the  constant  PL  numerical  lines  is 
seen  in  Fig.  41.  Note  that  the  agreement  of  the  frozen-gas  results  with  the 
experiments  is  poor  in  both  and  Ow  comparisons.  Because  of  the  large 
disagreement  with  the  frozen-gas  model,  the  discrepancy  lines  are  not  drawn 
here. 


Since  6  >  90°  is  the  additional  necessary  condition  beside  M2y  >  1  for 
the  transition  from  SMR  to  CMR,  it  is  essential  to  check  the  agreement 

between  the  measured  values  of  5  and  the  type  of  reflection  that  occurs.  A 
plot  of  the  experimental  values  of  6  with  incident  shock  wave  Mach  number  M 
for  the  various  types  of  reflections  is  shown  in  Fig.  42.  There  are  two  SMR 
points  at  ow  =  42°,  Ms  =  1.52  and  1.54  with  6  =  90.8°,  and  another  two  SMR 

points  at  r>w  =  20°,  Ms  =  2.08  and  2.11  which  have  6  of  90.1°  and  90.6°, 

respectively,  lie  above  the  6  -  90°  line.  Any  SMR  points  in  the  6  >  90° 
region  is  not  counterevidence  of  the  new  transition  condition  since  it  is 
only  a  necessary  condition  for  the  transition  and  not  a  sufficient 

condition.  On  the  other  hand,  all  experimental  CMR  but  one  at  -  4°,  M$  = 
3.47  and  6  =  8/°  lies  below  the  6  =  90°  line.  Note  that  a  discrepancy  of  3° 
is  considered  of  the  similar  magnitude  as  an  error  in  measurement.  This 
means  that  the  experimental  results  satisfy  the  condition  6  >  90°  and 

provide  good  evidence  that  this  is  a  necessary  condition.  It  is  worth 
mentioning  that  the  results  presented  in  Fig.  42  do  not  involve  any 
assumptions  or  ca 1 cul at i ons  .  They  illustrate  only  the  experimental  relation 
between  the  transition  from  SMR  to  CMR  and  the  angle  6.  Therefore  this  is  a 

direct  comparison  between  the  new  necessary  condition  and  the  available 

experimental  data. 


4.4.2  Angle  Between  Reflected  Shock  Wave  and  Triple-point  Trajectory  m1 

The  angle  V  is  sometimes  called  the  reflection  angle  as  shown  in  Fig. 
37  and  is  given  for  both  RR  and  MR  as 


£ 


> 


where  is  the  incident  flow  angle  to  the  reflected  shock  wave,  and  "'i  is 
the  flow  deflection  angle  behind  the  incident  shock  wave  (see  Fig.  3a,  4). 
livozdeva  et  al  [21]  calculated  the  reflection  angle  m'  as  a  function  of  the 
incident  shock  velocity  UQ ,  and  the  angle  of  incidence  f0  for  several  gases. 
They  discovered  that  the  increasing  of  incident  shock  velocity  leads  to  the 
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decrease  of  the  reflection  angle  and  its  sliding  into  the  negative  area 
under  certain  conditions.  Recently,  Lee  and  Glass  [25]  studied  the 
variation  of  and  concluded  that  it  is  always  positive  for  Rk,  SMk  and 

CMR.  The  shift  of  from  a  positive  to  a  negative  value  with  higher  Mach 
number  occurs  only  in  the  DM R  domain  for  imperfect  air  and  gases  with  low 
Y. 

Figure  43  shows  the  ( M.  -  ■  •>' )  plot  with  as  the  parameter.  The 
measured  value  of  is  indicated  beside  each  experimental  point  and  the 
solid  lines  are  the  calculated  o'  for  equilibrium  SF6.  The  discrepancy  of 
analytical  value  from  experimenta r  value  is  shown  by  the  extended  vertical 
line  from  the  data  point.  No  more  than  a  discrepancy  of  3.6°  is  found  in 
this  plot.  Thereby  showing  that  the  experiments  agree  very  well  with  the 
equilibrium-gas  results. 

A  similar  plot  is  given  in  Fiy.  44  with  0  used  as  the  comparing 
parameter.  Each  solid  line  corresponds  to  a  fixed  experimental  wedge  angle 
for  an  equilibrium  gas.  In  general  ,  the  experimental  data  of  high  wedge 
angles  lie  slightly  above  the  analytical  line;  whereas  for  low  wedge  angles, 
experimental  results  lie  under  the  analytical  line.  The  maximum  discrepancy 
obtained  by  using  ow  rather  than  0^  is  11°.  It  seems  from  the  figure  that 

u)'  for  high  o  agree  well  with  numerical  lines.  However,  numerical  lines 

are  so  close  to  one  another  that  dm'/dr^  is  small  and  the  discrepancy  is 
not  obvious.  It  is  evident  from  the  figure  that  w'  is  a  decreasing  function 
of  Ms .  At  low  Ms,  the  numerical  lines  are  close  to  vertical  which  implies 

that  m 1  is  a  very  sensitive  function  of  .  Thus  a  small  error  in  Ms  may 

cause  a  large  discrepancy  in  o'  at  low  Ms . 

As  discussed  earlier,  is  not  an  adequate  parameter  to  be  used  in  a 

comparison  until  y  can  be  predicted  accurately.  A  comparison  in  Fig.  44 
should  use  m'  as  the  parameter  instead  of  ow  [46],  As  an  interpretation  of 
an  ex  per imen t a  1  point  at  M  =  5.5  and  nw  =  4°  (w1  =  0.0°),  for  example, 
should  not  be  compared  with  numerical  result  which  has  o  of  4°,  out  with 
that  which  has  v'  of  U.U°.  The  discrepancy  for  the  point  is  not  Am*  =  11°, 
but  \o  ~  2°.  When  comparing  the  results  in  this  way,  it  cannot  be 
concluded  absolutely  that  the  discrepancy  is  large  at  low  Ow  and  small  at 
hiyh  °w.  However,  the  discrepancy  of  =  2°  can  be  read  at  low  ^w,  since 
dm'/dr>w  is  large,  but  becomes  difficult  at  high  o  This  is  one  of  the 
reasons  why  6  has  advantage  over  m*  in  a  comparison.  The  same  method  of 

comparison  can  be  applied  to  5  in  Fig.  39. 

It  is  seen  in  both  Figs.  43  and  44  that  at  higher  M$ ,  m'  shifts  into 
the  negative  region.  The  higher  the  wedge  angle  is,  the  sooner  is  the 
shift.  At  low  M  ,  the  value  of  m1  is  smaller  for  higher  or  o^.  After 

the  shift  of  m1  into  the  negative  region,  the  high  or  o'  numerical  lines 
begin  to  merge  and  cross  over.  Hence  m'  has  a  minimum  value  at  certain 
or  V  tor  some  range  of  M$.  That  is,  for  a  given  M$,  there  exist  two  values 
of  1  or  which  have  tne  same  value  of  m' .  Figure  45  and  46  show  the  (M$ 
-  plots  for  a  frozen-gas  analysis  using  and  ^  ds  parameter, 

respectively.  The  experimental  values  are  not  marked  in  Fig.  45  but 
they  are  identical  to  those  shown  in  Fig.  43.  Note  that  the  numerical  lines 
in  both  figures  are  shifted  up  and  to  the  right  from  the  equ  i  1  i  br  i  urn- gas 
results.  The  reflection  angle  ..>*  stays  always  positive  and  no  cross-over  of 
the  numerical  lines  is  observed  with  ''  in  Fig.  45.  The  numerical  lines  in 
Fig.  4b  all  lie  in  the  positive  m'  region  but  merging  and  crossing  over  of 


the  lines  is  observed  with  0W  as  parameter.  There  are  many  experimental 
data  in  the  negative  to'  region,  however,  no  frozen-gas  analytical  line 
exists  in  the  negative  w1  region.  Thereby,  the  frozen-gas  analysis  gives 
poor  agreement  with  experiments  and  the  discrepancy  lines  are  therefore  not 
indicated  here. 

4.4.3  First  Triple-Point-Trajectory  Angle  y 

When  applying  the  three-shock  theory  to  a  MR  it  is  necessary  to  give  a 
value  of  x  as  an  initial  condition  in  order  that  the  system  of  equations  be 
closed.  Law  and  Glass  [15]  proposed  an  empirical  method  for  predicting  the 
value  of  x  based  on  experimental  observations  that  the  Mach  stem  is  straight 
and  normal  to  the  wedge  surface,  and  introduced  an  additional  independent 
geometrical  relation 

4> 3  =  90°  -  x  (4.4) 

where  <j>3  is  the  incident  flow  angle  to  the  Mach  stem  (see  Fig.  4).  At  Q,  = 
0°,  x  nas  its  maximum  value  for  a  fixed  M$ .  If  M$  is  fixed,  x  is  a 
decreasing  function  with  0W  increasing  and  approaches  zero  at  the  critical 
wedge  angle  for  the  transition  from  MR  to  RR.  The  plot  of  x  as  a  function 
of  Ms  ,  with  the  actual  wedge  angle  ew  as  a  parameter  for  equilibrium  SF6  is 
shown  in  Fig.  47.  For  a  given  wedge  angle  0W,  x  is  a  decreasing  function 
with  increasing  Mach  number  M$ ,  except  for  high  0  at  low  M$.  With  0W  >  37° 
and  Ms  >  2.0,  the  value  of  x  becomes  almost  constant  over  the  entire  region 
of  M$ ,  implying  that  x  is  a  simple  function  of  Gw.  The  experimental  points 
are  also  plotted  in  the  figure.  The  agreement  between  the  analytical 
results  and  experiments  is  good  in  general,  with  discrepancies  of  the  same 
order  of  magnitude  as  the  error  in  measurement.  However,  the  experimental 
points  of  0  =  10°,  20°  and  30°  have  the  same  trend  as  the  higher  Ow  results 
that  they  level  off  and  become  independent  of  the  Mach  number  Ms .  The 
maximum  discrepancy  found  in  x  is  3.4°  at  M  =  3.76  and  0W  =  20s.  The 
relation  between  the  angles  6,  w1  and  x  is  given  Dy 

6+w1  =  90°  +  0W  +  x  (4.5) 

For  a  fixed  Ow,  the  following  applies 

A6  +  Aw1  =  Ax  (4.6) 

Consider  the  experimental  point  at  Ms  =  3.76  and  Qw  =  20°.  The 
differences  between  the  analytical  results  from  the  experimental  results  are 
A6  =  -5.4°  (Fig.  39),  Aw'  =  +2.0°  (Fig.  44)  and  Ax  =  -3.4°  (Fig.  47).  This 
illustrates  the  point  that  if  a  discrepancy  exists  in  x  for  a  data  point, 
each  of  the  ( Ms  -  6)  and  ( M,  -  to'  )  plots  will  not  display  the  same 
discrepancy,  but  their  sum  will.  In  the  above  example,  the  (M,  -  6)  plot 
shows  a  larger  discrepancy  in  6  than  in  x,  while  the  (Ms  -  w' )  plot  shows  a 
smaller  discrepancy  in  «'  than  in  x. 

Figure  48  shows  the  experimental  data  in  the  (M$  -  x)  plot  with  as 
the  parameter  for  the  frozen-gas  analysis.  Note  that  the  analytical  lines 
are  shifted  up  from  the  equilibrium  values  and  x  becomes  nearly  independent 
of  M$  at  high  Mach  numbers,  even  for  low  ow.  In  all  cases,  the  analysis 
gives  a  much  laryer  value  than  the  experimental  value.  Therefore  there  is 
poor  agreement  between  the  frozen-gas  analysis  and  the  experiments. 


As  described  in  Section  4.4,  experimentally,  the  Mach  stem  at  the 
triple  point  is  not  perpendicular  to  the  wedge  surface  because  it  has 
curvature.  This  discrepancy  affects  the  analytical  three-shock  theory 
solutions,  especially  the  value  of  y.  Besides,  the  bow  shock  wave  has  an 
effect  on  the  position  of  the  triple  point.  The  reflected  shock  wave  is  the 
result  of  two  processes;  one  is  a  reflected  shock  wave  emanating  from  the 
triple  point  and  the  other  is  the  bow  shock  wave  caused  by  the 
flow-deflection  process  of  the  flow  induced  by  the  incident  shock  wave  over 
the  wedge  corner  [13].  In  the  three-shock  theory,  the  interaction  of  these 
two  processes  has  not  been  taken  into  account.  At  low  Mach  numbers  or  wedge 
angles,  the  reflected  shock  wave  from  the  triple  point  becomes  weaker  more 
rapidly  than  the  bow  shock  wave.  Hence,  the  position  of  the  entire 
reflected  shock  wave  is  governed  by  the  bow  shock  wave.  In  the  limit  when 

the  bow  shock  is  a  Mach  wave,  the  intersection  point  of  the  bow  shock 

wave  and  the  incident  shock  wave  is  at  [16].  Thus  at  low  M$  and  small 

,  the  intersection  point  of  the  extension  of  the  bow  shock  wave  and  the 
incident  shock  wave  is  closer  to  a  point  where  it  has  o'NR.  This  results  in 
a  smaller  value  of  y  than  that  is  given  by  the  three-shock  theory. 
Therefore  a  better  method  for  the  prediction  of  y  which  takes  into  account 
the  curvature  of  the  Mach  stem  and  the  effects  of  the  bow  shock  wave  is 

necessary. 

4.4.4  Second  Triple-Point-Trajectory  Angle  y' 

The  best  available  analytical  method  to  date  for  the  prediction  of  y' 
wa-  developed  by  Ben-Dor  [11].  He  assumed  in  the  analysis  that  the  second 
triple  point  T*  with  respect  to  the  first  triple  point  T,  moved  with  the 

same  horizontal  velocity  as  the  induced  flow  behind  the  incident  shock  wave 
I  and  deduced  an  empirical  relation  based  on  an  approximation  for  the 
relative  motion  of  the  second  triple  point  as 


-  tan-  I _ LlUi!* _ 

w  couo  -  P0/n1cot(<f1  -  <t>0  -  Oj) 


(4.7) 


This  method  of  predicting  y'  is  applicable  for  a  DMR  as  well  as  for  a 
CMR  in  which  y1  is  defined  as  the  kink  trajectory  angle. 


Unlike  the  angles  !  and  .>*  ,  the  comparison  of  y'  based  on  the  same 
value  of  ow  is  affected  by  the  calculated  wedge  angle  o  ^  ,  since  y'  has  to 
be  measured  from  the  wedge  surface.  From  Section  4.4.3,  it  has  been 
described  that  the  present  method  of  predicting  y  is  not  accurate.  In  order 
to  obtain  the  same  o'  as  the  experiment,  the  calculated  solution  will  give 
values  of  ''wcajand  ycal  which  are  different  from  the  measured  lw  and  y. 
Consequently,  if  the  value  of  is  compared  with  the  measured  y1  from 
experiment,  effectively,  this  means  that  the  comparison  is  done  on  two 
different  wedge  angles.  To  avoid  this  ambiguity,  the  comparison  should  be 
done  on  the  second  confluence  angle  t  which  is  measured  from  the  wedge 
base  instead  of  the  wedge  surface  and  is  defined  as 


■  °w  *  >' 


(4.8) 


Figure  49  "hows  the  (M$  -  plot  with  as  the  parameter  for 
(.qui  1  ibrium-gas  analysis.  The  experimental  data  are  shown  with  their 
measured  indicated  also.  As  revealed  by  the  solid  lines,  ^  is  generally 
a  decreasing  function  with  increasing  Mach  numbers  and  decreasing  ^  , 


except  at  high  Mach  numbers  and  0^  where  a  slight  increase  in  0^  is 
observed,  and  0“  becomes  nearly  independent  at  high  Mach  numbers.  Note  the 
large  discrepancies  of  the  analytical  result  from  the  experimental  result  at 
relatively  low  Mach  numbers  for  a  set  of  fixed  wedge  angle  experiment.  The 
analytical  0"  are  much  larger  than  the  measured  value  at  lower  Mach 
numbers. 

The  disagreement  is  due  to  the  assumption  that  the  distance  j i*  between 
the  two  triple  points  T  and  T'  remains  finite  (see  Fig.  9a).  However,  the 
length  n*  should  decrease  from  a  certain  CMR  value  with  decreasing  M$  and 
eventually  when  SMR  occurs,  f*  =  0.  That  is,  experimentally ,  the  second 
triple  point  T'  approaches  and  merges  with  the  first  triple  point  T  with 
A*/L*  =  0,  x'=  x  and  %  -  %•  This  will  be  discussed  later  in  more  detail. 
The  agreement  improves  significantly  as  M&  is  increased.  Thus  it  can  be 
concluded  that  the  prediction  of  x‘  works  well  for  strong  CMR  and  DMR,  and 
the  experiments  agree  well  with  the  equilibrium-gas  results. 

Figure  50  shows  a  (Ms  -  x‘)  plot  using  Gw  for  the  equilibrium-gas 
model.  It  can  be  seen  that  x'  is  a  decreasing  function  with  increasing  Mach 
number  and  wedge  angle.  At  high  wedge  angles,  ew  =  37°,  40°  and  42°,  the 
values  of  x'  are  shown  as  predicted  by  the  analysis  that  they  decrease  with 
increasing  Mach  number  and  eventually  equal  to  zero. 

The  discrepancies  in  this  figure  can  be  due  to  two  causes.  One  is  that 
analytical  solution  does  not  meet  the  requirement  that  i*/l*  approach  zero 
as  the  Mach  number  decreases.  The  other  is  that  even  though  X*/L* 
approaches  zero  is  taken  into  account,  and  x'  approaches  x,  x'  still  cannot 
be  predicted  accurately  since  x  itself  is  not  well  known. 

The  (M$  -  0^)  plot  using  0^  as  parameter  and  the  (Ms  -  x’)  plot  using 
0W  for  the  frozen-gas  model  are  shown  in  Fig.  51  and  52,  respectively.  In 
both  of  these  figures,  the  frozen  analytical  lines  are  seen  to  be  all 
shifted  up  from  their  equilibrium  values.  The  measured  value  of  for  each 
data  point  is  identical  to  that  given  in  Fig.  49  and  is  not  repeated  in  Fig. 
51.  Since  the  disagreement  between  experiments  and  frozen-gas  analysis  is 
large,  the  discrepancy  lines  are  not  drawn  in  Fig.  51.  In  the  (M$  -  x') 
plot  of  Fig.  52,  the  frozen  gas  analytical  results  are  all  too  much  larger 
than  the  measured  experimental  values  and  the  agreement  is  poor. 

The  other  angle  that  is  worth  discussing  is  the  difference  between  x‘ 
and  x,  that  is  ( x '  -  x) ,  since  it  indicates  the  position  of  the  second 

triple  point  relative  to  the  first  triple  point.  A  [M$  -  (x'  -  x)J  plot  of 
the  experimental  results  is  shown  for  five  wedge  angles  =  42°,  37°,  20°, 
10°  and  4°  in  Fig.  53.  A  dashed  line  is  used  to  approximate  the  profile  of 

( x '  -  x)  for  each  fixed  wedge  angle.  As  SMR  just  terminates  (e.g.  at  Gh  = 

42°  and  Ms  =  1.55),  (x'  -  x)  =0,  and  w'  is  positive  indicating  that  the  two 
triple  points  merge  as  one  at  the  first  triple  point.  When  CMR  and  DMR 
begin  to  form,  the  value  of  ( x ‘  -  x)  increases  with  Mach  number  Ms.  It  then 
reaches  a  maximum  and  decreases  in  value.  In  the  positive  region  of 
( x '  -  x)  »  w’is  positive  as  shown  in  Fig.  54(a).  As  M$  increases  further, 

( x 1  -  x)  crosses  zero  and  goes  into  the  negative  region.  In  this  case  of 
(x‘  -  x  )  =  0,  w1  is  zero  as  shown  in  Fig.  54b  meaning  that  the  second 

triple  point  lies  at  a  finite  distance  away  from  the  first  triple  point  and 
both  have  the  same  trajectory  direction.  It  is  seen  in  Fig.  53  that  the 
higher  the  wedge  angle,  the  sooner  (x'  -  x)  goes  negative.  When  ( y ‘  -  x) 
goes  negative,  <*>'  is  also  negative  and  the  reflection  pattern  is  illustrated 


in  Fig.  54c.  At  higher  wedge  angle  and  Mach  number,  x  =  0,  (x1  -  x)  =  -x 
and  a)'  is  negative  as  shown  in  Fig.  54d.  Lee  and  Glass  [25]  studied  this 
last  case  analytically  for  a  perfect  gas  with  y  =  1.093  and  reported  that 
X1  =  0  may  only  be  hypothetical  and  it  does  not  exist  in  practice.  However, 
in  the  present  study,  the  experimental  results  confirmed  the  equilibrium- 
gas  analysis  that  the  case  of  x’  =  0  does  occur.  This  means  that  the  second 
triple  point  T*  attaches  to  the  wedge  surface  and  the  reflected  shock  wave  R 
strikes  the  wedge  surface  and  is  reflected  like  an  RR.  A  shadowgram  of  a 
OMR  with  x 1  =  0  taken  at  Ms  =  7.96  and  Ow  =  37°  is  shown  in  Fig.  55.  Note 
the  small  triangular  region  bounded  by  the  Mach  stem,  the  reflected  shock 
wave  and  the  wedge  surface.  The  second  reflected  shock  wave  is  lying  low  on 
the  wedge  surface  and  terminates  at  the  corner  as  an  attached  straight  bow 
shock  wave  showing  that  the  state  behind  the  second  reflected  shock  wave  is 
supersonic.  Consequently,  this  pattern  of  UMR  is  different  from  the  other 
three  patterns  shown  in  Fig.  54(a)  to  (c).  A  DMR  with  x'  =  0  has  two 
confluence  points,  two  reflected  shock  waves,  but  only  one  Mach  stem,  one 
slipstream  and  the  region  behind  the  second  reflected  shock  wave  is 
supersonic  instead  of  subsonic.  Figure  56  shows  an  interferograin  of  a  DMR 
with  x 0  at  M$  =  4.20,  ow  =  37°.  The  second  slipstream  is  clearly  seen  in 
this  picture. 


The  behaviour  of  x‘  relative  to  x  can  be  explained  by  the  following 
physical  interpretation.  For  a  given  wedge  angle,  when  SMR  just 

terminates,  the  kink  K  is  at  the  first  triple  point  T  and  (x*  -  x)  -  0. 
With  a  slight  increase  in  M  ,  CMR  occurs.  Now,  5  is  slightly  larger  than 
90°,  the  kink  K  moves  to  a  finite  distance  on  the  reflected  shock  R  away 

from  the  first  triple  point  T,  and  in  such  circumstances  (x‘  -  x)  >  0*  CMR 
continues  to  occur  with  further  increase  in  Ms ,  K  moves  further  away  from  T, 
the  value  of  5  gets  larger  and  (x*  -  x)  also  increases  slightly.  At  the 
onset  of  DMR,  the  portion  of  the  reflected  shock  R  between  T  and  K  has  to 
bend  further  down,  ie.  the  angle  6  has  to  increase,  to  have  greater 
density  change  across  the  kink.  When  the  densities  p2  >  n 2  as  shown  in  Fig. 
9(b)  become  large  enough,  DMR  occurs.  Because  the  reflected  shock  R  drops 
(6  increases)  and  the  pressure  in  regions  (4)  and  (b)  pushes  the  second 

triple  point  T*  upstream,  x*  is  getting  closer  to  x  and  (x’  -  x)  decreases. 
If  Ms  is  increased  further,  x'  will  overlap  x  and  (x'  -  x)  =  0  aga'n, 

however,  in  this  circumstance,  T  and  T'  are  at  a  distance  apart  and  not 

merged  at  a  point.  At  higher  M$ ,  the  flow  in  region  (1)  is  getting  stronger 

and  stronger,  driving  the  second  reflected  shock  R1  down  closer  to  wedge 
surface  and  the  pressure  in  regions  (4)  and  (5)  pushing  the  second  triple 
point  T‘  further  upstream,  thus  making  (x‘  -  x)  <  0.  For  certain  wedge 

angles,  it  is  possible  to  have  the  flow  in  region  (1)  be  so  strong  that  it 

presses  T*  to  the  wedge  surface  and  x’  =  0  or  (x‘  -  x)  =  -x.« 

4.5  Comparison  of  Analytical  and  Experimental  Results  of  Physic' 


Quantities  in  Regular  Reflection 

For  regular  reflection,  there  is  no  need  for  the  prediction  of  x 
because  x  =  Q  and  =  nw,  thus  only  the  comparison  of  the  angles  6  and  w 
witn  ow  as  a  parameter  for  frozen  and  equilibrium-gas  analyses  is  necessary 
In  fact,  a  comparison  of  6  is  essentially  equivalent  to  a  comparison  of  w 
since  the  two  are  related  by 


A  +  (O1 


(4.9) 


as  seen  from  Fig.  37.  Comparisons  of  both  angles  are  shown  here.  Table  5 
presents  the  measured  and  analytical  results  of  ft  and  <j'  in  RR. 

A  plot  of  (M^  -  ft )  plane  with  the  experimental  results  is  shown  in  Fig. 
57.  The  solid  fines  are  the  calculated  analytical  g  for  equilibrium  SF6. 
The  dashed  line  is  tne  transition  boundary  between  RRand  MR  below  which  no 
analytical  solution  of  6  in  RR  exists.  The  angle  ft  is  seen  as  an  increasing 
function  with  increasing  Mach  number  and  wedge  angle,  apart  from  the  low 
Mach  number  range  at  gw  =  45°.  Except  for  a  bit  of  scattering  of  data  point 
at  ow  =  40°  and  Ms  about  4.1,  the  agreement  between  experiment  and 
equilibrium-gas  analysis  is  very  good.  It  is  seen  that  for  ow  =  45°,  as  the 
Mach  number  decreases  from  M$  =  2.38,  the  experimental  ft  drop  more  gradually 
than  the  analytical  values.  There  are  four  RR  persistent  data  points  shown 
in  the  range  of  1.5  <  M$  <  1.7  of  the  same  wedge  angle.  Note  that  these 
persistent  RR  points  exist  in  the  range  where  no  two-shock  theory  solution 
can  be  obtained. 

Figure  58  shows  a  plot  of  the  experimental  results  in  the  (Ms  -  </) 
plane  with  the  equilibrium  analytical  lines  drawn  in  solid.  The  dashed 
line  is  the  transition  boundary  of  RR  and  MR.  On  the  contrary,  here  the 
angle  u‘  is  a  decreasing  function  with  increasing  Mach  number  and  wedge 
angle,  besides  the  low  Mach  number  range  for  <~>w  =  45°  where  an  increase  in 
(*)'  is  observed.  The  same  good  agreement  is  found  here  between  experiments 
and  equilibrium  gas  analytical  results.  At  higher  Mach  numbers,  the 
analytical  lines  of  w'  are  getting  closer  and  begin  to  level  off.  It  is 
expected  that  with  ever  increasing  Mach  number  or  raising  the  wedge  angle 
gw  closer  to  90°,  the  values  of  ft  and  <o'  will  approach  their  limits,  and  the 
reflected  shock  will  be  seen  lying  closer  and  closer  to  the  wedge  surface. 
For  M$-* *> ,  the  pressure  in  region  (1)  becomes  so  large  that  it  forces  the 
reflected  shock  R  to  lie  on  the  surface  of  the  wedge  in  steady  flow  as  well 
as  nonstat  ionary  flow,  thus  u>'*0.  This  illustrates  the  fact  that  the 
reflection  angle  V  is  always  positive  and  can  be  zero  for  RR,  whereas  to'  is 
positive  for  all  gases  and  can  be  zero  or  even  negative  for  some  low  y  gases 
in  the  case  of  MR  [25]. 

Consider  the  transition  boundary  map  for  frozen  SF6  in  Fig.  32(a).  All 
the  experimental  RR  results  lie  well  inside  the  MR  domain.  Therefore  no 
two-shock  theory  solution  for  frozen-gas  case  can  be  obtained  and  no 
agreement  between  experiments  and  the  frozen-gas  analysis  is  seen. 

It  can  be  observed  from  both  Figs.  57  and  58  that  the  behaviour  of  the 
analytical  ft  and  ,a‘  for  r  =  45°  at  M$  <  1.4  do  not  have  the  same  trends  as 
in  the  high  Mg  range.  Refer  to  Fig.  3  for  the  flow  system  and  consider  Ow 
is  fixed  at  45°  ($„  also).  The  flow  deflection  angle  Q,  and  the  flow 
velocity  UL  in  region  (1)  are  both  increasing  functions  of  II  or  M-.  If  M$ 
and  <t>0  are  given,  Oj  and  U,  can  be  determined.  Now,  the  reflected  shock  R 
needs  to  be  positioned  so  that  the  flow  IL  is  deflected  by  an  angle  02  =  0, 
across  the  reflected  shock  R  and  is  parallel  to  the  wall.  At  gw  =  45°  ana 
Ms  <  1.4,  increasing  M$  will  not  be  able  to  set  up  Uj  large  enough  to  meet 

the  boundary  condition  if  4>0  remains  constant.  To  obtain  02,  one  can  either 
fix  <t> j  and  change  or  fix  Uj  and  change  $[.  However,  with  M$  and  <t>0  being 
given,  Uj  is  fixed  and  so  the  angle  <t>j  has  to  be  increased.  R  becomes  more 
normal  to  Ui  and  can  now  be  achieved.  As  the  value  of  <t>,  increases, 
physically,  R  propagates  further  away  from  the  wedge  surface,  thus  decreases 
ft  and  increases  <,>' . 


At  M  >  2.2  for  the  same  Ow  =  45°,  Uj  increases  at  a  much  faster  rate 
than  it  does  at  low  M$.  Uj  is  getting  large  and  R  needs  to  be  more  tangent 
to  it.  As  the  wave  angle  4»1  decreases,  R  lies  closer  to  the  wall  and  as  a 
result  increases  6  and  decreases  o>* . 


4.6  Pressure  Behind  Second  Triple  Point 

Three-shock  theory  has  been  used  to  analyze  flows  around  the  vicinity 
of  the  first  triple  point  and  for  the  prediction  of  the  transition  lines 
between  RR  and  MR.  Ben-Dor  [11]  was  the  first  to  analyze  the  second 
triple-point  system  and  set  up  an  analytical  formulation  using  three-shock 
theory  and  by  applying  an  analogy  with  the  first  triple-point  system  as 
shown  in  Fig.  59.  However,  no  known  publication  has  ever  shown  any 
numerical  results  for  the  physical  quantities  at  the  second  triple-point 
system.  It  is  well-known  from  the  transition  maps  that  the  most  probable 
type  of  Mach  reflection  which  will  occur  at  high  M.  when  RR  terminates  is  a 
DMR,  since  the  RR-^-OMR  boundary  is  the  dominant  line,  especially  for  y  < 
1.4,  that  covers  a  large  range  of  Mach  numbers  in  the  RR<->-MR  region.  It  is 
known  from  experiments  that  the  flow  regions  behind  the  first  triple  point  T 
are  uniform,  thus  the  assumption  of  two  uniform  states  behind  T  in  the 
analysis  of  the  first  triple-point  system  is  well  justified.  However,  the 
flow  regions  behind  the  second  triple  point  T*  are  not  uniform  owing  to  the 
fact  that  the  second  Mach  stem  M'  is  nonuniform  (for  example,  see  Figures  26 
and  28).  In  order  to  analyze  the  second  triple-point  system  numerically, 
regions  (4)  and  (5)  behind  T*  have  to  be  assumed  as  two  uniform  states.  Of 
course,  if  a  numerical  solution  of  the  Euler  equations  is  obtained  for  the 
entire  flow,  such  an  assumption  is  not  necessary.  Since  the  growth  of 
second  Mach  stem  M'  contradicts  the  assumption  of  uniform  states  behind  T1, 
tnerefore  the  results  can  only  be  applied  in  the  vicinity  of  the  second 
triple  point.  A  method  of  calculating  the  second  triple  point  flow  fields 
is  given  in  Appendix  A. 

Since  RR  occurs  at  higher  wedge  angles  than  MR  for  any  given  flow  Mach 
number  in  region  (0),  after  the  passage  of  incident  I  ana  reflected  R 
shocks,  P2  acquires  a  high  pressure.  In  MR  the  highest  pressure  would  be 
achieved  by  a  DMR  since  the  flow  in  region  (2)  becomes  supersonic  relative 
to  the  second  triple  point  T1  and  a  shock  wave  (second  Mach  stem  M')  is 
required  to  turn  the  flow  subsonical ly .  Behind  the  second  triple  point  T‘, 
the  pressures  are  raised  to  P4  and  P5.  Hence,  it  is  only  necessary  to 

compare  the  highest  pressure  attained  by  RR  and  DMR,  and  determine  which 

type  of  reflection  would  give  the  highest  pressure. 

The  highest  pressure  in  RR  occurs  behind  the  reflected  bow  shock  R  at 
point  E  immediately  behind  the  reflection  point  P  as  illustrated  in  Fig. 
60(a).  In  DMR  the  highest  pressure  occurs  in  region  F  immediately  behind 
the  second  Mach  stem  M'  and  the  primary  slipstream  S  as  illustrated  in  Fig. 
60(b).  Two-shock  theory  can  be  used  to  calculate  pressure  P2  in  the 
vicinity  of  the  reflection  point  P  in  RR,  thus  the  highest  pressure  at  E 
takes  the  value  of  P2.  However,  three-shock  theory  can  only  be  used  to 
calculate  pressures  around  the  two  triple  points  and  not  elsewhere.  Due  to 
this  limitation,  only  the  pressure  behind  the  second  triple  point  in  DMR, 
shown  as  point  G  in  Fig.  60(b),  can  be  compared  with  the  highest  pressure  in 
RR.  Although  such  comparison  is  not  ideal,  it  does  provide  a  general  trend 
and  an  approximation  of  the  highest  pressure  in  these  two  types  of 

reflection.  The  maximum  pressure  in  DMR  occurs  on  the  wedge  surface  in 


rigion  F  as  indicated  by  the  numerical  results  of  Glaz  et  al  [63],  Booen  and 
Needham  [64],  and  Fry  et  al  L 65 J -  The  location  of  maximum  pressure  on  the 
wedge  surface  is  also  consistent  with  the  experimental  density  contours  from 
the  present  work. 

To  study  the  behaviour  of  the  pressure  behind  T'  in  RR  and  DMR  a  series 
of  computational  runs  at  a  fixed  Mach  number  of  M  =  6  (M$  =  6  is  chosen 
arbitrarily  as  it  cuts  across  the  RR^-*DMR  transition  boundary  line;  other  Ms 
can  also  be  used  instead)  over  a  wide  range  of  Ow  was  done  in  perfect  Ar  to 
avoid  any  real-gas  effects.  The  results  are  plotted  in  Fig.  61.  When 
9w+90°,  P2  approaches  the  value  given  by  the  relation  for  a  normal  shock 
wave  reflection  P2,  =  (a  +  2  -  P01)/(l  +  a  PQ1)  where  a  =  (y+1)/(y-1)  and 
P01  is  the  strength  of  the  incident  shock  I.  As  ^  is  decreased  from  90°, 
P2  decreases  and  reaches  a  minimum  value  before  it  increases  to  much  higher 
values  near  the  transition  angle.  Such  behaviour  can  be  explained  as 
follows:  since  the  incident  shock  Mach  number  is  fixed,  the  velocity 

component  normal  to  the  incident  shock  1  (with  respect  to  the  reflection 
point  P)  and  the  pressure  P,  are  constant  regardless  of  G^.  As  0W  is 
decreased,  the  flow  across  tne  incident  shock  is  deflected  to  a  lesser 
extent  and  U1  gets  smaller.  To  satisfy  the  boundary  condition  that  the  flow 
behind  the  reflected  shock  must  be  parallel  to  the  wall,  ie.  o2  =  O,,  the 
wave  angle  <t> L  increases  and  causes  the  flow  U  t  to  become  more  normal  to  the 
reflected  shock.  The  pressure  P2  depends  on  ll'l  -  the  normal  component  of  Uj^ 
to  the  reflected  shock.  Decreasing  Uj  will  reduce  Un1,  whereas  increasing 
the  wave  angle  <t>,  will  increase  the  value  of  U1] .  With  U,  decreasing  and  ^ 
increasing,  a  minimum  value  of  Ur|  and  P2  would  be  obtained. 

The  highest  pressure  of  P2  in  the  RR  domain  is  obtained  at  the 
transition  angle  (point  X  in  Fig?  61).  Note  that  this  P2  is  even  higher 
than  the  head-on  collision  value.  With  a  slight  decrease  of  at  the 
transition  angle,  DMR  begins  to  form.  The  pressure  P^  and  P5  behind  T' 
right  after  transition  (point  Y  in  Fig.  61)  is  slightly  higher  than  P2  in  RR 
just  before  transition.  However,  that  is  not  always  the  case.  The 
pressure  P2  in  RR  and  P4  and  Pb  at  point  G  in  DMR  (see  Fig.  60)  at  the 
transition  angles  for  4  <  M^  <  10  are  presented  in  Table  6.  It  is  shown 
that  at  lower  M$,  P2  in  RR  is  larger  than  P4  and  P5  in  DMR;  at  higher  M$,  P4 
and  Pb  in  DMR  become  greater  than  P2  in  RR?  Once  DMR  is  formed,  P4  and  P5 
are  decreasing  functions  of  ,  as  shown  in  Fig.  61.  It  is  surprising  that 
the  flow  problem  is  non-linear,  yet  P4  and  P5  in  DMR  behave  almost  linearly 
with  v3w.  The  bow  shock  which  terminates  on  the  shock-tube  wall  near  the 
wedge  corner  C  moves  outward  at  a  very  much  slower  rate  than  the  incident 
shock  I.  By  assuming  that  the  bow  shock  wave  near  the  corner  C  is 
stationary,  the  stagnation  pressure  at  the  corner  C  can  be  computed  and  the 
results  are  also  presented  in  Table  6.  Although  the  stagnation  pressure  at 
corner  C  is  high,  its  effects  are  localized.  However,  the  pressure  P2  in 
RR,  and  the  pressures  P4  and  Pb  in  DMR  propagate  up  the  wedge  along  with  the 
incident  shock  wave.  Therefore,  the  compressive  pressure  they  exert  is  not 
a  local  effect  and  becomes  of  major  interest.  There  is  a  certain  range  of 
pressure  in  which  a  required  pressure  can  be  obtained  with  more  than  one  0^ 
at  point  E  in  RR  or  at  point  G  in  DMR,  e.g.,  point  A  and  B  in  the  RR  domain, 
and  point  C  in  the  DMR  domain  all  give  the  same  pressure  for  three  different 

V 

The  following  conclusions  can  be  made.  Although  at  lower  M$  the 
pressures  P4  and  P,  in  DMR  are  lower  than  P2  in  RR,  the  highest  pressure  in 
region  F  is  higher'  than  P4  and  Pb  at  point  G  as  indicated  by  experimental 
and  numerical  results.  Therefore,  the  highest  pressure  in  DMR  can  be  higher 


than  the  highest  pressure  in  RR.  At  higher  M  ,  P4  and  Pb  at  point  G  in 
DMR  are  already  higher  than  P.  in  RR,  thus  the  highest  pressure  in  region  F 
in  OMR  is  even  higher  than  tne  h.ghest  pressure  in  RR.  It  can  be  stated 
with  confidence  that  with  a  strong  incident  shock  wave  (Ms  >  5  for  y  =  5/3) 
the  highest  pressure  is  achievable  by  a  DMR.  With  weaker  incident  shock 
wave  (M$  <  5  for  y  =  5/3)  it  can  only  be  surmised  that  the  highest  pressure 
can  be  attained  by  a  DMR  as  well.  Data  from  complete  numerical  simulations 
will  provide  better  answers.  Due  to  the  uncertainty  of  how  good  the 
analysis  of  the  second  triple-point  system  predicts  the  flow  properties 


behind  T‘,  comparisons  of  the  experimental  pressure  history  of  DMR  in  air 
[26]  with  the  numerical  simulation  by  Glaz  et  al  [63]  and  with  the  numerical 
results  of  this  work  are  given  by  Hu  and  Shirouzu  [31].  From  a  practical 
point  of  view  more  field  testing  and  shock  tube  experiments  should  be 
directed  at  studying  the  transition  regime  between  RR  and  DMR  to  verify  the 
numerical  results  obtained  from  the  computer  program  developed  during  the 
course  of  this  work. 

4 . 7  Application  of  Reflections  in  Pseudo-Stati onary  Flow  to  the 

Interaction  of  Spherical  Blast  Waves  with  Planar  SuTTace 

With  the  knowledge  of  blast  wave  flows  and  obi i que-shock -wave 
reflections  gained  in  recent  years,  it  is  now  possible  to  trace  the  path  of 
a  spherical  blast  wave  that  is  detonated  in  free  air  on  a  transition  map  in 
the  (M$  -  P>w)  plane.  It  is  true  that  in  an  explosion,  the  flow  field  near 
the  charge  is  very  complicated  and  very  difficult  to  analyze,  but  when  the 
blast  wave  front  outruns  the  fireball  and  starts  to  interact  with  the  ground 
surface,  the  flow  can  be  reasonably  treated  as  two-dimensional  at  any 
instant  of  time.  Hence,  the  analysis  used  here  for  pseudo-stationary  flow 
can  be  applied  to  the  interaction  of  a  blast  wave  with  a  planar  surface.  It 
becomes  possible  to  predict  the  type  of  reflection,  the  incident  shock  Mach 
number  and  the  corresponding  wedge  angle  that  a  target  at  a  known  distance 
away  from  the  explosion  may  encounter.  This  information  is  of  great  value 
to  researchers  and  experimenters  in  the  field. 

The  flow  field  associated  with  a  free-air  explosion  is  shown  in  Fig. 
62.  To  demonstrate  how  the  study  of  oblique-shock-wave  reflection  can  be 
applied  to  the  interaction  of  a  spherical  blast  wave  with  a  planar  surface, 
a  typical  charge  of  TNT  detonated  in  free  air  off  a  perfect  reflecting 
planar  surface  in  a  standard  atmosphere  of  760  torr  (101.33  kPa)  and  288  K 
is  considered.  Sadek  and  Gottlieb  [66]  obtained  an  equation  to  describe  the 
spherical  blast-wave  front  in  a  surface  explosion  for  a  1  kg  TNT.  Their 
results  were  reduced  to  the  case  of  a  1  kg  TNT  hei  ght-of-burst  (HOB) 
explosion  using  Sach's  scaling  laws  with  a  geometric  factor  of  2  for  a 
perfect  reflecting  surface.  The  scaled  equation  would  then  give  the  shock 
radius  r  from  the  charge  C  as  a  function  of  the  shock  arrival  time  t$ .  For 
a  given  HOB  h,  the  ground  range  x  and  the  corresponding  wedge  angle  are 
related  geometrical  ly  with  the  radial  distance  r  in  the  case  of  RR. 
However,  in  the  case  of  MR,  the  computation  of  the  ground  range  x  and  the 
cor  respond i ng  wedge  angle  o  are  complicated  by  the  introduction  of  the 
first-triple-point  trajectory.  Since  the  Mach  stem  M  usually  occurs  bulged 
out  in  real  explosions  or  large-scale  field-tests,  the  ground  range  x  is 
simply  estimated  by  assuming  that  the  Mach  stem  M  is  an  arc  of  a  circle  with 
the  radius  being  the  distance  between  zero  ground  range  and  the  first  triple 
point  T,  which  is  shown  as  r^  in  Fig.  62(b).  Details  of  the  method  of 
calculating  the  blast-wave-trajectory  path  are  given  in  Appendix  A.  Since 


the  spheri cal -blast-wave  flow  can  be  treated  as  two-dimensional  at  any 
instant  of  time,  the  merging  of  the  two  triple  points  at  the  RR<->DMR 
transition  boundary  applies,  and  there  exists  a  blast-wave-trajectory  path 
which  crosses  the  point  where  all  three  transition  lines  meet  on  the 
transition  map. 

The  computed  results  of  the  blast-wave-front  trajectory  path  for 
various  HOB  cases  are  presented  graphically  in  Fig.  63  on  the  transition  map 
in  the  (M  -  Qw)  plane.  Contours  of  ground  range  x  and  height  of  the  first 

triple  point  Yy  are  also  plotted.  Height  of  burst  h  and  ground  range  x  in 

the  figure  are  scaled  to  the  weight  of  the  charge  W  and  the  ambient  pressure 
P0  by  introducing  a  non-dimensi onal i zed  scaling  factor  of  (W/W0  •  Po/P)-1/3 
where  W0  and  P0  are  taken  as  1  kg  and  1  atm,  respectively.  Height  of  the 
first  triple  point  Yy  is  also  shown  as  a  fraction  of  the  HOB  h. 

Assume  that  a  charge  of  1  kg  TNT  is  detonated  at  a  HOB  of  1  m  with  the 

ambient  pressure  being  1  atm.  When  the  blast  wave  front  just  hits  the 
planar  surface  directly  under  the  charge,  x  =  0  m,  Gw  =  90°  and  =  4.4. 
As  the  blast  wave  front  propagates  outward  further,  the  incident  shock  wave 
collides  with  and  is  reflected  off  the  planar  surface  resulting  in  regular 
reflection  HR,  which  occurs  up  to  a  ground  range  of  x  =  0.82  m  and  Ms 
diminishes  to  3.5.  When  DMR  just  occurs,  the  height  of  the  first  triple 

point  Yy  is  zero  at  the  transition  line  and  DMR  continues  in  the  range  of 

0.82  m  <  x  <  1.19  m.  At  the  transition  point  from  DMR  to  CMR,  M$  decreases 

to  2.9,  Yy  increases  to  0.022  m,  and  CMR  occurs  in  the  ground  range  of  1.19 
m  <  x  <  1.9  m.  At  the  termination  of  CMR,  M$  drops  to  2.3,  Yy  grows  to 
almost  0.1  m,  and  SMR  begins  to  form  at  a  ground  range  x  >  1.9  m. 

With  a  given  charge  weight  W,  ambient  pressure  P,  HOB  h  and  ground 
range  x,  one  can  read  off  from  Fig.  63,  M$ ,  Qw,  Yy  and  the  type  of 
reflection  that  is  going  to  occur.  With  M$  and  0W,  the  flow  properties 
behind  the  blast  wave  front  at  the  ground  range  x  can  be  predicted  using 
the  numerical  results  presented  in  tabular  and  graphical  forms  by  Hu  and 
Shirouzu  [31].  These  results  are  essential  in  predicting  the  dynamic 

response  of  structures  in  the  vicinity  of  a  blast  wave.  For  a  different 
explosive  charge  and  height  of  burst,  a  different  blast-wave-front  path  has 
to  be  constructed  on  the  transition  map  in  order  to  predict  what  type  of 
reflection  would  occur  for  a  certain  ground  range.  Figure  63  also  shows 
the  effect  of  changing  the  height  of  burst  and  charge  weight  on  the  blast 
wave  front  trajectory  path.  The  trajectories  of  12  different  HOB  cases  are 
calculated  and  plotted  on  the  same  figure  to  show  the  effect  of  changing  the 
HOB. 

There  are  several  interesting  points  that  can  be  made  about  Fig.  63. 
First,  increasing  the  HOB  for  a  given  charge  will  decrease  the  incident 

shock  wave  Mach  number  M$  of  the  blast  wave  as  it  interacts  with  the  planar 
surface.  Second,  let  Ax  be  defined  as  the  difference  between  the  ground 
ranges  when  a  specific  type  of  reflection  begins  and  terminates.  For 

example,  on  the  scaled  line  h  =  1.0  m,  x  increases  from  0  m  for  0^  =  90°  to 
0.82  m  when  Gws50°  at  the  RR*->DMR  line  and  Ax  =  0.82  m.  In  RR  region.  Ax 
increases  with  HOB.  For  example,  at  a  scaled  h  of  0.7  m.  Ax  =  0.58  m;  at  h 
=  1.00  m.  Ax  =  0.82  m;  at  h  =  2.00  m,  Ax  =  1.7  m  and  at  h  =  3.00  m.  Ax  =  2.8 

m.  However,  in  both  DMR  and  CMR  regions,  there  is  a  HOB  h  which  gives  an 

optimum  Ax.  This  is  because  at  small  h,  Ax  increases  with  HOB.  As  h 
becomes  large,  the  domains  of  both  DMR  and  CMR  diminish  and  so  their  Ax  also 


decrease.  For  instance,  in  the  DMR  region,  at  a  scaled  h  of  l).  5  m.  Ax  = 
(0.83-0.41)  =  0.42  m;  at  h  =  0.80  m.  Ax  =  (1.10-0.67)  =  0.43  m;  at  h  =  0.9 
m,  Ax  =  (1.16-0.  76)  =  0.40m;  at  h  =  1.00  m.  Ax  =  (1.19-0.82)  =  0.37  m  and  at 
n  =  1.26,  Ax  =  (1.28-1.01)  =  0.27  m.  Due  to  the  increasing  and  decreasing 
behaviour  of  Ax,  a  maximum  Ax  exists  in  the  DMR  domain  and  this  also  applies 
to  the  CMR  domain.  In  the  SMR  region,  Ax  is  an  increasing  function  of  HOB 
since  the  SMR  domain  enlarges  at  low  M  .  Third,  theoret i ca 1 1 y ,  the 

interaction  of  a  blast  wave  front  with  tne  ground  may  result  in  three 
possible  series  of  reflections:  (1)  RR-*DMR->CMR-*SMR,  (2)  RR->CMR>5MR,  or  (3) 
RR+SMR.  However,  experiments  have  shown  that  the  CMR  to  DMR  transition  line 
approaches  the  SMR  to  CMR  transition  line  and  merges  at  a  single  point  on 
the  RR  termination  boundary  (see  Sec.  4.3)  as  shown  by  the  dotted  line  in 
Fig.  63.  Therefore,  only  two  possible  series  of  reflection  can  occur:  (1) 
RR>  DMR+  CMR  -»SMR  or  (2)  RR+SMR.  The  second  possibility  occurs  if  the 

incident  shock  Mach  number  of  the  blast  wave  Ms  <  1.62  when  RR  terminates. 
Fourth,  from  Fig.  63  it  can  be  seen  that  owing  to  the  steepness  of  the 
RR*->SMR  transition  line  compared  to  the  blast-wave  paths  in  the  (M$  -  C^) 
plane,  RR  cannot  reoccur  when  M  >1.  Lastly,  increasing  the  charge  weight  W 
has  the  same  effect  as  decreasing  the  HOB,  whereas  an  increase  in  the 

ambient  pressure  P  has  results  similar  to  increasing  the  HOB. 

The  method  used  in  this  work  to  study  the  interaction  of  a  blast  wave 
with  a  planar  surface  for  the  explosive  TNT  can  be  used  for  other  types  of 
explosives  in  predicting  the  blast-wave-front-trajectory  path  on  the 

transition  map,  and  in  calculating  the  correspond1! ng  ground  range  and  height 
of  the  first  triple  point  for  the  HOB  case.  For  a  spherical  flow,  the 

growth  of  the  reflected  spherical  shock  cannot  be  determined  from  a 

two-dimensional  analysis.  Consequently,  the  second  triple  point  T1 
trajectory  is  indeterminate  using  this  analysis.  It  would  be  important  to 
verify  these  results  experimentally  in  the  labortory  or  in  field  tests. 

5.0  CONCLUSIONS 

An  experimental  and  numerical  investigation  was  made  of 
pseudo-stationary  obi i que-shock -wave  reflections  in  polyatomic  SF6  gas.  The 
domains  and  transition  boundaries  between  the  various  types  of  reflection 
(RR,  SMR,  CMR  and  DMR)  were  established  in  the  (M$  -  0W)  and  (M$  -  R^) 
planes  for  both  frozen  and  vibrational  equilibrium  SF&.  The  transition 

boundaries  predicted  on  the  basis  of  an  equilibrium  flow  agreed  well  with 

the  experimental  results.  The  new  criterion  [lb]  6  >  90°  was  verified  by 
the  present  work  that  it  is  an  additional  necessary  condition  for  the 
transition  from  SMR  to  CMR.  However,  RR  persists  beyond  the  boundary  line 
determined  by  the  detachment  criterion.  A  start  has  been  made  to  explain 
the  persistence  of  RR  in  terms  of  the  induced  viscous  boundary  layer  on  the 
wedge  surface  [16,23].  A  continuation  of  the  investigation  of  the 
persistence  of  RR  with  its  boundary-layer-di spl acement  effect  is  presently 
under  way  [43]. 

Experimental  results  show  that  the  CMR  to  DMR  transition  line 
approaches  the  SMR  to  CMR  transition  line  near  the  RR  termination  boundary 
and  they  eventually  merge  at  the  RR  boundary.  The  fact  that  the  ratio  7*/L* 
approaches  zero  as  the  wedye  angle  comes  near  to  the  RR  termination  boundary 
was  also  verified  in  this  experimental  study.  A  better  analytical  method  to 
determine  the  position  of  the  kink  or  the  second  triple  point  near  the  RR 
boundary  is  required  in  order  to  predict  accurately  the  behaviour  of  the 
merging  of  the  three  transition  lines  at  a  sinyle  point  on  the  RR 
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termination  line.  It  is  also  clear  from  the  present  study  that  improved 
analytical  criteria  must  be  found  for  the  SMR->CMR->OMR  transitions  in  order 
to  have  better  agreement  between  analysis  and  the  experimental  data.  Since 
it  is  difficult  to  distinguish  CMR  from  DMR  in  an  interferogram,  a 
shadowgram  is  required  in  the  transition  region  of  CMR  and  DMR  to  provide 
additional  evidence. 

Several  infinite-fringe  interferograms  of  the  density  fields  were 
evaluated  and  their  wall-density  distribution  were  plotted.  It  has  been 
shown  that  the  use  of  higher  initial  pressure  enhanced  the  details  and 
accuracy  of  the  isopycnics  significantly.  The  present  experimental  results 
provide  invaluable  data  for  computational  fluid  dynamicists  in  testing  and 
improving  the  numerical  codes  that  predict  various  flow  properties  in 
nonstationary  oblique  shock-wave  reflections. 

The  comparison  of  experimental  results  in  Ar,  air,  C02  and  SF6  at  GL  = 
20°  and  M$  =  3.4±0.2  substantiate  the  claims  of  Lee  and  Glass  [25]  that  the 
transition  boundaries  are  shifted  downwards  and  the  reflection  angle  u>' 
decreases  with  decreasing  y.  There  are  higher  density  jumps  across  the  Mach 
stem,  slipstream  and  the  bow  shock  or  lambda  shocks  at  the  wedge  corner  for 
smaller  y.  The  shock-wave  boundary-layer  interaction  (lambda- shock 
interaction),  rolling  up  of  the  slipstream  and  curving  of  the  Mach  stem  are 
all  prominent  with  a  lower  y  gas,  because  a  stronger  shock  can  be  produced 
by  a  lower  y  gas  as  given  by  the  strong  shock  wave  limit  p/ p0>(  y+1 )/(  y-1 ) 
and  the  density  ratio  is  large  as  y  >  1. 

The  large  differences  in  the  analytical  physical  quantities  between 
frozen  and  equilibrium  SF6  permit  a  definitive  conclusion  regarding  which 
model  is  valid.  Comparisons  of  the  fundamental  angles  6  and  u>' ,  and  rr  (  Cfc 
-  Ow  +  x' )  in  Mach  reflection  using  as  the  parameter,  and  the 
comparisons  of  6  and  w'  in  RR  using  as  the  parameter,  justified  the 
validity  of  the  equilibrium  SF6  analysis.  The  measured  angles  6,  a>' ,  y  and 
X1  in  MR  are  also  compared  with  the  analytical  results  using  ^  as  the 
parameter.  Nonetheless,  the  agreement  with  the  equil ibrium-gas  analysis 
using  Ow  is  secondary  to  that  of  using  0^,  due  to  the  inadequate  prediction 
of  x»  by  assuming  that  the  Mach  stem  is  straight  and  perpendicul ar  to  the 

wedge  surface.  An  improved  method  of  predicting  x  which  takes  into  account 

the  Mach-stem  curvature  and  the  interaction  of  the  bow  shock  wave  with  the 
reflected  shock  wave  is  necessary.  It  has  been  verified  that  the  reflection 

angle  <*>'  is  positive  in  RR  and  can  be  positive  or  negative  in  MR  L25J.  The 

angle  between  the  first  and  second-triple-point  trajectories  ( y  -  \‘ )  was 
discussed  and  a  DMR  with  x'  =  0  was  obtained  experimentally. 

The  second  triple-point  system  was  analyzed  numerically  for  argon  (y  = 
5/3)  at  Ms  =  6  over  a  large  range  of  G^.  When  the  transition  of  RR  to  DMR 

occurs,  the  pressure  behind  the  second  triple  point  T’  (regions  4  and  5)  in 

DMR  is  almost  equal  to  the  pressure  behind  the  reflection  point  P  (region  2) 
in  RR.  Numerical  results  show  that  for  a  given  Mach  number,  the  highest 
pressure  is  achieved  by  means  of  a  DMR,  and  not  by  a  RR.  This  may  be 

because  the  flow  in  region  (0)  has  to  experience  two  or  three  shock  waves 

before  it  can  reach  region  (4)  or  (5).  Pressure  histories  as  functions  of 
Ms  and  ©w  obtained  by  numerical  simulation  at  the  RR<->-DMR  transition 
boundary  are  required  to  substantiate  the  present  experimental  and 
analytical  results. 
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A  typical  explosive  charge  (TNT)  detonated  at  various  hei ghts-of-burst 
above  a  perfect  reflecting  planar  surface  in  air  was  considered.  The 
variations  of  the  incident  shock  Mach  number  of  the  spherical -blast-wave 
front  as  it  decays  and  the  corresponding  wedge  angle  were  plotted  on  a 
transition  map  in  the  (Ms  -  %}  plane.  It  shows  that  all  four  types  of 
shock-wave  reflection  can  occur  in  a  free  air  explosion.  However,  if  the 
height  of  burst  is  increased  past  a  certain  limit,  only  two  types  of 
shock-wave  reflection,  RR  and  SMR ,  can  occur. 


The  collision  of  a  Mach  stem  from  a  spherical  explosion  with  various 
types  of  buildings  is  of  practical  importance.  This  is  presently  being 
simulated  at  UTIAS  L62]  by  investigating  the  collisions  of  the  three  types 
of  Mach  reflection  with  a  90°  ramp  fixed  to  the  wedge  surface. 
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Table  1 

Dale  Constants  for  Several  Gases  [26] 

Gas 

K  (m3/kg) 

air 

2.274  x  10-4  (X  -  5896  A) 

Ar 

1.574  x  10~4  (X  -  6943  A) 

Ar 

1.629  x  10~4  (X  *  3472  A) 

go2 

2.266  x  10"4  (X  =*  5896  A) 

H2 

7.870  x  10-5  (A  «  5896  A) 

He 

2.017  x  10~4  (X  -  5896  A)  j 

N2 

2.376  x  10~4  (A  -  6943  A) 

n2 

2.460  x  lO'4  (A  =  3472  A) 

N 

3.280  x  10'4  (A  =  6943  A) 

-4  o 

N 

3.310  x  10  (A  -  3472  A  ) 

—  4  © 

S?6 

1.182  x  10  (A  *  5896  A ) 

—  4  ° 

so2 

2.344  x  10  (X  «=  5896  A  ) 

In it ial 

Table  2  (continued) 

Conditions  for  SF6  Experiments 

EXP 

6w 

(deg) 

Ms 

Type 

Po 

(torr) 

(kPa) 

(psi) 

P3 

(kPa) 

To 

(°c) 

Dr  iver 

29 

42 

2.74 

- 

14.00 

1.9 

20.0 

137.9 

24.2 

C02 

31 

42 

2.82 

RR 

23.85 

3.2 

35.0 

241.3 

24.3 

C02 

108 

42 

2.84 

DKR 

42.90 

5.7 

63.5 

437.8 

24.0 

C02 

107 

42 

3.00 

DMR 

31.97 

4.3 

64.5 

444.7 

24.4 

C02 

76 

42 

3.01 

DMR 

70.90 

9.5 

64.5 

444.7 

24.2 

C02 

25 

42 

3.02 

RR 

70.90 

9.5 

137.5 

948.0 

25.8 

C02 

28 

42 

3.10 

- 

109.00 

14.5 

145.0 

999.7 

24.7 

C02 

77 

42 

3.10 

- 

49.00 

6.5 

110.0 

758.4 

23.8 

C02 

27 

42 

3.12 

- 

78.70 

10.5 

125.0 

861.8 

24.8 

C02 

106 

42 

3.14 

DMR 

25.86 

3.4 

62.0 

427.5 

24.0 

C02 

74 

42 

3.14 

DMR 

50.00 

6.7 

115.0 

792.9 

23.2 

C02 

24 

42 

3.16 

RR 

55.80 

7.4 

136.0 

937.7 

25.8 

C02 

26 

42 

3.20 

- 

56.00 

7.5 

132.0 

910.1 

25.8 

C02 

105 

42 

3.21 

DMR 

23.30 

3.1 

62.0 

427.5 

23.9 

C02 

23 

42 

3.21 

RR 

50.00 

6.7 

136.0 

937.7 

25.2 

C02 

22 

42 

3.26 

RR 

49.00 

6.5 

136.0 

937.7 

25.6 

C02 

78 

42 

3.26 

DMR 

49.00 

6.5 

136.0 

937.7 

24.2 

C02 

75 

42 

3.27 

DMR 

49.00 

6.5 

135.0 

930.8 

23.5 

C02 

104 

42 

3.43 

DMR 

16.73 

2.2 

63.5 

437.8 

24.2 

C02 

79 

42 

3.44 

- 

37.00 

4.9 

136.5 

941.1 

22.8 

C02 

85 

42 

3.45 

DMR 

17.35 

2.3 

60.5 

417.1 

23.7 

C02 

,r*VT*i 


Table  2(continued) 


Initial  Conditions  for  SF6  Experiments 


EXP 

6w 

(deg) 

Ms 

Type 

Po 

(torr) 

(kPa) 

(psi) 

P3 

(kPa) 

To 

(°C) 

Driver 

86 

42 

3.45 

DMR 

15.95 

2.1 

61.5 

424.0 

23.4 

C02 

100 

42 

3.47 

RR 

15.95 

2.1 

61.5 

424.0 

22.8 

C02 

84 

42 

3.48 

DMR 

15.95 

2.1 

62.5 

430.9 

23.3 

C02 

81 

42 

3.57 

DMR 

23.87 

3.2 

105.0 

723.9 

24.7 

C02 

32 

42 

3.63 

RR 

11.70 

1.6 

59.5 

410.2 

24.5 

C02 

83 

42 

3.64 

RR 

11.65 

1.6 

62.5 

430.9 

26.0 

C02 

103 

42 

3.65 

DMR 

12.10 

1.6 

64.5 

444.7 

23.9 

C02 

80 

42 

3.66 

DMR 

27.00 

3.6 

136.5 

941.1 

22.9 

C02 

102 

42 

3.66 

DMR 

12.18 

1.6 

65.5 

451.6 

22.0 

C02 

109 

42 

3.67 

RR 

11.65 

1.6 

63.5 

437.8 

24.1 

C02 

82 

42 

3.73 

RR 

10.34 

1.4 

61.5 

424.0 

25.4 

C02 

16 

42 

5.63 

RR 

4.20 

0.6 

317.5 

2189.1 

23.6 

C02 

101 

42 

8.23 

RR 

3.70 

0.5 

59.5 

410.2 

24.6 

HE 

11 

40 

1.49 

SMR 

166.80 

22.2 

25.5 

175.8 

25.1 

C02 

12 

40 

1.53 

SMR 

138.00 

18.4 

23.0 

158.6 

24.8 

C02 

10 

40 

1.72 

CMR 

226.00 

30.1 

49.5 

341.3 

24.9 

C02 

143 

40 

1.90 

CMR 

46.18 

6.2 

12.5 

86.2 

24.6 

C02 

8 

40 

1.94 

DMR 

143.00 

19.1 

49.0 

337.8 

24.6 

C02 

2 

40 

2.07 

DMR 

100.00 

13.3 

47.0 

324.1 

23.7 

C02 

177 

40 

2.08 

DMR 

30.05 

4.0 

14.0 

96.5 

24.0 

C02 

176 

40 

2.15 

DMR 

56.20 

7.5 

27.0 

186.2 

24.2 

C02 

T2-3 


Table  2  (continued) 


Initial  Conditions  for  SFs  Experiments 


EXP  ew  Ms  Type  P0  P3 


Dr iver 
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Table  2  (continued) 


Initial  Conditions  for  SF6  Experiments 


EXP 

6w 

(deg) 

Ms 

Type 

Po 

(torr) 

(kPa) 

(psi) 

Pb 

(kPa) 

To 

(°C) 

Driver 

35 

37 

1.70 

CMR 

91.00 

12.1 

20.6 

142.0 

23.6 

C02 

36 

37 

1.86 

cm 

62.45 

8.3 

20.0 

137.9 

23.5 

C02 

37 

37 

2.07 

CMR 

42.65 

5.7 

19.8 

136.5 

23.3 

C02 

61 

37 

2.13 

DMR 

47.80 

6.4 

45.0 

310.3 

24.3 

C02 

41 

37 

2.18 

DMR 

35.20 

4.7 

19.0 

131.0 

23.7 

C02 

63 

37 

2.23 

DMR 

50.70 

6.8 

45.0 

310.3 

23.4 

HE 

42 

37 

2.31 

DMR 

28.90 

3.9 

19.3 

133.1 

24.0 

C02 

39 

37 

4.07 

- 

4.00 

0.5 

319.0 

2199.4 

23.3 

AIK 

67 

37 

4.14 

- 

17.60 

2.3 

177.5 

1223.8 

24.2 

C02 

68 

37 

4.14 

- 

17.60 

2.3 

179.0 

1234.2 

24.2 

C02 

62 

37 

4.14 

DMR 

17.60 

2.3 

179.0 

1234.2 

24.9 

C02 

69 

37 

4.15 

- 

17.60 

2.3 

180.5 

1244.5 

24.4 

C02 

70 

37 

4.18 

- 

11.95 

1.6 

132.5 

913.6 

24.6 

C02 

59 

37 

4.18 

- 

11.95 

1.6 

132.5 

913.6 

23.3 

C02 

64 

37 

4.20 

DMR 

11.97 

1.6 

133.0 

917.0 

23.7 

C02 

58 

37 

4.21 

- 

11.95 

1.6 

133.5 

920.4 

24.2 

C02 

71 

37 

4.22 

DMR 

11.95 

1.6 

136.0 

937.7 

24.4 

C02 

56 

37 

5.19 

DMR 

4.19 

0.6 

285.5 

1968.5 

23.8 

C02 

73 

37 

5.52 

DMR 

5.42 

0.7 

364.0 

2509.7 

23.7 

C02 

57 

37 

5.55 

DMR 

4.18 

0.6 

277.0 

1909.8 

23.4 

C02 

43 

37 

5.78 

DMR 

5.00 

0.7 

182.5 

1258.3 

24.3 

N2 

T2-5 


.  v 


% 
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Table  2  (continued) 


Initial  Conditions  for  SFe  Experiments 


EXP 

ew 

(deg) 

Ms 

Type 

Po 

(torr) 

(kPa) 

(psi) 

Pa 

(kPa) 

To 

(°C) 

Driver 

171 

37 

5.95 

DMR 

2.82 

0.4 

378.0 

2606.2 

24.1 

C02 

172 

37 

6.48 

DMR 

15.10 

2.0 

94.0 

648.1 

24.3 

HE 

40 

37 

6.73 

- 

3.95 

0.5 

310.5 

2140.8 

23.4 

AIR 

66 

37 

6.76 

- 

4.18 

0.6 

316.0 

2178.7 

25.1 

C02 

53 

37 

7.65 

- 

3.63 

0.5 

59.0 

406.8 

23.5 

HE 

54 

37 

7.66 

DMR 

4.17 

0.6 

288.0 

1985.7 

23.5 

C02 

55 

37 

7.72 

Dl® 

4.17 

0.6 

281.0 

1937.4 

22.7 

C02 

50 

37 

7.75 

- 

3.63 

0.5 

60.8 

419.2 

22.8 

HE 

65 

37 

7.89 

DMR 

4.21 

0.6 

306.0 

2109.8 

24.5 

C02 

72 

37 

7.96 

DMR 

11.05 

1.5 

184.0 

1268.6 

23.4 

HE 

173 

37 

7.99 

DMR 

6.36 

0.8 

92.0 

634.3 

24.8 

HE 

174 

37 

8.11 

- 

6.36 

0.8 

94.0 

648.1 

24.5 

HE 

52 

37 

8.34 

- 

3.00 

0.4 

59.5 

410.2 

23.0 

HE 

51 

37 

8.45 

DMR 

3.00 

0.4 

62.5 

430.9 

23.0 

HE 

38 

37 

8.81 

- 

3.95 

0.5 

316.0 

2178.7 

23.3 

AIR 

124 

30 

1.76 

SMR 

116.70 

15.6 

22.0 

151.7 

23.9 

C02 

126 

30 

1.81 

CMR 

101.20 

13.5 

22.0 

151.7 

23.6 

C02 

128 

30 

1.93 

CMR 

78.00 

10.4 

22.4 

154.4 

23.7 

C02 

127 

30 

2.19 

CMR 

44.10 

5.9 

21.5 

148.2 

23.5 

C02 

125 

30 

2.30 

CMR 

23.80 

3.2 

20.5 

141.3 

23.8 

C02 

142 

30 

2.35 

CMR 

17.95 

2.4 

12.5 

86.2 

22.9 

C02 

140 

30 

2.59 

DMR 

11.17 

1.5 

12.1 

83.4 

24.2 

C02 

T2-6 


Table  2  (continued) 


Initial  Conditions  for  SF6  Experiments 


EXP 

6W 

(deg) 

Ms 

Type 

rr 

O 

*1 

v-/ 

►d 

o 

(kPa) 

(psi) 

Ps 

(kPa) 

To 

(°C) 

Dr  iver 

141 

30 

2.72 

DMR 

8.70 

1.2 

12.0 

82.7 

24.4 

C02 

131 

20 

2.08 

SMR 

30.00 

4.0 

12.5 

86.2 

23.3 

C02 

135 

20 

2.11 

SMR 

29.40 

3.9 

12.1 

83.4 

23.8 

C02 

129 

20 

2.13 

CMR 

50.50 

6.7 

21.5 

148.2 

23.2 

C02 

130 

20 

2.83 

DMR 

13.30 

1.8 

20.6 

142.0 

23.6 

C02 

132 

20 

2.84 

DMR 

13.30 

1.8 

20.5 

141.3 

23.6 

C02 

133 

20 

2.97 

DMR 

10.70 

1.4 

20.7 

142.7 

23.8 

C02 

134 

20 

3.02 

DMR 

9.70 

1.3 

20.2 

139.3 

23.6 

C02 

136 

20 

3.15 

DMR 

7.90 

1.1 

20.2 

139.3 

24.0 

C02 

137 

20 

3.51 

DMR 

4.58 

0.6 

20.6 

142.0 

24.0 

C02 

138 

20 

3.76 

DMR 

7.39 

1.0 

46.5 

320.6 

23.6 

C02 

139 

20 

4.08 

DMR 

4.82 

0.6 

47.5 

327.5 

23.8 

C02 

175 

10 

2.60 

SMR 

11.17 

1.5 

12.1 

83.4 

25.1 

C02 

150 

10 

2.86 

CMR 

13.33 

1.8 

23.0 

158.6 

23.1 

C02 

151 

10 

3.00 

CMR 

9.70 

1.3 

21.6 

148.9 

23.8 

C02 

153 

10 

3.14 

CMR 

7.95 

1.1 

23.1 

159.3 

24.3 

C02 

154 

10 

3.45 

DMR 

11.43 

1.5 

46.5 

320.6 

23.9 

C02 

155 

10 

3.65 

DMR 

12.08 

1.6 

63.0 

434.4 

24.3 

C02 

152 

10 

3.73 

DMR 

10.42 

1.4 

62.5 

430.9 

24.2 

C02 

156 

10 

4.06 

DMR 

6.44 

0.9 

63.0 

434.4 

24.0 

C02 

157 

10 

4.25 

DMR 

4.95 

0.7 

64.5 

444.7 

23.9 

C02 

2- 
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Table  2  (continued) 

Initial  Conditions  for  SF6  Experiments 


EXP 

ew 

(deg) 

Ms 

Type 

Po 

(torr) 

(kPa) 

(psi) 

P  3 

(kPa) 

To 

(°c) 

Driver 

158 

10 

4.43 

nim 

10.18 

1.4 

159.0 

1096.3 

24.6 

C02 

169 

4 

3.37 

SMR 

13.55 

1.8 

49.5 

341.3 

24.1 

C02 

161 

4 

3.47 

CMR 

15.95 

2.1 

67.5 

465.4 

23.5 

C02 

160 

4 

3.74 

CMR 

10.40 

1.4 

64.0 

441.3 

23.5 

C02 

159 

4 

4.05 

CMR 

5.00 

0.7 

49.0 

337.8 

23.2 

C02 

170 

4 

4.28 

CMR 

4.98 

0.7 

67.5 

465.4 

24.1 

C02 

166 

4 

4.36 

CMR 

10.15 

1.4 

151.0 

1041.1 

24.3 

C02 

162 

4 

5.50 

CMR 

5.42 

0.7 

388.0 

2675.2 

23.8 

C02 

165 

4 

5.94 

- 

6.70 

0.9 

383.0 

2640.7 

24.1 

AIR 

164 

4 

5.94 

- 

2.78 

0.4 

386.0 

2661.4 

25.1 

C02 

163 

4 

6.01 

- 

13.52 

1.8 

63.5 

437.8 

24.9 

HE 

167 

4 

6.26 

- 

17.50 

2.3 

95.5 

658.4 

24.6 

HE 

168 

4 

6.65 

- 

15.10 

2.0 

101.0 

696.4 

23.9 

HE 

Comparison  for  MR  in  SFs  at  fixed 


177  40  2  08  DMR  1172  158  30  45  51174  155  30  471  1043  28 


Table  3  (continued) 


Experiments  Equilibrium  Analysis  Frozen  Analysis 


in  m 

> 0  o 


H  H  O 
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o 
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OJ 
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N 
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43 

in 

* 

n 

n 

n 

n 

n 

n 

m 

01 

>0 

n 

n 

CD 

O' 

4) 

CO 
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n 

oi 

OJ 

oi 

iH 

^4 

t4 

H 

w 

,rH 

H 

m 

Oi 

xf 

o 

o 

OJ 

o 

O' 

O' 

oi 

ri 

o i 

oi 

o* 

& 

43 

H 

m 

in 

in 

n 

n 

OJ 

OJ 

OJ 

OJ 

n 

H 

^4 

OI 

OJ 

n 

m 

CO 

N 

CD 

CD 

n 

CD 

4) 
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o 
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O' 

O' 
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01 
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Table  4  (continued) 


Table  4  (continued) 
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Table  5 


Comparison  for  RR  in  SFg 
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Table  5  (continued) 
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Table  6 


Numerical  Results  of  Pressure 
at  RR  t  DMR  Transition  Boundary 

for  Y  •  5/3  at  P0  ■  15  torr  (2  kPa)  ,  T0  -  300  K 


Type 

Region 

Pressure  * 

torr  (kPa) 

A  % 

Stagnation  Pressure 
at  Corner  C 

torr  (kPa) 

RR 

2 

1547.1 

(206.3) 

-5.51 

1283.6 

(171.1) 

DMR 

A, 5 

1461.9 

(194.9) 

RR 

2 

2570.0 

(342.6) 

-0.43 

2595.7 

(346.1) 

DMR 

A, 5 

2558.9 

(341.2) 

RR 

2 

3831.1 

(510.8) 

1.82 

4311.1 

(574.8) 

DMR 

A, 5 

3900.7 

(520.1) 

RR 

2 

5328.5 

(710.4) 

3.03 

6401.5 

(853.5) 

DMR 

A, 5 

5490.1 

(732.0) 

RR 

2 

7058.6 

(941.1) 

3.80 

8850.9 

(1180.0) 

DMR 

A, 5 

7326.9 

(976.8) 

RR 

2 

9018.6 

(1202.4) 

4.35 

11650.0 

(1553.2) 

DMR 

A, 5 

9410.8 

(1254.7) 

RR 

2 

11212.1 

(1494.8) 

4.72 

14793.6 

(1972.3) 

DMR 

A, 5 

11741.1 

(1565.4) 

a)  Regular  Reflection 


b)  Mach  Reflection 
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Fig.  6  Schematic  diagram  of  double-Mach  reflection. 


Fig.  8  Detachment  and  mechanical  equilibrium  criteria 
for  RR(MR)  termination,  y  =  1.4. 


Fig.  9(a)  Schematic  diagram  illustrating  assumption  of 
position  of  K  or  T*  and  flow  direction  in 
region  (1)  relative  to  K. 
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Fig.  10  Basic  instrumentation  setup  for  recording 
shock  wave  reflection  process. 


a)  Present  model 


b)  Elevated  wedge  model 


c)  Symmetric  wedge  model 


d)  Asymmetric  wedge  model 


Fig.  11  Schematic  diagram  of  possible  models  for  studying 
shock  wave  reflection  in  a  shock  tube. 
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Fig.  30(d)  Finite-fringe  interf erogram  of  an  attached 
I'Oinpl ex-Mach  reflection  in  COg  [24]. 

Ms  =  3.43,  0W  =  20.3°,  P0  =  30.0  torr, 

T3  =  294.0  K,  A  =  6943  X. 
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Fig.  31(e)  Transition  lines  for  RR  J  MK  and  CMR  Z  DMR  for 
perfect  gases  with  different  y  and  imperfect 
air  [25]. 
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Fig.  33(a)  Comparison  of  predicted  domains  of  various  types 
of  reflection  and  experimental  results  in  the 
plane  for  vibrational  equilibrium  SFf  . 
SMR  i  CMR  transition  boundaries  : 
solid  line  is  6  =  90°,  dashed  line  is  Mn.  =  1. 
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Comparison  of  predicted  domains  of  various  types 
of  reflection  and  experimental  results  in  the 
(Ms-,',w>  P 1  anc‘  for  vibrational  equilibrium  SF6  . 
SMR  CMR  transition  boundaries  : 


solid  line  is  <S  =  90°,  dashed  line  is  M- 


Fig.  34  Infinite-fringe  interf erogram  of  a  double-Mach 
reflection  with  X  <  1°  in  SF6 . 

Ms  =  1.81,  9W  =  43°,  P0  =  101.2  torr, 

T0  =  297.1  K,  X  =  3471.5  A, 
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Fig.  36  Measured  ratio  £,*/ L*  vs.  wedge  angle  for 

several  Mach  numbers  in  SF6 .  See  Fig.  9(a) 
for  the  definition  of  £*  and  L*. 
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Fig.  50  Variation  of  x'  vs  Mg  at  fixed  0W  for  vibrational 
equilibrium  SF6,  Po  =  15  torr,  To  =  300  K. 
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Fig.  52  Variation  of  x'  vs  Mg  at  fixed  0W  for  frozen  SF6 
with  y  =  4/3,  P0  =  15  torr,  T0  =  300  K. 
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Fig.  53 
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Plot  of  the  measured  values  for  (x'  -  x)  vs  Ms 
for  several  wedge  angles  in  SF6 . 


(b)  x'  =  x,  cu‘=  0  ,  el  --  8„ 


Fig.  54  Schematic  diagrams  of  four  patterns 
of  double-Mach  reflection. 
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Fig.  55  Shadowgram  showing  a  double-Mach  reflection 
with  (jj*  <  0  and  x*  =  0. 

Ms  =  7.96,  0W  =  37°,  Po  =  11.0  torr. 

To  =  296.6  K  in  SFs . 


m 


Fig.  56  Infinite-fringe  interf erogram  showing  a 
double-Mach  reflection  with  u)'  <  0  and 

X '  <  x  • 

Ms  =  4.20,  ew  =  37°,  Po  =  12.0  torr. 

To  =  296.9  K  in  SF6 . 


a)  Regular  Reflection 


b)  Double-Mach 
Reflection 


Fig.  60  Location  of  highest  pressure  in  RR  and  DMR  near 
the  RR  t-  DMR  transition  boundary. 

Point  C  -  stagnation  pressure  at  corner, 

Point  E  -  highest  pressure  in  RR, 

Region  F  -  highest  pressure  in  DMR, 

Point  G  -  pressure  behind  second  triple  point  T 


P2  in  RR  (point  E) 
and  P5  in  DMR  (point  G) 


Fig.  61  Pressure  behaviour  at  point  E  in  RR  and  at  point  G  in  DMR 
as  a  function  of  wedge  angle.  Ms  =  6.0,  Po  =  15  torr 
(2  kPa),  and  To  =  300  K  in  argon. 
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Method  of  Calculations 


Reflections  Inside  a  Shock  Tube 


Since  the  planar  incident  shock  wave  moves  with  constant  velocity  along 
the  wedge  surface,  the  entire  reflection  phenomenon  can  be  considered  as 
pseudo-stationary  in  a  frame  of  reference  attached  to  a  point  which  moves 
with  a  constant  velocity.  Therefore,  instead  of  three  independent  variables 
x,  y  and  t,  the  phenomenon  is  described  in  terms  of  x/t  and  y/t  and  the  flow 
is  self-similar  [19].  In  a  pseudo-stationary  flow,  the  shock  wave 
configuration  remains  self-similar  and  grows  linearly  with  time  from  the 
moment  the  incident  shock  wave  collides  with  the  wedge  corner.  The 
phenomenon  is  also  assumed  to  be  a  two-dimensional  inviscid  flow. 

The  calculations  for  obtaining  the  flow  properties  of  each  region  in 
the  angular  vicinity  of  a  shock  confluence  point  require  the  formulation  of 
the  oblique-shock-wave  relations  for  each  shock  wave.  The  reference  point 
is  the  reflection  point  P  in  RR,  and  the  triple  point  T  in  MR;  the  region 
which  are  divided  by  shock  waves  and  a  slipstream  are  designated  by  0  to  3 
as  shown  in  Fig.  3a  and  4.  The  details  of  the  analytical  formulations  can 
be  found  in  Law  [47].  The  physical  quantities  on  both  sides  of  each  shock 
wave  in  the  vicinity  of  the  reference  point  satisfy  the  following  equations 
of  motion: 

Conservation  of  tangential  velocity 

p^tan<t>j  -  pj tan( 4>-j  -  0j)  (A.l) 

Continuity 

PjU^si n<t>j  =  pjUjSinUf  -  9j)  (A. 2) 

Normal  momentum 

Pj  +  2sin2<t>1-  =  Pj  +  pjUj2sin2(<t'1*  -  Gj)  (A. 3) 

Energy 

h^  +  —  Ui2sin2d>i  =  hi  +  -  U.:  2s  in2(  <K-  -  Gj  (A. 4) 

i  2  *  i  j  2  J  *  j 

where  i  and  j  refer  to  the  states  ahead  and  behind  the  shock,  respectively. 
The  equations  of  state  applied  are 

P  =  p(P,  T),  h  =  h(P,  T)  (A. 5) 

In  RR,  or  the  two-shock  theory,  there  are  two  sets  of  equations: 

For  incident  shock  I:  i  =  0,  j  =  1  (A. 6) 

For  reflected  shock  R:  i  =  1,  j  =  2  (A.7) 


These  two  sets  of  equations  are  solved  under  the  boundary  condition 
that  the  flow  deflections  across  the  incident  and  reflected  shocks  be  equal 
and  opposite  with  the  direction  of  the  flow  behind  the  reflected  shock  wave 
being  along  the  wedge  surface,  i.e., 

ei  =  e2  (A. 8) 

The  calculating  procedures  for  the  solution  of  equation  (A. 5)  through  (A. 8) 
are  as  follows  [11]: 

1.  The  initial  conditions  of  the  pressure  PQ,  the  temperature  T0,  the 

Mach  number  Ms  and  the  wedge  angle  9W  are  given. 

2.  The  first  set  of  equations  (A. 6)  is  solved  for  P,,  T,,  U,,  and  8,. 

Note  that  U0  =  Us  sec9w  and  <t>0  =  90°  -  9W. 

3.  Some  initial  guess  for  <t>1  is  made. 

4.  The  second  set  of  equations  (A. 7)  is  solved  to  obtain  the  values  for 

P2,  T2,  U2  and  02. 

5.  The  calculated  results  are  checked  against  the  boundary  condition 

( A.8) . 

6.  A  new  value  of  <t>1  is  predicted. 

7.  The  procedure  is  repeated  until 

I  Si  “  02l  <  10'4 

In  the  case  of  MR,  or  the  three-shock  theory,  there  are  three  sets  of 
equations  for  the  first  triple-point  system: 

For  incident  shock  I:  i  =  0,  j  =  1  (A.9) 

For  reflected  shock  R:  i  =  1,  j  =  2  (A. 10) 

For  Mach  stem  M  :  i  =  0,  j  =  3  ,  ( <t> -j  is  referred  to  as  <i>3)  (A. 11) 

The  above  three  sets  of  equations  are  solved  under  the  boundary 
conditions  that  the  flow  deflection  across  the  incident  and  reflected  shocks 
be  equal  to  the  flow  deflection  across  the  Mach  stem  and  the  pressure  on 
both  sides  of  the  slipstream  be  identical,  i.e., 

e3  =  b1  -  e2  (A. 12) 

P2  =  P3  (A. 13) 

In  MR,  it  is  necessary  to  predict  a  value  of  the  triple-point 
trajectory  angle  x  in  order  that  the  system  of  equations  be  closed.  The 
Mach  stem  is  assumed  normal  to  the  wedge  surface  based  on  experimental 
observations,  and  Law  and  Glass  [15]  introduced  an  additional  geometrical 
rel at  ion 


w 


<t>3  =  90°  -  x  (A. 14) 

The  calculating  procedures  for  the  solution  of  equations  (A. 5)  and 

(A. 9)  through  (A. 14)  are  as  follows  [11]  : 

1.  The  initial  conditions  of  P0,  T0,  Ms  and  0W  are  given. 

2.  Some  initial  guesses  for  x  and  4^  are  made. 

3.  The  first  set  of  equations  (A. 9)  is  solved  for  P,,  T,,  U,  and  0,.  Note 

that  U0  =  Us  cosec<}>0  and  4>0=  90°  -  (6#  +  x)* 

4.  The  second  set  of  equations  (A. 10)  is  solved  for  P2,  T2>  U2  and  02. 

5.  Using  Eq.  (A. 14),  the  third  set  of  equations  (A. 11)  is  solved  for  P3, 

T3,  U3  and  03. 

6.  The  calculated  results  are  checked  against  the  boundary  conditions 

(A. 12)  and  (A. 13). 

7.  New  values  of  x  and  4>:  are  predicted. 

8.  The  solutions  are  iterated  until 


nra 


|03  -  el  +  02|  <  io-4 
|(P3  -  p2)/p2|  <  io-4 


Details  of  the  analytical  formulation  of  the  second  triple  point  T'  is 
given  by  Ben-Dor  Lll].  The  flow  fields  around  T'  are  shown  in  Fig.  59  with 
the  prime  to  denote  that  the  properties  are  measured  with  respect  to  T‘. 
The  three  set  of  equations  for  the  system  are  : 


For  reflected  shock  R: 


i  =  1,  j  =  2 


(A. 15) 


For  second  reflected  shock  R': 


i  =  1,  j  =  4 

(4>-j  is  referred  to  as  4>4) 


(A. 16) 


For  second  Mach  stem  M* 


i  =  2,  j  =  5 


(A. 17) 


The  above  three  sets  of  equations  are  solved  under  the  boundary 
conditions  that  the  flow  deflection  across  the  reflected  shock  R  and  the 
second  Mach  stem  M'  be  equal  to  the  flow  deflection  across  the  second 
reflected  shock  R‘,  and  the  pressure  on  both  sides  of  the  second  slipstream 
S'  be  identical ,  i .e. , 


04  "  02  "  05 


(A. 18) 


(A. 19) 


The  calculation  involves  a  transformation  of  the  already  calculated 

values  of  Ui  and  4>i  from  a  frame  of  reference  attached  to  T  to  a  frame  of 
reference  attached  to  T‘.  Since  thermodynamic  properties  do  not  depend  on 
the  frame  of  reference,  the  prime  is  omitted  from  the  thermodynamic 
variables  in  the  equations.  The  values  of  U'  and  4>!  are  given  as 


-  A. 3  - 


**  M  ‘j*  ’  J*  "  « 


- u*  D 


sin‘!(<i)1  +  4>0  -  ©i)  sin(  <t>!  +  4>0  -  eL 


*i  ~  ^  1  - 


(A. 21) 


The  calculating  procedures  for  the  solution  of  equations  (A. 5)  and 
(A. 15)  through  (A. 21)  are  as  follows: 

1.  The  initial  conditions  of  Plt  Tp  M^,  <f0,  4^  and  01  are  given. 

i  t 

2.  Some  initial  guesses  for  4>2  and  4>4  are  made. 

3.  Using  equations  (A.2(j)  and  (A. 21),  the  first  set  of  equations  (A. 15)  is 
solved  for  P2,  T2,  U2  and  02. 

I  I 

4.  The  second  set  of  equations  (A. 16)  is  solved  for  P4,  T4,  U4  and  04. 

I  i 

5.  The  third  set  of  equations  (A. 17)  is  solved  for  P5,  T5,  U5  and  ©5. 

6.  The  calculated  values  are  checked  against  the  boundary  conditions 
(A. 18)  and  (A. 19). 

i  i 

7.  New  values  of  <t>2  and  <t>4  are  predicted. 

8.  The  solutions  are  iterated  until 

|<t,4  -  02  +  0j|  <  10-4 

and 

|(P„  -  P5)/P5|  <  10-" 


Ground  Reflection  of  a  Blast  Wave 

The  curve  fit  that  was  reported  by  Sadek  and  Gottlieb  [66]  to  represent 
the  shock  trajectory  for  a  1  kg  TNT  ground  explosion  in  a  standard 
atmosphere  is 


a,  a0t 


rS  =  iVaV  +  aits  +  a3^(1  +  “^s)]1'1 


(A. 22) 


where  rs  is  the  shock  radius,  ts  is  the  shock  arrival  time,  a,  =  0.72 
m,  a2  =  1.872  m/s,  a3  =  1.29  m,  a4  =  0.3403  m/s  and  a^  =  0.3403  m/ms. 
Mach  number  of  the  blast  front  can  be  obtained  by  differentiating  the  above 
equation  and  then  dividing  it  by  the  speed  of  sound. 

In  the  calculation  of  the  blast-front-trajectory  path  on  the  (Ms  - 
0W)  transition  map,  an  incident  shock  wave  Mach  number  Ms  is  first 
chosen.  The  time  of  arrival  ts  is  then  determined  and  substituted  into 
equation  (A. 22)  to  obtain  the  shock  radius  rs.  To  obtain  an  equivalent 
shock  radius  r  for  the  case  of  HUB  explosion  above  a  perfect  planar  surface, 
r<-  is  multiplied  by  the  factor  21/3. 


For  RR,  at  a  given  HOB  h  and  Ms,  the  radius  of  the  blast  front  r  can 
be  calculated  as  mentioned  above.  Ground  range  x  and  the  corresponding 
wedge  angle  9W,  as  shown  in  Fig.  A-l(a),  can  be  determined  from  the 
relations 


c2  =  r2  -  h2 


sin0w  =  h/r 


(A. 23) 
(A. 24) 


For  MR,  at  a  given  HOB  h  and  Ms,  the  radius  of  the  blast  front  r  is 
calculated  the  same  way  as  for  RR.  To  find  the  ground  range  x  and  the  wedge 
angle  9W,  it  is  essential  to  know  the  height  of  the  first  triple  point 
Yj.  As  illustrated  in  Fig.  A-l(b),  x,  h,  Yj  and  ow  are  related  by 

r2=  x2  +  ( h  -  Yj) 2  (A. 25) 

tanew  =  (h  -  Yj)/x  (A. 26) 

and  the  slope  of  the  first-triple-point  trajectory  is  given  by 

=  tany  (A. 27) 

dx 

If  equation  (A. 26)  is  first  differentiated,  then  substituted  into 
equation  (A. 27)  and  finally  the  result  is  integrated,  the  following  two 
relations  are  obtained  [67]: 


Yt  .Ww  tanx  dew 

_L  =  1  -  exp  /  - - - 

h  ed  sin2ew(l  +  tanx/tanew) 


for  ew  >  ed 


for  ew  <  ec 


(A. 28) 


(A. 29) 


where  9d  is  the  RR«->MR  transition  angle  according  to  the  detachment 
criterion.  Equation  (A. 28)  is  applicable  to  -RR  which  states  that  no  triple 
point  exists  if  9W  is  above  9q.  Equation  (A. 29)  is  applicable  to  MR  and 
can  be  integrated  if  the  relation  between  9W  and  j;  in  a  two-dimensional 
shock  tube  flow  is  assumed  valid  here. 

The  following  fix  point  iterative  scheme  is  used  to  determine  the  value 
of  Yj  : 

1.  The  HOB  h  and  the  incident  shock  Mach  number  Ms  are  given. 

2.  Calculate  the  shock  radius  r  with  the  given  Ms. 

3.  Some  initial  guess  of  Yj  is  made. 

4.  Calculate  wedge  angle  9W  from  the  relation  obtained  by  substituting 

equation  (A. 25)  into  (A. 26) 


9w  ■ 


(h  -  Yt)2 


-  if1'2! 


(A. 30) 


A  -  A  .v 


5.  With  9W  calculate  Yj  using  equation  (A. 29). 

6.  Check  if  the  difference  of  the  new  Yj  and  the  last  Yj  is  within  the 
tolerance  of  10“4. 

7.  Steps  4  and  5  are  repeated  until  the  tolerance  in  step  6  is  achieved. 

Once  Yt  is  determined,  x  and  ew  can  be  calculated  using  equation 
(A. 25)  and  (A. 26).  In  real  explosions,  the  Mach  stem  M  is  not  entirely 
straight  from  the  first  triple  point  to  the  ground,  but  curved  as  shown  in 
Fig.  A-l(b).  Such  bulging  effect  of  the  Mach  stem  is  simply  simulated  by 
assuming  that  the  Mach  stem  M  is  an  arc  cutting  the  first  triple  point  at  T 
and  the  ground  at  x'  with  its  radius  being  r^.  The  ground  range  x  in  MR 
then  takes  the  value  of  x'. 


J 


Triple-point 

Trajectory 


Thermodynamic  Properties  of  Sulfur  Hexafluoride 


Sulfur  hexafluoride  is  a  gaseous  dielectric  of  great  importance  in  the 
electrical  and  communications  industries  as  an  insulating  and  quenching 
medium  in  circuit  breakers  [48].  It  alone  or  in  a  mixture  with  other  gases 
has  the  possibilities  as  a  refrigerant.  Its  sysmmetric  shape  and  the 

presence  of  the  six  fluorine  atoms  make  it  of  theoretical  interest  in  the 
study  of  intermol ecul ar  forces.  Plasma  environment  studies  of  a  body  during 
atmospheric  reentry  have  created  an  interest  in  the  effect  of  electrophilic 
gases  on  high  temperature  plasmas.  In  particular,  when  SF6  is  injected  into 
a  supersonic  plasma  stream,  a  large  eletron  density  reduction  is  observed 
[49].  Recently,  in  the  course  of  finding  a  method  of  simulating 

blast-intercept  loadings  on  full  scale  reentry  vehicle  models,  the  heavy  gas 
SF6  was  used  since  it  would  produce  high  flow  Mach  numbers  and  high  dynamic 
pressures  [50]. 

The  polyatomic  SF&  molecule  has  one  sulfur  atom  surrounded  by  six 
fluorine  atoms  in  the  form  of  a  regular  octahedron  with  15  vibrational 
degrees  of  freedom  (see  Appendix  C) .  The  characteristic  modes  of  vibration 
of  SF6  with  corresponding  degeneracy  factors  and  wave  numbers  [51]  are 
shown  in  Fig.  Bl. 

The  thermodynamic  quantities  are  given  by 

P  =  p  £  T,  a  *  /yP/p,  U  =  Ma 
m 

In  the  frozen-gas  model,  only  translational  and  rotational  degrees  of 
freedom  are  excited  to  the  new  equilibrium  state  and  the  other  degrees  are 
frozen  at  the  initial  state,  thus 


h  = 


m 


where 

Y0  -  4/3 
and 

y  =  y0 

In  the  equilibrium-gas  model,  all  vibrational,  translational  and  rotational 
modes  are  excited  to  the  new  equilibrium  state  immediately  behind  the  shock 
wave.  The  enthalpy  and  specific  heats  ratio  are  given  by 

h  =  — Y°—  i  T  O'  nk  _ I* _ i 

y0  -  1  m  k  exp(Tk/T-l)  m 

y/(y0-1)  +  V 
1/(y0-1)  +  V 


Y 


where 


Tl/T 


V  =  y  nkexp(Tk/T)  f. 

k  exp(Tk/T)  -  1 

and 


Tk  =  hvk/K 

with  h  being  the  Planck  constant  and  K  the  Boltzmann  constant. 

The  composition  of  SF6  has  beeen  computed  by  Frie  [52]  at  atmospheric 
pressure,  including  the  chemical  species  of  SF6,  SF4,  SF2,  S2,  F2,  S,  F,  S+, 
S++,  F+,  F"  and  e“.  The  variation  of  particle  density  with  temperature  for 
SF6  is  shown  in  Fig.  B2.  Minute  amounts  of  SF4  and  F  begin  to  occur  around 
a  temperature  of  1000  K  and  the  SF6  line  drops  from  the  total  number  of 
particles  (ini)  line  at  a  temperature  about  1400  K,  with  the  dissociation 
reaching  a  peak  slightly  below  2000  K.  The  ionization  of  sulfur  and 
fluorine  peaks  at  12000  K  and  17000  K,  respectively.  Frost  and  Liebermann 
[48]  have  extended  the  computation  of  the  composition  of  SF6  up  to  a 
pressure  of  16  atm  and  a  temperature  of  45000  K.  Their  results  agree  well 
with  the  earlier  calculation  of  Frie  [52]  at  1  atm.  The  compressibility 
factor  Z,  which  is  an  indicator  of  the  number  of  moles  of  gas  formed  due  to 
dissociation  and/or  ionization  from  one  mole  of  gas  at  room  temperature,  of 
SF6  at  various  pressures  and  temperatures  is  reproduced  from  reference  [48] 
in  Table  Bl.  The  value  of  Z  in  the  temperature  range  from  1000  K  to  2300  K 
for  several  pressures  is  plotted  in  Fig.  B3.  It  is  seen  in  this  figure  that 
the  dissociation  process  is  intensified  with  increasing  temperature  and 
decreasing  pressure  as  expected. 


Using  Table  Bl  and  Fig.  B3,  the  range  of  Mach  number  in  which  no 
dissociation  of  SF6  occurs  can  be  determined.  In  the  discussions  which 
follow,  larger  than  5%  of  dissociation  of  the  gas  is  considered  significant. 
According  to  the  vibrational  equilibrium-gas  analysis  at  Ms  =  11  with  an 
initial  pressure  and  temperature  P„  =  15  torr  (2  kPa)  and  T0  =  300  K,  the 
pressure  P:  =  2.53  atm,  the  temperature  T1  =  1470  K  with  Z  *  1.055  in  region 
(1)  behind  the  incident  shock  front.  Thus  the  vibrational  equilibrium-gas 
analysis  predicts  that  dissociation  will  become  significant  at  Ms  >  11. 

Since  the  comparison  of  the  analytical  and  experimental  results  is  done  for 
cases  with  Ms  <  5.5,  it  is  necessary  to  check  if  dissociation  occurs 
behind  the  reflected  shock  wave  as  well.  As  discussed  in  Sec.  4.6,  for  a 
given  Ms,  the  highest  pressure  obtained  from  reflection  occurs  at  the  RR 
«->  DMR  transition  boundary.  For  Me  =  5.5,  just  before  RR  terminates,  P2  = 
5.11  atm  and  T2  =  790  K  in  region  (2).  However,  this  T2  is  not  high  enough 
to  cause  significant  dissociation  to  occur.  Hence,  in  the  range  of  Mach 
number  Ms  <  5.5,  no  dissociation  of  the  gas  comes  into  play  and  only 
vibrational  excitation  needs  to  be  taken  into  account  in  the  analysis. 
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B3  Compressibility  factor  Z  for  SF6  as  a  function 
of  temperature  at  different  pressures  [48]. 


Table  B1 


Compressibility  Factor  of  Sulfur  Hexafluoride  [48] 


T 

( 103  K) 

1  atm 

i  . . 

2  atm 

Z 

4  atm 

8  atm 

16  atm 

1.0 

1  1.001 

1 

1.0001 

1.0000 

1.0000 

1.0000 

1.4 

1.031 

1.020 

1.0133 

1 . 0088 

1.0059 

1.7 

1.669 

1.4070 

1.2490 

1.1550 

1.098 

2.0 

3.005 

2.9050 

2.6920 

2.3040 

1.8660 

2.3 

3.668 

3.3340 

3.1500 

3.0410 

2.9400 

2.6 

5.490 

4.8100 

4.1500 

3.6280 

3.3010 

3.0 

6.480 

6.3000 

5.9900 

5.4810 

4.8140 

3.5 

6.830 

6.7360 

6.6200 

6.4720 

6.2530 

4 

6.960 

6.9220 

6.8640 

6.7790 

6.6650 

5 

6.996 

6.9920 

6.9850 

6.9700 

6.9430 

6 

7.004 

7.0020 

6.9990 

6.9950 

6.9890 

7 

7.019 

7.0130 

7.0090 

7.0050 

7.0005 

8 

7.064 

7.0450 

7.0320 

7.0230 

7.0150 

9 

7.164 

7.1190 

7.0850 

7.0610 

7.0430 

10 

7.335 

7.2500 

7.1840 

7.1340 

7.0960 

11 

7.555 

7.4370 

7.3340 

7.2500 

7.1830 

12 

7.780 

7.6500 

7.5200 

7.4040 

7.3060 

13 

8.020 

7.8700 

7.7200 

7.5850 

7.4590 

14 

8.380 

8.1400 

7.9500 

7.78400 

7.6340 

15 

1 

8.950 

8.5500 

8.2550 

8.0250 

7.8340 

16 

9,750 

9.1600 

8.7000 

8.3500 

8.0860 

18 

11.710 

10. 8b00 

10.0700 

9.4000 

8.8660 

Table  B1  (continued) 


Compressibility  Factor  of  Sulfur  Hexafluoride  [48] 


T 

(103  K.) 

1  atm 

2  atm 

Z 

4  atm 

8  atm 

16  atm 

20 

13.150 

12.5000 

11.7100 

10.8600 

10.0650 

22 

13.940 

13.5100 

12.9500 

12.2400 

11.4200 

26 

14.770 

14.5300 

14.2200 

13.8200 

13.3100 

30 

15.780 

15.3000 

14.9400 

14.6100 

14.2600 

35 

19.010 

17.8600 

16.8000 

15.9500 

15.3000 

40 

21.230 

20.5700 

19.6600 

18.5500 

17.4000 

v  v.  ■  v .  -I v-" v  ({cv's! 
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Vibrational  Relaxation  in  Sulfur  Hexafluoride 


When  a  shock  wave  propagates  through  a  gas,  the  translational  and 
rotational  degrees  of  freedom  of  the  gas  molecules  are  excited  to  the  new 
equilibrium  state  in  a  few  mean-free-path  lengths.  The  other  internal 
degrees  of  freedom,  on  the  other  hand,  require  a  much  longer  length  to  reach 
their  equilibria.  Relaxation  lengths  play  an  important  role  in  the  analysis 
of  gasdynamic  phenomena,  In  the  case  that  the  internal  degree  of  freedom 
whose  relaxation  length  is  much  longer  than  a  characteristic  length  of  the 
phenomenon,  the  internal  degrees  of  freedom  can  be  treated  as  frozen  at  the 
initial  state.  However,  if  the  relaxation  length  of  the  internal  degrees 
of  freedom  is  shorter  than  a  characteristic  length  of  the  phenomenon,  it  can 
be  assumed  to  be  in  equilibrium  behind  the  shock  wave.  If  the  relaxation 
length  and  the  characteristic  length  lie  between  the  extreme  cases  of  frozen 
and  equilibrium  flow,  the  gas  is  in  nonequi  1  ibrium.  In  the  problem  of 
shock-wave  reflection,  a  character! st  ic  length  of  1  mm  is  chosen  [16]. 

Vibrational  relaxation  of  the  polyatomic  molecule  SFfe  has  been 
extensively  studied  by  ultrasonic  methods  [53,  54,  55]  in  the  temperature 
range  of  200  K  to  400  K.  Vibrational  relaxation  times  were  determined  by 
time-resolved  measurements  of  the  post-shock  density  gradient  using  a  laser 
beam  deflection  technique  over  the  temperature  range  of  450  K  to  1050  K  by 
Breshears  and  Blair  [56].  The  vibrational  relaxation  time  of  SFg 
attributable  to  SF6  -  SFb  collisions  is  expressed  as 

Px  =  2.92  x  10" 3  exp( 38.0  T~ 1 / 3) usec-atm 

where  P  is  the  pressure  in  atm  and  T  is  the  temperature  in  degree  Kelvin. 
Breshears  and  Blair  also  indicated  that  the  ultrasonic  data  in  the  low 
temperature  range  of  200  K  to  400  K  are  in  good  agreement  with  their 
relaxation  time  equation. 

Using  flow  properties  behind  a  shock  wave  and  the  vibrational 
relaxation  time  equation,  a  vibrational  relaxation  length  can  be  obtained. 
The  vibrational  relaxation  lengths  behind  a  normal  shock  wave  versus  Mach 
number  are  plotted  in  Fig.  Cl  for  the  diatomic  gases  N2  and  0 2,  triatomic 
gas  C0?  and  polyatomic  gas  SFb.  These  curves  are  for  the  initial  conditions 
of  P„  =  15  torr  (2  kPa)  and  T0  =  300  K,  and  were  obtained  on  the  basis  of 
the  theoretical  and  empirical  analyses  [56,  57,  58,  59].  It  is  clear  from 
Fig.  Cl  that  for  a  given  Ms ,  N2  has  the  longest  vibrational  relaxation 
length,  which  is  then  followed  by  0,.  CO^  and  SF&.  Among  the  four  gases, 
SF6  has  the  shortest  vibrational  relaxation  length.  This  ordering  can  be 
explained  in  terms  of  the  type  of  molecule  (no.  of  vibrational  degree  of 
freedom)  and  the  vibrational  character i st ic  temperature.  To  describe  the 
motion  of  N  nuclei  in  a  molecule,  it  needs  3N  coordinates;  there  are  3N 
degrees  of  freedom.  For  a  diatomic  and  a  linear  triatomic  molecule,  there 
are  three  translational  and  two  rotational  degrees  of  freedom.  The  number 
of  vibrational  degrees  of  freedom  becomes  3N-5.  For  a  polyatomic  molecule, 
there  are  three  rotational  degrees  of  freedom  instead  of  two,  and  thus  there 
are  3N-6  vibrational  degree  of  freedom.  Consequently,  diatomic  gases  N2  and 
02  have  each  one  degree,  linear  triatomic  gas  C02  has  four  degrees  and 
polyatomic  SFh  has  15  degrees  of  freedom  in  vibration.  The  vibrational 


characteristic  temperatures  T|<  of  these  four  gases  are  listed  also  in 
Fig.  Cl.  N2  has  a  much  higher  than  02,  therefore  it  needs  a  very  high 
temperature  to  get  its  vibrational  mode  excited.  02  can  be  excited  at  a 
lower  temperature,  and  it  reaches  equilibrium  sooner  than  N2,  thus  for  a 
given  Ms,  its  vibrational  relaxation  length  is  shorter  than  N2.  C02  does 

not  only  have  more  vibrational  modes  than  02,  but  also  three  of  its  four 
Tjj's  are  much  smaller  than  that  of  02.  This  implies  that  the  first  three 
vibrational  modes  of  C02  are  excited  earlier,  reaching  equilibrium  sooner 
and  have  a  shorter  relaxation  length  than  02.  Moreover,  the  last  mode  has 

a  Tj<  that  is  higher  than  that  of  N2.  It  is  expected  that  vibrational 

excitation  in  C02  will  reach  complete  equilibrium  at  a  temperature  higher 
than  the  other  four  gases. 

The  average  T^  for  SF6  is  lower  than  for  the  other  four  gases  and  it 
reaches  vibrational  equilibrium  much  faster  than  the  others.  Raising  the 
pressure  at  a  given  temperature  enhances  particle  collisions  and  vibrational 
equilibrium  is  reached  more  quickly.  From  Fig.  A4,  it  is  expected  that 
vibrational  relaxation  in  SF6  at  Ms  <  1.7  with  an  initial  pressure  of  15 
torr  (2  kPa)  and  a  temperature  of  300  K  will  take  place  in  less  than  1  mm. 
However,  at  such  a  low  Mach  number,  a  relatively  higher  initial  pressure  is 
required  for  the  experiment.  Increasing  the  pressure  will  shift  the 

vibrational  relaxation  length  curve  down,  therefore  no  vibrational 

relaxation  length  will  be  observed.  Since  the  vibrational  relaxation  length 
is  less  than  1  mm  over  the  entire  range  of  Mach  number  studied  in  this  work, 
it  can  be  safely  assumed  that  SF6  is  in  vibrational  equilibrium. 


